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ABSTRACT 


We  propose  to  study  the  use  of  the  meteor  burst  channel  in  communication 
systems  by  investigating  the  following  important  aspects: 

1)  The  development  and  enhancement  of  an  accurate  and  reliable  channel 
model  based  on  recently  available  empririal  data.  Analysis  of  this  model  results 
in  analytical  expressions  for  communication  parameters  such  as  channel  duration 
and  throughput  to  be  used  as  design  and  analysis  tools. 

2)  The  optimization  of  throughput  for  fixed  transmission  rate  under  the 
constraint  of  a  given  maximum  bit  error  rate.  This  will  demonstrate  the  room 
for  improvement  in  existing  systems  using  constant  transmission  rate. 

3)  The  feasibility  of  efficiently  communicating  over  the  M3C  using  variable 
bit  rate  and  employing  a  feedback  protocol  to  monitor  the  channel.  This  approach 
will  dramatically  improve  the  throughput  in  comparison  with  constant  transmission 
rate  systems. 

4)  Analysis  of  Automatic-Repeat-Request  (ARQ)  Transmission  over  MBC  by 
studying  performance  measures  such  as  duty  cycle,  throughput  and  waiting  time 
as  a  function  of  packet  length,  coding,  data  rate  and  modulation  technique. 
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1  STATEMENT  OF  THE  PROBLEM 

1.1  Introduction 

There  is  currently  an  urgent  need  for  alternate  channels  for  communication. 
The  large  variety  and  quantity  of  users  throughout  the  world  have  caused  all  of 
the  existing  channels  to  become  extremely  congested.  A  world  wide  search  for 
such  alternate  channels  is  now  being  conducted  for  both  present  day  and  future 
communication  systems  through  the  use  of  different  types  of  channels  such  as 
optical,  laser,  cable,  satellite,  etc.  [14,  15,  17] 

The  particular  application  of  Beyond-Line-of-Sight  (BLOS)  communications , 
which  uses  the  High  Frequency  (HF)  spectral  band  from  3  to  30  MHz,  is  very 
important  to  many  current  and  potential  users.  However,  HF  is  very  sensitive  to 
solar  disturbances  such  as  sunspot  activity,  solar  storms  and  other  galactic 
phenomena,  as  well  as  multipath  returns  from  both  ground  and  atmosphere,  weather 
related  attenuation,  and  other  degrading  factors  [18,19]  .  Because  of  these  factors 
a  channel  for  communication  is  needed  which  has  low  congestion,  is  robust  and 
relatively  indestructible,  has  very  little  outside  interference  and  does  not 
compete  for  bandwidth  with  existing  communication  systems.  The  research  proposed 
herein  represents  such  an  alternative  mode  for  communication. 

Meteors  are  small  aggregates  of  matter  which  upon  entering  the  Earth’s 
atmosphere  burn  up  and  form  long  columns  or  trails  of  ionized  particles  in  the 
upper  atmosphere  at  altitudes  of  80  to  120  km.  The  ionized  columns  diffuse  within 
seconds  yet  are  able  to  support  radio  communication  over  a  range  of  distance  up 
to  2000  km  by  reflecting  a  radio  signal  and  thus  provide  us  with  the  so-called 
meteor  burst  channel  (MBC) .  The  frequency  range  used  is  from  30  to  100  MHz.  The 
lower  limit  is  dictated  by  the  need  to  be  above  the  HF  range  and  avoid 
ionospherically  reflected  signals.  The  upper  limit  is  set  primarily  by  equipment 
sensitivity  limitations  since  reflections  at  the  higher  frequencies  are  weaker 
than  those  at  lower  frequencies. 

Excluding  showers  (annual  Pleiades  shower)  the  occurrence  of  a  meteor  burst 
appears  in  time  to  be  random  with  a  Poisson  arrival  and  an  average  interval  time 
on  the  o. der  of  several  seconds  to  minutes  depending  upon  time  of  day,  season, 
and  global  location,  and  communication  system  parameters. 

Two  types  of  trails  are  described  in  the  literature:  Underdense  trails 
which  are  more  frequent  and  have  an  electron  line  density  below  10l*  e/m;  and 
overdense  trails  which  are  less  frequent  and  have  electron  line  densities  of 
more  than  101*  e/m.  The  power  received  from  an  underdense  trail  is  approximated 
as  an  exponential  time  function  whose  initial  amplitude  is  proportional  to  the 
square  of  the  electron  line  density  and  the  exponential  decay  time  constant 


varies  randomly  from  one  trail  to  the  next.  The  power  received  from  an  overdense 
trail  is  modeled  by  a  complex  time  function  whose  shape  for  a  given  set  of 
communication  system  parameters  is  determined  by  the  value  of  the  electron  line 
density. 

In  the  late  fifties  and  early  sixties,  a  considerable  amount  of  research 
was  conducted  on  the  use  of  meteor  burst  communications  for  beyond  line-of-sight 
(BLOS)  transmission  of  digital  data.  In  recent  years,  interest  in  meteor  burst 
has  been  renewed  because  of  several  factors.  The  development  of  microprocessors 
(enabling  inexpensive  implementation  of  sophisticated  system  control  procedures) 
and  inexpensive  solid-state  memories  (for  the  buffering  required  to  interface 
a  burst  transmission  scheme  to  constant  data  sources  and  sinks)  is  partly 
responsible.  In  addition,  the  nuclear  survivability  of  the  meteor  burst  medium 
is  superior  to  other  BLOS  media  such  as  satellite  and  HF  skywave ,  resulting  in 
considerable  interest  in  the  military  community. 

1.2  Identification  of  the  Problem 


The  research  proposed  herein  represents  such  an  alternative  mode  for  BLOS 
communication.  It  relates  a  novel  app voach  for  BLOS  communication.  The  approach 
is  aimed  at  efficaciously  using  a  previously  under  utilized  and  inefficiently 
used  method  of  communication,  namely,  the  Meteor  Burst  Channel  (MBC)  [1,  9], 
This  medium  possesses  all  of  the  qualities  mentioned  above  for  a  new  communications 
channel.  The  channel  operates  in  the  relatively  unused  lower  portion  of  the 
Very  High  Frequency  (VHF)  band  ranging  from  30  to  100  MHz.  It  avoids  the 
degradations  exhibited  by  HF  on  the  low  end  and  ionospheric  phenomena  at  the 
high  end.  This  region  of  the  spectrum  provides  a  sufficiently  large  bandwidth 
for  efficient  data  communication.  It  is  not  easily  destroyed.  It  has  a  privacy 
feature  inherent  in  its  structure  [1,9]  and  by  virtue  of  its  "newness"  does  not 
interfere  with  existing  systems.  For  MBC,  the  antennas  have  a  high  elevation 
angle  and  rarely  interfere  with  present  use  of  the  frequency  band  when  such  is 
present . 

This  study  seeks  to  define  systematically  and  analytically  the  limits  of 
transmission  of  digital  information  over  the  meteor  burst  channel  by  using  an 
improved  model  to  develop  the  mathematical  relationships  among  link  geometry, 
communication  circuitry  parameters,  physical  meteor  characteristics,  modulation 
technique  and  a  given  communication  protocol.  Furthermore,  the  research  will 
be  aimed  at  the  determination  of  optimal  average  throughput  for  both  constant 
and  variable  transmission  rates  under  the  employment  of  an  information  feedback 
scheme  for  reliable  channel  monitoring. 
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These  results  will  serve  as  a  vital  design  and  analysis  tool  for  the 
communication  system  engineer  using  the  MBC.  Furthermore,  they  will  demonstrate 
the  viability  and  potential  improvement  over  past  and  currently  implemented 
systems.  Finally,  they  will  identify  the  necessary  and  challenging  areas  for 
future  research. 

1 . 3  Significance  of  the  Problem 

Meteor  Burst  communications  has  wide  ranging  effects  for  modern  society. 
It  offers  simple  and  reliable  communications  and  can  operate  in  the  simplex, 
half-duplex  and  full-duplex  mode  over  distances  of  200  to  1200  miles.  [9,22], 
Efficient  use  of  this  channel,  will  result  in  increased  data  throughput  for 
applications  such  as  facsimile  transmission  (FAX) ,  TELEX,  transmission  of  computer 
data,  radio  amateur  usage,  voice  burst  transmission  using  phonomes  (i.e., 
transmission  in  which  the  voice  is  digitally  encoded  using  rates  of  50  to  250 
bits  per  second)  and  military  communication  systems.  These  types  of  communications 
and  data  transmission  facilities  involve  transactions  of  banks,  universities  and 
medical  centers,  corporate  business  concerns  such  as  online  systems,  and  gov¬ 
ernment  agencies.  In  addition,  recreational  and  pleasure  or  leisure  activities 
could  make  use  of  this  mode  of  communication. 

If  shown  to  be  feasible  and  commercially  cost  effective  business  as  well 
as  home  users  could  provide  a  potentially  large  market.  The  advent  of  the 
microprocessor  will  reduce  the  cost  of  such  systems  thus  enabling  Inexpensive 
implementation  of  the  control  requirements  to  be  met  [9,22,23].  This  potential 
for  low  cost  will  open  up  the  home  and  business  market  which  has  already  invested 
millions  of  dollars  in  personal  computers  and  modems.  Furthermore,  the  number 
of  these  types  of  users  is  increasing  daily. 

Some  of  the  government  agencies  interested  in  alternate  communications 
means  are  the  Department  of  Agriculture,  the  Department  of  Energy  (DOE),  the 
National  Oceanic  and  Atmospheric  Administration  (NOAA) ,  the  Department  of  Defense 
(DOD),  the  National  Aeronautics  and  Space  Administration  (NASA)  and  the  Defense 
Communications  Agency-(DCA)  [9,22]. 

1.4  PhysUfll  Properties  of  Meteor  Trails 

Nearly  all  of  the  present  knowledge  regarding  the  physical  properties  of 
meteors  has  been  obtained  from  visual,  optical,  and  radio  observations  of  the 
trails  formed  by  the  meteoric  particles  as  they  enter  the  atmosphere  of  the 
earth.  In  this  section  the  properties  of  the  meteoric  particles  themselves  and 
the  trails  which  they  form  will  be  reviewed.  The  mathematical  model  and  conceptual 
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framework  reflecting  physical  and  radio  properties  of  the  trail  will  be  presented 
in  chapter  2.  The  review  of  meteoric  particles  and  the  trails  they  form  in 
sections  1.4. 1-1. 4. 4  is  taken  directly  from  Sugar's  [1]  1964  review  paper. 

1.4.1  Meteoric  rertlcjAa 

For  the  purposes  of  this  discussion  the  term  "Meteors"  will  be  used  to 
apply  to  those  particles  entering  the  earth's  atmosphere  that  are  burned  by 
frictional  heating.  This  definition  includes  the  very  small  particles,  the 
micrometeorites,  that  slowly  settle  through  the  atmosphere  without  being 
destroyed.  It  also  includes  the  large  meteors  which  manifest  themselves  as 
fireballs  or  the  even  larger  ones  which  reach  the  earth's  surface  as  meteorites. 
The  micrometeorites  are  not  of  concern  because  even  though  they  are  the  most 
numerous  of  the  various  types,  they  enter  the  atmosphere  too  slowly  to  produce 
any  significant  ionization.  The  large  meteoroids,  although  they  produce  sub¬ 
stantial  ionization,  are  of  little  concern  here  because  their  rate  of  occurrence 
is  extremely  low. 

Some  of  the  properties  of  meteoric  particles  are  summarized  in  Table  1 
[24],  [25].  The  particles  of  interest  in  meteor  propagation  are  those  with 
masses  In  the  range  103  to  10"7  g  and  dimensions  in  the  range  8  cm  to  40  microns. 
Before  being  trapped  by  the  gravitational  field  of  the  earth  these  particles 
move  in  orbits  around  the  sun.  Their  composition  is  uncertain;  however,  they 
appear  to  be  almost  entirely  of  cometary  origin.  A  substantial  fraction  of  them 
are  not  single  solid  particles  but  fragile  loosely-bound  agglomerates,  sometimes 
called  "dustballs", 

The  meteors  can  be  divided  into  two  classes,  the  shower  meteors  and  the 
sporadic  meteors.  The  shower  meteors  are  collections  of  particles  all  moving 
at  the  same  velocity  in  fairly  well-defined  orbits  or  streams  around  the  sun. 
Their  orbits  intersect  the  orbit  of  the  earth  at  a  specific  time  each  year  and 
at  these  times  the  well-known  meteor  showers  are  observed.  In  the  cases  where 
the  particles  are  uniformly  distributed  around  the  orbit,  the  size  of  the  shower 
varies  little  from  year  to  year.  If,  on  the  other  hand,  there  is  a  concentration 
of  particles  within  the  orbit,  the  extent  of  a  shower  can  vary  substantially  in 
successive  years.  (The  Leonid  shower  of  1833  is  an  example  of  an  unusual  display 
of  meteor  activity.  During  the  peak  of  the  shower  as  many  as  20  meteor  trails 
were  often  seen  by  a  single  observer  in  one  second  .)  The  shower  meteors,  while 
the  most  spectacular,  account  for  only  a  small  fraction  of  all  meteors.  It  is 
the  nonshower  or  sporadic  meteors  that  comprise  nearly  all  of  the  meteors  of 
interest  in  radio  propagation.  These  meteors  are  those  which  do  not  have 
well-defined  streams  but  rather  seem  to  move  in  random  orbits.  Thus,  whereas 
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shower  meteors  appear  to  be  coming  from  a  specific  point  in  the  sky- -the  radiant 
point  for  the  shower- -sporadic  meteors  have  radiants  that  appear  to  be  randomly 
distributed  over  the  sky. 

The  relation  between  sporadic  and  shower  meteors  is  not  clear  at  the  present 
time.  It  is  possible  that  the  sporadic  meteors  represent  the  final  stages  in 
the  decay  of  meteor  streams.  On  the  other  hand  there  is  some  indication  that 
the  sporadic  meteors  are  in  fact  distributed  in  a  large  number  of  relatively 
small  groups  and  that  the  earth  is  immersed  in  about  11  such  groups  at  a  time. 
Further  work  is  needed  to  resolve  the  origins  of  the  sporadic  meteors  [1], 

Shower  meteors  are  most  easily  recognized  in  terms  of  their  radiants, 
velocities,  and  time  of  occurrence,  since  these  parameters  are  relatively  fixed. 
The  radiants  and  times  of  occurrence  of  sporadic  meteors  are  random.  Their 
radiant  points  are  not,  however,  uniformly  distributed  in  the  sky  but  are  for 
the  most  part  concentrated  toward  the  ecliptic  plane  (the  plane  of  the  earth's 
orbit)  and  move  in  the  same  direction  around  the  sun  as  the  earth  moves  [26], 
The  orbits  are  not  uniformly  distributed  along  the  earth's  orbit  but  are  con¬ 
centrated  so  as  to  produce  a  maximum  incidence  of  meteors  at  the  earth  in  July 
and  a  minimum  in  February.  This  variation  in  the  space  density  of  meteors  is 
shown  in  Fig.  1  which  has  been  adapted  from  the  data  of  Hawkins  [26], 

The  rate  of  incidence  of  sporadic  meteors  at  the  earth  is  further  modified 
by  two  factors.  The  first  of  these- -resulting  in  a  regular  diurnal  variation 
in  meteor  rate  is  illustrated  in  Fig.  2  .  On  the  morning  side  of  the  earth, 
meteors  are  swept  up  by  the  forward  motion  of  the  earth  in  its  motion  around 
the  sun.  On  the  evening  side  the  only  meteors  reaching  the  earth  are  those  which 
overtake  it.  This  results  in  a  maximum  occurrence  rate  around  6  a.m.  and  a 
minimum  rate  around  6  p.m.  The  ratio  of  maximum  to  minimum  depends  on  the 
latitude  of  the  observer.  A  further  minor  seasonal  variation  is  introduced 
because  of  the  tilt  of  the  earth's  axis  relative  to  the  ecliptic  plane  [27]. 
This  variation,  also  dependent  on  the  latitude  of  the  observer,  can  change  the 
expected  hourly  rates  by  factors  as  large  as  1.4:1. 

The  mass  distribution  of  sporadic  meteors  is  such  that  there  are  approx¬ 
imately  equal  total  masses  of  each  size  of  particle.  There  are,  for  example, 
10  times  as  many  particles  of  mass  10'*  g  as  there  are  particles  of  mass  10~3 
g.  This  approximate  relation  between  particle  mass  and  number  is  given  in  Table 
1.  The  mass  distribution  of  shower  meteors  is  somewhat  similar  to  that  for  the 
sporadic  meteors  with  the  important  difference  that  there  are  more  large  particles 
relative  to  the  number  of  smaller  ones  than  for  the  sporadic  meteors  . 
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The  velocities  of  meteors  approaching  the  earth  are  in  the  range  11.3  to 
72  km/sec.  The  lower  limit  is  the  escape  velocity  for  a  particle  leaving  the 
earth  and  is  therefore  the  lowest  velocity  that  a  particle  falling  toward  the 
earth  can  have.  The  upper  limit  is  the  sum  of  two  components,  a  30  km/sec 
component  associate  j  with  the  earth  in  its  orbit  around  the  sun,  and  a  42  km/sec 
component  associated  with  the  meteor  itself.  This  latter  velocity  is  the 
escape  velocity  for  a  particle  leaving  the  solar  system.  Nearly  all  observations 
have  indicated  that  meteor  velocities  fall  in  the  above  range  and  that  the  meteors 
are  thus  members  of  the  solar  system. 

1.4.2  Meteor  Trails 

1)  Formation:  As  a  meteoric  particle  enters  the  earth's  atmosphere  it 
collides  with  air  molecules  which  then  become  trapped  in  its  surface.  The  impact 
energy  produces  heat  which  evaporates  atoms  from  the  meteor  and  these  move  off 
with  a  velocity  which  is  substantially  equal  to  that  of  the  meteor.  Collisions 
between  these  high  velocity  atoms  and  the  surrounding  air  result  in  the  production 
of  heat,  light,  and  ionization,  distributed  in  the  form  of  a  long,  thin,  paraboloid 
of  revolution  with  the  particle  at  the  head  of  it.  The  electron  line  density 
in  the  trail  is  proportional  to  the  mass  of  the  particle.  In  the  evaporation 
process  each  original  impact  frees  many  meteoric  atoms  and  thus  the  total  mass 
of  air  molecules  striking  the  meteor  is  small  compared  to  the  meteoric  mass 
itself.  As  a  consequence  the  velocity  of  the  meteor  remains  quite  constant  until 
the  m-steor  is  nearly  completely  evaporated.  (Those  meteors  that  appear  to  be 
dustballs  rather  than  solid  particles  break  up  or  "fragment"  upon  entering  the 
atmosphere  and  proceed  as  a  group  of  smaller  particles.  For  this  case  the  above 
description  appears  to  be  correct  if  it  is  applied  to  the  individual  particles 
after  fragmentation  ocurs.) 

2)  Heights :  As  a  meteoric  particle  approaches  the  earth  no  appreciable 
ionization  if  formed  until  the  particle  enters  the  relatively  dense  air  at  heights 
below  about  120  km.  Above  that  height  collisions  of  the  particles  with  air 
molecules  are  not  frequent  enough  to  be  of  significance.  As  the  particle  traverses 
the  region  below  about  120  km  it  vaporizes  rapidly  and  most  particles  are 
completely  evaporated  before  reaching  80  km.  The  relatively  small  thickness  of 
the  meteor  region  is  a  result  of  the  rapid  change  in  air  density  which  occurs 
over  the  height  range  quoted.  At  120  km  the  mean  free  path  is  5.4  m  and  this 
decreaes  in  an  approximately  logarithmic  manner  so  that  at  80  km  it  is  only  3.8 


mm. 
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The  height  distribution  of  trails  varies  somewhat  with  particle  charac¬ 
teristics.  The  higher  velocity  particles  produce  trails  at  the  higher  heights 
with  the  mean  trail  height  increasing  to  about  10  km  as  velocity  increases  from 
less  than  15  km/sec  to  greater  than  60  km/sec  [1]  .  The  particles  of  higher  mass 
produce  maximum  trail  ionization  at  lower  heights.  Over  the  mass  range  of  103 
to  10~7  g  the  height  variation  is  around  44  km.  There  is  also  a  height  variation 
with  the  zenith  angle  of  the  trail  orientation,  with  the  larger  angles  corre¬ 
sponding  to  greater  heights.  A  variation  of  about  13  km  is  associated  with 
nearly  the  whole  range  of  zenith  angles. 

3)  Lengths:  The  lengths  of  meteor  trails  are  primarily  dependent  on 
particle  mass  and  zenith  angle.  Typical  lengths  range  up  to  50  km  with  the  most 
probable  trail  length  for  sporadic  meteors  being  15  km.  (Several  definitions 
of  "length"  can  be  used  and  the  one  chosen  here  uses  as  end  points  the  points 
with  a  threshold  value  chosen.) 

4)  Initial  Radius:  Until  recently  it  has  been  assumed  that,  at  the  time 
of  their  formation,  trails  had  an  initial  radius  of  the  order  of  the  mean  free 
path  or  at  most  about  14  times  this  radius.  However,  photographic  and  radio 
measurement"  have  suggested  that  the  initial  radius  is  significantly  larger  than 
indicated  above  and  that  this  increase  is  probably  associated  with  the  dustball 
and  fragmentation  hypothesis. 

The  trail  radii  indicated  by  these  measurements  are  in  the  range  0  to  1.2 
m  (with  a  mean  value  of  0.65  m)  for  the  photographic  work,  and  0.55  to  4.35  m 
for  the  radio  work.  The  0.55  m  was  for  a  height  of  81  km  where  the  mean  free 
path  is  only  5  x  10  m.  The  results  at  121  km  are  more  nearly  consistent  with 
Manning's  hypothesis  since  the  initial  radius  observed  is  4.35  m  and  the  mean 
free  path  is  5  m. 

5)  Dissipation:  Once  the  trail  is  formed  it  expands  by  diffusion  [21] 
at  a  relatively  low  rate,  producing  a  radial  distribution  of  material  that  is 
approximately  Gaussian.  The  quantity  (4  Dt  +  r  )  may  be  taken  as  the  approximate 
radius  of  the  trail  after  a  time  t  where  D  is  the  diffusion  coefficient  and  r 
is  the  initial  radius  of  the  trail.  D  varies  from  1  m2/sec  at  an  85  km  height 
to  140  m2/sec  at  115  km.  Thus  after  one  second  a  trail  will  have  a  radius  in 
the  range  2  to  20  m. 

The  practical  lifetime  of  a  trail  is  of  course  dependent  on  the  means  for 
detecting  it.  Most  trails  detected  by  radio  means  are  those  resulting  from  small 
particles  and  these  last  only  a  fraction  of  a  second.  The  larger  particles 
produce  more  densely  ionized  trails  and  trails  with  durations  of  the  order 
of  a  minute  are  observed  several  times  per  day.  Trails  with  durations  of 
the  order  of  an  hour  or  more  are  extremely  rare. 


The  dissipation  of  trails  is  further  complicated  by  the  presence  of  winds 
in  the  meteor  region.  At  the  time  of  formation  trails  are  auite  straight  but 
they  are  rapidly  deformed  by  these  winds  which  have  typical  velocities  of  the 
order  of  25  m/sec  and  vertical  gradients  with  mean  relative  maxima  of  about  100 
m/sec/km.  A  wind  shear  of  this  magnitude  can  rotate  part  of  a  trail  through  an 
angle  of  5°  in  one  second. 

1.4.3  Reflections  Properties  of  Trails 

The  distribution  of  energy  reflected  by  a  meteor  trail  is  a  function  of 
many  variables.  The  ionization  density  distribution  across  and  along  the  trail, 
the  orientation  of  the  trail,  the  radio  wavelength,  the  polarization  of  the 
incident  wave  relative  to  the  trail,  motion  of  the  trail  either  as  part  of  the 
process  of  formation  or  due  to  ionospheric  winds,  and  the  straightness  of  the 
trail  are  all  significant.  In  discussing  the  reflection  properties  it  is 
convenient  to  divide  the  trails  into  two  classes  underdense  trails  and  overdense 
trails  and  to  examine  the  properties  of  each  class  independently.  Underdense 
trails  are  those  wherein  the  electron  density  is  low  enough  so  that  the  incident 
wave  passes  through  the  trail  and  the  trail  can  be  considered  as  an  array  of 
independent  scutterers.  Overdense  trails  are  those  wherein  the  electron  density 
is  high  enough  to  prevent  complete  penetration  of  the  incident  wave  and  to  cause 
reflection  of  waves  in  the  same  sense  that  the  ordinary  ionospheric  reflections 
occurs.  A  rough  sorting  of  trails  into  these  two  categories  can  be  done  on  the 
basis  of  trail  lifetime  or  duration.  At  long  wavelengths  the  underdense  trails 
have  durations  of  less  than  about  one  second  while  the  overdense  trails  have 
longer  durations.  At  long  wavelengths  the  effective  duration  of  a  trail  is  large 
compared  to  the  time  it  takes  the  trail  to  form,  and  the  trail  may  be  considered 
to  have  a  cylindrical  shape. 

The  solutions  to  be  considered  here  are  useful  approximations  to  the 
physical  problems.  They  will  apply  quite  well  to  some  of  the  trails  observed 
and  rather  poorly  to  others.  A  complete  analysis  of  the  reflections  from  even 
a  relatively  simple  trail  would  be  far  too  complex  to  be  of  any  practical  use 
or  interest. 


1.4.4  Other  Aspects 

In  this  section  some  additional  practical  aspects  will  be  discussed  from 
a  qualitative  viewpoint.  The  discussion  is  directed  toward  the  long  wavelength 
cases  since  relatively  few  short  wavelength  observations  have  been  made  and 
analyzed. 
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1)  Long  Wavelength  Reflections  during  Trail  Formation :  The  transient 
state  associated  with  trail  formation  is  of  interest  since  it  accounts  for  some 
of  the  observed  characteristics  of  trail  reflections.  As  a  trail  is  being  formed, 
but  before  the  meteor  reaches  the  first  Fresnel  zone,  a  weak  reflection  is 
obtained  from  the  incomplete  trail.  This  comes  primarily  from  the  part  of  the 
trail  corresponding  to  the  shortest  transmission  path  at  the  instant;  and  when, 
as  is  usually  the  case,  this  is  the  head  of  the  trail,  the  reflected  signals 
are  shifted  in  frequency  because  of  the  motion  of  the  effective  reflecting  point. 
As  the  meteor  approaches  the  first  Fresnel  zone  for  the  trail  this  frequency 
shift  approaches  zero,  and  thus  the  received  frequency  decreases  with  time.  The 
observed  frequercy  will  of  course  depend  on  trail  orientation,  meteor  velocity 
and  observing  wavelength.  A  maximum  shift  of  the  order  of  5  kHz  is  possible 
for  50  MHz  forward- scatter  observations  over  a  1000  km  path. 

2)  Trail  Drift  and  Distortion:  The  effects  of  ionospheric  winds  are 
appreciable  for  trails  which  last  for  the  order  of  a  second  or  more.  A  Doppler 
shift  of  the  received  frequency  will  be  associated  with  the  average  wind  velocity 
of  the  trail.  For  a  velocity  of  25  m/sec  this  "body  Doppler"  can  be  as  large 
as  18  Hz  for  backscatter  observations  at  50  MHz  and  will  be  somewhat  less  for 
forward- scatter  observations . 

The  trail  distortion  resulting  from  wind  shears  can  lead  to  the  formation 
of  several  local  first  Fresnel  zones  for  the  trail  since  a  distorted  trail  can 
have  several  points  where  the  transmission  path  length  has  a  local  minimum. 
These  local  minima,  or  "glints"  as  they  have  been  called  [41] ,  are  strong  scatters 
and  the  received  signal  is  a  composite  of  their  contributions.  Since  they  are 
moving  at  different  velocities  the  signals  from  each  have  different  Doppler 
shifts  and  the  resultant  composite  signal  fades  in  an  irregular  manner  [24], 
(At  50  MHz  the  fading  rate  observed  for  forward  scatter  are  of  the  order  of  1-10 
Hz.)  In  addition  to  producing  the  fading  observed  for  long-enduring  trails  the 
wind  shears  can  rotate  a  trail  sufficiently  to  produce  reflections  when  the 
initial  orientation  has  not  been  suitable.  Thus  these  trails  lose  their  aspect 
sensitivity  as  time  goes  on,  and  if  their  life  is  of  the  order  of  10  seconds 
they  will  scatter  in  all  directions . 

3)  Polarization  Effects:  Thus  far  the  only  indication  that  the  polarization 
of  the  incident  wave  can  be  of  significance  was  the  inclusion  of  the  sin2  (a) 
term  in  each  of  the  transmission  equations  for  the  forward- scatter  cases.  This 
term  accounts  for  the  foreshortening  of  the  electric  vector  of  the  incident  wave 
which  occurs  when  it  is  viewed  from  the  receiver.  In  addition  to  this,  a  kind 
of  electron  reasonance  can  occur  which  increases  the  reflection  coefficient  for 
the  trail.  This  is  associated  with  the  restoring  force  that  the  electrons 
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experience  as  the  incident  electric  field  displaces  them  from  their  equilibrium 
positions  in  a  direction  normal  to  the  trail.  This  reasonance  can  occur  only 
for  underdense  trails  and  only  then  for  trail  diameters  much  less  than  the 
wavelength.  This  reasonance  at  most  doubles  the  reflection  coefficient  and  thus 
can  increase  the  received  power  by  a  factor  of  4. 

4)  Diversity  Effects:  Underdense  trails  act  as  small  coherent  sources, 
and  therefore  CV  signals  scattered  by  them  exhibit  good  space  and  frequency 
correlation.  Correspondingly,  very  little  pulse  broadcasting  is  observed  when 
pulse  signals  are  used.  The  limited  available  data,  suggest  that  for  a 
forward-scatter  path  about  1000  Van  in  length,  operating  near  50  MHz  the  correlation 
coefficent  observed  for  single  underdense  trails  will  fall  to  0.5  for  antenna 
spacings  of  the  order  of  150  km  along  the  path,  for  antenna  spacings  of  the  order 
of  30  km  across  the  path,  and  for  frequency  separations  greater  than  5  MHz. 
Three  micro-sec  pulses  show  no  appreciable  broadening  under  these  conditions. 

Overdense  trails,  in  contrast  to  the  underdense  trails,  tend  to  act  like 
relatively  large  sources  and  therefore  exhibit  much  poorer  space  and  frequency 
correlation  properties.  Again,  as  in  the  underdense  case,  specific  data  on  these 
properties  are  not  available.  Present  results  suggest  that  the  correlation 
observed  for  single  overdense  trails  will  fall  to  0.5  for  antenna  spacings  of 
the  order  of  50  [41].  Measurements  using  short  pulses  over  a  1000-km  path 
indicate  that  the  received  composite  signal  can  have  a  total  time  spread  as  large 
as  10  micro-  sec  [44]. 

i-s  sm  Samldfeifttlsm 

There  are  three  potential  modes  of  operation  for  a  meteor  burst  system: 
point-to-point,  netted,  and  broadcast.  While  a  meteor  burst  system  could  be 
operated  in  any  or  all  of  these  modes,  all  known  meteor  burst  systems  that  have 
currently  been  implemented  have  been  designed  for  point-to-point  applications. 
In  this  section,  we  discuss  some  of  the  possibilities  for  configuring  meteor 
burst  systems  in  point  to  point  mode.  The  general  discussion  in  the  following 
sections  1.5. 1-1. 5. 3  is  taken  directly  from  Oetting's  [9]  1980  review  paper. 

1.5.1  Point  to  Point  Mode 

The  point-to-point  mode  is,  of  course,  straightforwardly  implemented.  The 
only  requisite  for  effective  system  control  is  the  ability  of  the  transmitting 
terminal  to  discern,  as  accurately  as  possible,  the  beginning  and  end  of  a  useful 
burst.  If  transmission  begins  too  late  or  ends  too  soon,  valuable  burst  time 
will  be  wasted.  If  the  transmission  extends  past  the  useful  portion  of  the 
burst,  a  high  error  rate  will  result. 
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In  all  point* co-point  systems,  we  assume  that  a  feedback  path  is  available. 
In  the  half -duplex  case,  the  forward  and  return  links  share  the  same  frequency, 
while  in  the  full-duplex  case,  they  use  two  different  frequencies.  For  achieving 
data  transmission  between  point  A  and  point  B,  perhaps  the  simplest  strategy  is 
to  assign  a  different  frequency  to  each  transmitter  and  to  continuously  transmit 
a  probe  signal  from  the  master  station  (say,  point  A).  When  point  B  detects 
the  probe,  he  transmits  a  preamble  followed  by  his  information.  Point  A  uses 
the  preamble  to  synchronize  his  receiver  and  is  subsequently  able  to  receive 
data.  It  is  usually  required  that  A  acknowledge  the  reception  of  the  message 
or  block  of  data. 

Several  variations  of  this  basic  point-to-point  approach  have  been 
developed.  Early  meteor  burst  systems  used  a  signal  amplitude  threshold  to 
determine  if  the  probe  signal  were  being  received.  The  JANET  system  [1] - [13] 
improved  upon  this  procedure  by  substituting  a  signal-to-noise  ratio  measurement 
for  the  fixed  threshold  test,  thereby  reducing  the  false  alarm  rate.  After 
extensive  testing  of  the  JANET  B  system,  the  SHAPE  Technical  Centre  concluded 
that,  even  with  the  SNR  approach,  a  suitable  compromise  between  efficient  use 
of  the  channel  and  low  error  rate  was  not  achievable.  Thus,  they  developed  their 
own  approach  using  ARQ  with  one  teletype  character  per  seven-bit  block.  This 
apporach  is  referred  to  as  the  COMET  system  (5]  and  [9], 

The  possibilities  for  netted  meteor  burst  systems  have  not  been  well 
explored  to  date.  The  only  known  system  capable  of  netted  operation  is  the 
Western  Union  hardware  currently  being  installed  for  the  Department  of  Agriculture 
[28].  This  system  will  consist  of  511  remote  sites  that  communicate  with  two 
master  stations.  Although  the  Department  of  Agriculture  application  involves 
one-way  transmission  of  data  from  the  remote  to  the  master  stations,  the  system 
can  be  configured  to  provide  two-way  data  transmission.  A  network  could  then 
be  achieved  by  having  a  master  station  operate  in  a  store -and- forward  mode. 

1.5.2  Low  toMlUfl  2f  Intercept  (LPI1  ftng 
AnUJffm  iMl  Consideration 

One  of  the  primary  attractions  of  meteor  burst  communications  is  its 
inherent  privacy  resulting  from  the  restricted  footprir.t  of  reflections  from 
meteor  trails.  This  property  has  received  considerable  attention  in  the  meteor 
burst  literature  and  has  prompted  a  large  number  of  industry  proposals  to  provide 
LPI  meteor  burst  systems  for  numerous  Department  of  Defense  (DoD)  and  non-DoD 
agencies . 
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Signals  propagating  between  a  specific  transmitter  and  receiver  can  normally 
be  detected  by  an  eavesdropping  receiver  only  if  the  latter  is  located  in  an 
elliptically  bounded  "footprint"  in  the  vicinity  of  the  primary  receiver.  By 
reciprocity,  interfering  signals  from  a  remote  transmitter  will  propagate  to  a 
meteor  burst  receiver  only  if  the  interfering  transmitter  is  within  the  footprint 
of  the  desired  transmitter.  Measurements  indicate  that  the  major  and  minor  axes 
of  the  footprint  are  on  the  order  of  2000  and  25  mi,  respectively.  Clearly, 
the  footprint  is  not  a  well-defined  boundary  within  which  the  signal  is  always 
received  and  beyond  which  it  is  never  received.  Instead,  the  situation  is  better 
depicted  by  Fig.  3  which  shows  theoretical  and  experimental  results  for  the 
correlation  between  the  signals  detected  by  two  different  receivers  as  a  function 
of  the  distance  between  them.  When  the  unintended  receiver  is  moved  in  a  direction 
perpendicular  to  the  transmission  path,  the  correlation  falls  off  more  rapidly 
than  when  the  unintended  receiver  is  moved  along  the  propagation  path.  [9] 

While  the  theoretical  along- the -path  correlation  falls  off  at  an  alarmingly 
slow  rate,  the  experimental  results  illustrate  two  key  points:  the  theoretical 
performance  is  difficult  to  achieve  in  practice,  and  the  correlation  decreases 
at  higher  operating  frequencies.  In  addition,  it  has  been  shown  [29]  that  fading 
reflections  are  responsible  for  a  significant  part  of  the  correlation  at  large 
separations  (trails  that  have  been  distorted  by  winds  exhibit  broader  reradiation 
patterns  than  newly  formed  trails).  By  communicating  only  during  the  nonfading 
portion  of  each  burst,  the  size  of  a  footprint  can  be  kept  to  a  minimum. 

LPI  is  best  achieved  by  operating  in  a  burst  mode  with  a  low  duty  cycle. 
While  interception  by  a  beyond -the -horizon  transmitter  is  then  governed  by  the 
probabilities  shown  in  Fig.  3,  a  meteor  burst  system  is  vulnerable  to  line-of -sight 
interception  unless  spread  spectrum  techniques  are  employed. 

As  for  AJ  considerations,  line-of -sight  jamming  can  be  easily  accomplished 
in  view  of  the  losses  incurred  by  the  desired  signal  in  the  course  of  being 
reflected  from  a  meteor  trail.  In  the  case  of  beyond- the -horizon  jamming,  the 
jammer  can  distrupt  only  a  small  fraction  of  the  traffic  if  he  relies  on  reflections 
from  meteor  trails  to  propagate  his  jamming  signals.  A  much  better  stategy  is 
to  use  extremely  high  power  (perhaps  in  a  pulsed  model)  and  rely  on  ionoscatter 
propagation.  This  approach  can  be  defeated  by  operating  at  the  higher  meteor 
burst  frequencies,  since  ionoscatter  falls  off  much  faster  with  increasing 
frequency  than  does  meteor  trails  reflections.  Finally,  a  jammmer  can  attempt 
to  spoof  the  system  by  transmitting  a  simulated  probe  signal,  but  spoofing  can 
be  defeated  by  judicious  system  design. 
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1.5.3  System  Performance 

The  performance  parameters  of  greatest  Interest  in  commercial  applications 
are  the  average  throughput  (in  bit  per  second),  the  average  waiting  time  (to 
transmit  a  messge  of  specific  length),  and  the  probability  of  error.  These 
parameters  will  be  influenced  to  varying  degrees  by  the  environmental  variables : 
background  noise  level,  external  interference,  geographic  location,  season,  and 
time  of  day.  The  system  designer,  on  the  other  hand,  can  exert  control  over 
system  performance  by  manipulating  his  design  variables  -  transmitted  power,  RF 
frequency,  instantaneous  data  rate,  modulation  technique  and  antenna  gain. 

The  throughput  of  a  meteor  burst  system  depends  on  factors  such  as  transmitted 
power,  operating  frequency,  communications  range,  received  noise  and  instan¬ 
taneous  data  rate.  The  increase  in  throughput  that  can  ultimately  be  attained 
is  limited  by  intersymbol  interference  due  to  multipath  effects,  practical  bit 
synchronization  problems,  and,  of  course,  the  available  transmitter  power  and 
antenna  gains. 

In  addition  to  the  factors  under  the  control  of  the  system  designer,  the 
throughput  also  depends  on  such  environmental  variables  as  time  of  day  and  season 
of  the  year.  These  variations  typically  involve  a  factor  of  five  diumally  and 
three  seasonally. 

1.6  Past  Systems 

A  landmark  meteor  trail  communication  system  was  operational  in  the 
1950s- -the  Canadian  JANET  system.  It  was  used  for  teletype  communications  between 
Toronto  and  Port  Arthur,  a  distance  of  approximately  1000  km  (1] , [4] .  A  double 
sideband  AM  duplex  channel  was  used  to  transmit  and  receive  cosine- squared  shaped 
pulses  in  pulse  position  modulation  (ppm)  at  a  rate  of  650  bits  per  second  (bps) . 
Punched  paper  tape,  magnetic  tape,  and  a  toriodal  magnetic  matrix  storage  core 
were  used  as  storage  media  to  compress  messages  in  time  for  transmission  and  to 
store  and  expand  the  messages  at  the  receiving  station.  Despite  its  relative 
antiquity,  the  JANET  system  embodied  many  of  the  features  of  today's  systems: 
stored  digital  data  which  is  transmitted/received  in  bursts  after  carrier 
detection,  duplex  operation  at  VHF  frequencies  around  50  MHz,  and  low  effective 
duty  cycles  (about  0.1  for  the  JANET  system).  Other  links,  such  as  a  one  way 
Bozeman,  MONTANA -STanford,  California  meteor  burst  link  were  also  operational 
during  the  1950s,  but  the  JANET  system  represents  the  first  mature,  complete, 
practical  hardware  implementation  of  a  meteor  burst  system. 
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Subsequent  experimental  tests  were  conducted  during  the  1960s  and  1970s. 
Notable  among  these  efforts  Is  the  experimental  work  summarized  by  Sugar  [4]  and 
the  Implementation  of  the  COMET  (Communication  by  MEteor  Trails)  system  which 
has  operated  between  the  Netherlands  and  Southern  Prance.  COMET  made  use  of 
frequency  diversity  and  Automatic  Repeat  for  reQuest  (ARQ)  which  allowed  for 
signal  repetition  in  the  event  of  no  radio  path.  The  COMET  system  used  a  2000 
baud  signaling  rate  and  FSK  with  a  6  KHz  deviation.  The  storage  devices  were 
similar  to  those  of  the  JANET  system.  COMET  demonstrated  the  practicality  of 
meteor  burst  communications  under  a  variety  of  conditions.  Typical  results 
indicated  successful  transmission  of  a  50  baud  (average)  message,  150  characters 
long  with  a  delay  of  about  a  minute.  Worst  case  delays  of  3  or  4  minutes  were 
typical  [5).  Studies  of  COMET-type  systems  have  continued  into  the  1980s  [5] [9]. 


15 


2  MATHEMATICAL  MODEL  AND  CONCEPTUAL  FRAMEWORK 

In  this  chapter  the  basic  equations  that  model  the  meteor  burst  channel 
are  presented.  In  addition,  basic  assumptions,  constraints  and  empirical  factors 
are  stated. 

2.1  Pover  Received  and  Related  Parameters 

The  power  received  from  the  trail  model  used  here  is  from  the  classic  paper 
by  Eshelman  and  Manning  [6]  and  Hines  and  Forsyth  [16]  . 

2.1.1  Power  Received  from  an  Pnderdense  Trail 

Meteor  trails  with  fewer  then  101*  electrons  per  meter  of  length  are  referred 
to  as  underdense.  Eshelman  has  shown  that  a  simple  Fresnel  description  of  the 
process  is  normally  sufficient.  The  scattering  then  is  specular  for  the  majority 
of  recevied  signals  that  is  the  incident  and  scattered  rays  make  equal  angles 
6 with  the  axis  of  the  trail  and  the  received  signal  comes,  in  effect,  from  the 
principal  Fresnel  zone  alone.  See  Fig.  4 

The  principal  Fresnel  zone  is  centered  on  a  point  P  in  the  trail  at  a 
distance  RT  and  R*  from  the  transmitter  T  and  receiver  R  respectively  such  that 
RT  +  Rr  is  a  minimum  for  the  trail  in  question.  The  received  power  rises  rapidly 
in  the  fraction  of  a  second  taken  by  the  meteor  to  traverse  the  principal  zone 
and  reaches  a  peak  volume  at  time  t-0. 

Although  the  meteor  trail  is  formed  as  a  narrow  column,  a  few  centimeters 
in  diameter,  it  immediately  begins  to  expand  by  diffusion.  As  the  diameter 
increases,  the  scattered  signal  suffers  increasingly  from  a  destructive 
interference  of  the  fields  scattered  by  individual  electrons.  In  ideal  cir¬ 
cumstances,  neglecting  Fresnel  ripples  and  initial  diameter,  the  received  power 
decays  exponentially  in  time,  see  fig.  7.  The  above  is  summarized  in  the  next 
equation. 


no- 


PrCTG,risin2(,a)\3exp(-^-~) 
l6n2RTR,(RT*  *,)(!  - cos^sin2*)9  6XP 


32n 2P  \ 

\zsecz$  J 


'where. ' 


PT  Transmitter  power. 

Gt,Gr  Transmitter  and  receiver  antenna  gains,  respectively, 
r.  Classical  radius  of  an  electron  (2.82  10*ls  m) . 


(2. 1.1-1) 
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a  Angle  between  the  electron  field  vector  E  and  R R. 

q  Electron  line  density  of  the  trail. 

X  Carrier  wavelength. 

r0  Initial  radius  of  the  trail  (.65  m.). 

4>  Angle  of  incidence/reflection  of  the  transmitted  plane  wave. 

Rt  Distance  from  transmitter  to  trail. 

R%  Distance  from  receiver  to  trail. 

D  Diffusion  coefficient  of  the  atmosphere. 

|3  Angle  between  the  trail  and  the  propagation  plane  formed  by  Rr  and  Rr. 


Figures  4,  5,  and  6  depict  the  geometrical  parameters. 

A  minimal  set  of  assumptions  is  applied  in  the  literature  [7,11,21]  for 
analysis : 

1. )  The  trail  occurs  at  midpoint  between  receiver  and  transmitter  i.e.  Rt-Rr 

2. )  The  burst  occurs  at  an  altitude  of  100  km. 

3. )  The  trail  is  travelling  at  a  plane  perpendicular  to  the  plane  formed  by  RT 

and  Rr  i.e.  p-n/2.  In  addition,  we  set  a-n/2 

Under  these  assumptions  and  after  substituting  for  the  physical  constants 
we  rearrange  the  above  equation  to  arrive  at  the  power  received  equation  for 
the  underdense  case . 


where , 


C u  -  2.5 1  79581  ♦  10 


-  _  _  %  q  (  33  339  \ 

-  32  PrCTC,\ ’exp^-T— ) 


(2.1  .1  -  2) 


B- 3.166-  I0'3 


X.2sec2$ 

D 


It  should  be  emphasized  that: 

1. )  P(0) ,  the  maximum  power  from  the  underdense  trail  is  proportional  to  q2  ; 
for  a  specific  trail  q  is  fixed.  (The  random  character  of  q  is  discussed  later) . 

2. )  C0  the  constant  of  proportionality,  incorporates  the  effects  of  link  geometry, 
transmitter  power,  antenna  gains  and  carrier  wavelength.  For  a  given  communication 
system  Cj,  is  constant  since  link  distance,  transmitter  power,  antenna  gains  and 
carrier  frequency  are  set. 
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3. )  The  decay  time  constant  B  given  above  is  found  to  be  randomly  varying  from 
trail  to  trail.  The  expression  in  the  classic  model  corresponding  to  a  decay 
time  constant  is  taken  as  the  average  decay  time  constant  in  accordance  with 
observations.  To  date,  however,  researchers  have  assumed  all  underdense  trails 
to  share  the  same  decay  time  constant  to  fascilitate  their  analysis.  It  is  in 
this  added  complexity  to  the  model  that  our  treatment  deviates  from  published 
analysis  and  enhances  its  accuracy. 

4. )  The  decay  time  constant  B  and  the  electron  line  density  q  are  functionally 
independent  [1]  [7]  and  hence  statistically  independent. 

5. )  For  a  given  communication  system  (i.e.  Cu  fixed)  knowledge  of  the  electron 
line  density  ,q,  and  decay  time  constant  ,B,  completely  specify  the  underdense 
trail  behavior  in  time. 


2.1.2  Electron  Line  Density  Statistics 
i  Pnderdense  Case 

Based  on  recent  empirical  data  [33]  the  pdf  for  q  is: 


where , 


f„{q)-Qqp  qml„^q$q*- io"e/m 


Q’(p-nqL' 


1 


(2.1  .2-2) 


p  m  1.6 

qm)11- Minimum  q  as  explained  below. 

q  Bln  is  the  minimum  electron  line  density  q  a  trail  must  posses  to  be 
’seen’  by  the  communication  system  i.e.  q2nln  *-s  the  minimum  power  level  Pmln 
detected  by  the  communication  system.  In  the  context  of  constant  transmission 
rate  systems  (  as  disscussed  later  )  Pmln  is  the  power  level  that  corresponds  to 
the  maximum  allowed  bit  error  rate.  See  fig.  8 
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2.1.3  D«crg  ZlBft  Constant 
-  Pndcrdensc 


Case 


We  start  by  stating  that  the  expected  value  of  the  decay  time  constant  B 
is  well  approximated  by: 

2sec2$ 


f(S)-  5-3.16610' 


D 


sec 


(2.1  .3-  1  ) 


as  defined  in  section  2.1.1 

The  random  decay  time  constant  B  is  assumed  to  be  exponential  or  Rayleigh. 
Both  assumptions  have  been  shown  to  be  consistent  with  recent  experimental 
observation  [33],  In  general  the  Rayleigh  distribution  will  be  used  since  for 
very  small  B  (  extremely  fast  decay  )  the  probability  is  very  near  zero.  This 
is  appealing  from  a  physical  viewpoint  and  consistent  with  observation.  It 
should  also  be  noted  that  even  if  such  exceedingly  fast  decaying  trails  occur 
they  are  of  no  utility  for  communication. 

Exponential : 

(2. 1.3-2) 

B 

where  Sis  defined  above. 

Rayleigh: 

(2. 1.3-3) 

a 

where , 


2.1.4  Joint  PDF  fif  a  and  fi  -Underdense  Case 

Since  the  electron  line  density  and  the  decay  time  constant  are  statistically 
independent  the  Joint  pdf  is  given  by  their  product  i.e. 

/,C«>)  (2.1.4-  1) 

where  the  individual  densities  are  given  in  the  previous  sections. 

2-1-5  Power  received  from  Overdense  Trails 

Trails  with  more  101*  electrons  per  meter  are  often  termed  "overdense." 
When  treating  them,  it  is  no  longer  satisfactory  to  think  of  the  incident  wave 
as  penetrating  the  trail  without  serious  modification.  Instead,  coupling  between 
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Individual  electrons  plays  a  prominent  part,  and  the  whole  of  the  scattering 
process  is  conceived  more  simply  in  terms  of  reflections  from  a  cylindrical 
surface . 

The  surface  in  question  surrounds  the  axis  of  the  trail  at  a  radius  r, 
where  the  electron  volume  density  has  the  critical  value  normally  associated 
with  reflection  of  the  incident  wave  (cf.  reflection  from  the  ionosphere).  The 
volume  density  within  the  trail  exceeds  the  critical  value,  and  cannot  support 
propagation  in  the  usual  sense  at  the  pertinent  frequency.  As  the  electrons 
diffuse  outwards,  r  increases  from  an  initially  small  value,  passes  through 
a  maximum,  and  falls  to  zero  as  the  volume  density  at  the  axis  of  the  trail 
decreases  below  the  critical  value.  In  the  very  latest  stages  the  scattering 
reverts  to  the  underdense  type,  but  only  after  some  time,  and  so  only  when  the 
underdense  scattering  is  extremely  small. 

The  scattering  process  is  again  specular,  the  angles  of  incidence  and 
reflection  being  equal.  The  received  power  variation  is  given  by: 


where , 


f(O«C'0 


I r\+  4Dt  (  r,\* sec2$ 
V  sec2$  V  rtz(ro  ♦  4Z>() 


C'o 


PTOjC,  sin2(a)X.2 

32n2 R T R „{R 7  *  R -  cos^sin2*) 


(2.1.5-  1) 


and  all  other  parameters  are  as  defined  for  the  underdense  burst. 

Using  the  same  assumptions  as  for  underdense  case  (  a-p-n/2  and  RT-Rt) 
and  some  algebra  we  rewrite  the  above  as: 

P(.0m  CoyJ 1  (*  *  O  •  j  l<aq-k  (2. 1.5-2) 

where , 

_ PrCTC,\z 

C0- 3.166286 \0'3fD ■-* — 

R3r 

k  -  .105625/D 


a 


7. 1  43 1 43  1  O' ,f 


\2sec2$ 

D 


As  in  the  underdense  case  we  have  rearranged  the  classic  equation  for  ease  of 
handling.  Note  that  C0  contains  the  transmitted  power,  antenna  gains,  diffusion 
coefficient,  carrier  frequqncy  and  the  effects  of  link  distance.  Fig.  7 
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In  addition,  from  eq.  (2. 1.1-2)  we  relate  <:„  and  C0  by 


rr  (  33.33»  \ 

1. 2574816  1029^  • 


(2.1 .5-3) 


2.1.6  Electron  Line  Density  Statistics 
^  Overdense  Case 

Based  on  recent  empirical  data  [33]  the  pdf  for  q  is: 


f  Q(.q)-Qq~P  10l6e/m 


(2.1  .6-2) 


where , 


Q- (P 


-(£)•" 


p-  i  .s 

q„,„  -  Minimum  q  as  explained  below. 

qaino  is  the  minimum  electron  line  density  q  an  overdense  trail  must  posses 
to  be  'seen'  by  the  communication  system.  Note  that  q^,,,,  >  qLo  the  (physically) 
lowest  possible  electron  line  density  for  an  overdense  trail,  q^-lO14  e/m.  We 
define  another  constant  which  we  refer  to  later  as: 

Q'-i/[i-(q»i„o/qLo)p'1]. 


2.2  Link  Geometry 

In  this  section  we  present  the  relationship  between  link  distance,  L,  and 
station  to  trail  distance,  RT  (-Rr)  and  angle  of  incidence/reflection  <t>  See  Fig. 

6 

The  station  to  trail  distance  is  given  by 

*r  "  J(/i+  L*/8R,)2+  L*/4  (2.2-1) 


And  the  incidence/reflection  angle  <$>  is  found  thru 
secJ(6)  -  1  ♦ 


(2.2-2) 
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where , 

h  trail  altitude  (100  km). 

L  great  circle  distance  between  stations. 

Rr  distance  from  transmitter  to  trail. 

Rx  distance  from  receiver  to  trail. 

Rt  radius  of  the  earth,  6400  km. 


The  derivation  of  the  above  is  done  in  section.  4.1 
Since  it  is  often  convenient  to  express  L  and  the  resulting  Rg  in  thousands  of 
kilometers  for  ease  of  computation  we  then  have: 


tfr-  /?,- V(  1  +  ( .01 953) I2)2 *  L2 / 4 


(2.2-3) 


sec2(4>) 


1  + 


L2 

(0.2  +  ( .039)I2)2 


(2.2-4) 


s.t.  for  L-1000  km.  we  use  Lr-1  to  quickly  yield  Rj— Rr  “  514  km.  and  sec2  <}>-  18.5 


2.3  Trail  Arrival 

In  this  section  we  are  concerned  with  quantifying  the  trail  arrival  process 
in  terms  of  communication  system  parameters.  Since  the  arrival  is  Poisson  we 
need  to  estimate  the  average  time  between  bursts  or  equivalently  the  number  of 
bursts  per  second.  Oetting,  [9]  summarized  the  results  of  the  COMET  system  showing 
an  average  time  between  bursts  of  4  to  20  seconds. 

It  is  often  neccessary  to  estimate  the  number  of  arrivals  of  bursts  per 
second  with  an  electron  line  density  of  q  between  some  qj  and  q2.  We  let  A' (q) 
be  the  number  of  trails/s  with  an  electron  line  density  between  q  and  q+dq  e/m. 
Using  Weitzen  [10]  we  write 
19.1667-  ♦sin*- (tf 


A(q)- 


101 


(2.3-1) 


where  i; is  the  idealized  beamwidth  of  the  antenna  and  all  other  parameters  have 
been  previously  defined.  The  number  of  trails/s  with  q  in  a  given  range  [q1,q2], 
A,  is  found  by  integrating  the  above  density  over  q  yielding  the  following 
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A  19.1667-»8in»(*«»1f1  <7.1 

9i  L 

A  19.1667-»sin»(*,02r  q ,  ~j 

<?i  L  qj 

note  that  for  overdense  trails  (qi  >  101*)  we  scale  the  constant  19.1667  in  the 
above  Eq .  to  be  3 . 

In  addition,  care  must  be  exercised  in  setting  qmin  since  this  is  an  empirical 
relationship.  For  example,  for  a  1000  km  link  with  antenna  angle  of  46°  qnin  must 
be  above  1013  e/m  for  the  formula  above  to  yield  reasonable  rates. 

Finally,  we  mention  again  the  dependence  of  the  number  of  bursts/s,  A,  on 
other  factors  such  as  global  location,  time  of  day,  and  season.  These  factors 
should  appropriately  scale  any  value  of  A  used  for  analysis. 

2.4  Power  Spectral  Density  of  the  Received  Noise 

For  signal  frequencies  appropriate  to  the  MBC  the  dominant  noise  sources 
are  galactic  and  man-made  in  origin.  Consequently,  the  received  noise  power 
spectral  density,  N0,  is  a  function  of  the  galactic  noise  picked  up  by  the 
receiver  antenna  and  the  receiver  thermal  noise. 

N0  is  modeled  as  follows : 


(2.3-2) 

(2.3-2) 


N0 


104(  A  y 3 
L t  l  15/ 


(2.4-1) 


where , 

k-  Boltzmann  constant.  1 .3805  •  1  O'23  J/K 
T0-  290°  K 

Lt -  power  loss  between  the  antenna  and  receiver 
F  -  receiver  noise  figure 
X.-  wavelength  in  meters 


We  choose  typical  values  for  Lr  and  F  as  1.3  and  2.5,  respectively  to  yield 
the  following  model  for  N0 

>2  3 


Nr 


4-  10' 


W/Hz 


(2.4-2) 
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2.5  Modulation  and  Bit  Error  Rate 

The  bit  error  rate,  BER,  denoted  here  as  Pb,  is  In  general  related  to  the 
bit  energy  E*,  and  power  spectral  density  of  the  received  noise  N0  as 

(2.5-1) 

where  g()is  determined  by  the  modulation  technique. 

Since  E^t)-  P(t)/R(t)  where  R(t)  the  bit  transmission  rate  may  be  time  varying 
or  fixed  we  have  from  Eq .  (2.5-1) 

f‘(°-9(S))  (2s-2) 

or  alterntively ,  we  can  solve  for  the  power 

PtO-No  •*(/)•  g''(P*(/))  (2.5-3) 

where 

&"'(•)  is  the  inverse  of  g(  ). 


For  BPSK  we  have: 


(2.5-4) 


(2.5-5) 


(2.5-6) 


and 


P(O«N0  ^(,)  ln(~^j) 


(2.5-7) 
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3  ANALYSIS  AND  RESULTS 

This  chapter  presents  and  explain  our  approach  and  reasoning  to  the  use 
and  analysis  of  the  MB  channel.  The  significant  analytical  expression  to  be 
used  as  design  and  analysis  tools  are  given. 

3.1  Introduction  and  Approach 

As  stated  earlier,  our  aim  in  the  analysis  is  to  arrive  at  a  set  of 
analytical  relationships  for  important  communication  parameters  such  as  burst 
duration  and  throughput  using  as  accurate  a  model  as  possible.  These  relationships 
will  serve  as  analysis  and  design  tools  where  the  effects  of  system  parameters 
and  constraints  (power,  geometry,  BER)  on  throughput  can  be  assessed.  Our  goal 
is  to  estimate  the  optimal  average  throughput  for  a  given  communication  system. 

Two  types  of  communication  systems  are  considered: 

1)  A  system  with  a  constant  transmission  rate  and  a  time  varying  bit  error 
rate  which  is  constrained  not  to  exceed  some  maximum  allowable  value. 

2)  A  system  employing  a  variable  transmission  rate  which  mimics  the  time 
varying  power  in  order  to  maintain  a  practically  constant  bit  energy  and  BER 
for  all  transmitted  bits. 

For  the  constant  bit  rate  system,  since  the  power  is  time  varying,  the 
probability  of  bit  error  will  be  time  varying  in  a  fashion  dictated  by  the 
relevant  modulation  function.  Specifically,  for  the  underdense  burst  we  have 
an  exponentially  decaying  power  in  time  and  consequently  the  BER  is  monotonically 
increasing.  Whether  underdense  or  overdense  we  see  the  need  then  to  impose  a 
ceiling  on  the  BER.  For  such  systems,  where  the  probability  of  bit  error  is 
time  varying,  we  specify  a  maximum  allowable  BER,  P]„„,  such  that  all  tranmitted 

bits  shall  have  a  BER  less  than  or  equal  to  Pi _  This  constraint  on  the 

probability  of  bit  error  for  constant  transmission  rate  systems  implies  the 
existence  of  a  minimum  power  level  Pain  such  that  all  transmitted  bits  have  a 
received  power  greater  than  or  equal  to  PBin.  Such  is  obvious  when  remembering 
that  the  modulation  function  relating  power  to  BER  is  an  increasing  function  for 
decreasing  argument.  In  short,  Pb(t)  -  P^  when  P(t)  -  P^.  The  time  for  which 
the  power  received  from  the  meteor  burst  exceeds  a  prescribed  threshold  (Poln) 
is  the  burst  duration  tB  (also  duty  cycle  [1]). 

Clearly  the  duration  of  time  for  which  the  received  power  P  is  greater 
than  Poin  varies  from  one  meteor  burst  to  the  next.  For  a  constant  bit  rate 
system  the  burst  duration  is  crucial  in  determining  the  number  of  bits  transmitted 
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during  a  given  burst  and  ultimately  in  estimating  throughput.  We  shall  therefore 
analyse  its  relationship  to  system  parameters,  derive  its  statistics  and  compare 
it  with  current  dogma. 

As  we  increase  our  set  constant  bit  rate  the  corresponding  has  to  be 
increased  resulting  in  shortening  the  effective  duration  of  transmission.  This 
clear  tradeoff  between  transmission  rate  and  duration  of  transmission,  tB, 
suggests  an  optimal  choice  for  transmission  rate  (or  equivalently  an  optimal  tB) 
so  as  to  maximize  the  throughput.  The  link  between  the  number  of  bits  transmitted 
per  burst  and  throughput  is  done  through  the  average  arrival  rate  statistics  and 
trail  shape  statistics  (variation  in  q  and/or  B) .  The  judicious  choice  of 
transmission  rate  for  the  given  system  provides  us,  therefore,  with  the  optimal 
throughput . 

The  second  type  of  communication  system  considered  here  is  a  variable 
transmission  rate  system.  Under  such  a  system  the  bit  rate  mimics  the  time 
behavior  of  the  received  power  (or  equivalently  stated  the  bit  duration  varies 
in  a  reciprocal  manner  to  the  received  power  variation)  so  as  to  maintain  a 
constant  bit  energy.  The  direct  outcome  of  such  approach  is  a  constant  probability 
of  bit  error.  More  importantly,  we  note  that  we  are  not  restricted  by  the  need 
to  operate  above  a  prescribed  threshold  of  received  power;  for  as  the  power 
diminishes  we  simply  compensate  by  increasing  the  bit  duration.  This  adaptive 
approach  implies,  theoretically  at  least,  for  the  underdense  burst  an  infinitely 
long  time  of  transmission.  This  of  course,  is  offset  by  the  practical  limitation 
on  bit  duration  and  the  need  to  maintain  a  fairly  constant  BER.  The  assumption 
here  is  that  the  hardware  is  fast  enough  to  vary  the  transmission  rate  (bit 
duration)  from  one  bit  to  the  next.  Since  the  power  in  most  cases  diminishes 
during  the  duration  of  the  bit,  for  long  bit  duration  we  deviate  from  the 
requirement  to  maintain  a  constant  bit  energy  and  constant  BER.  It  is  therefore 
necessary  to  take  this  fact  into  account  in  determining  the  duration  of 
transmission.  This  criterion  of  allowing  only  a  small  drop  in  the  assumed 
constant  bit  energy  imposes  a  restriction  on  the  maximum  bit  duration  which  in 
turn  limits  the  duration  of  transmission.  As  in  the  previous  case  we  estimate 
the  throughput  by  averaging  over  the  ensamble  of  trails  and  incorporating  the 
effects  of  trail  arrival  statistics.  Finally,  we  express  the  correct  packet 
probability  in  terms  of  the  modulation  function,  packet  length  and  communication 
system  parameters. 

3.1.1  Protocol 

Both  stations  have  receiving  and  transmitting  capabilities.  The  station 
which  acts  as  transmitter  of  data  listens  continuously  for  a  continuous  tone 


26 


which  is  being  sent  from  the  receiver.  Upon  detection  of  this  tone  both  the 
presence  and  the  strength  of  the  channel  is  known  to  the  transmitter  which 
immediately  commences  transmission.  The  situation  is  completely  symmetrical 
with  respect  to  the  receiver.  Clearly,  the  beginning  of  the  burst,  i.e. 
availability  of  the  channel,  is  known  to  the  transmitter  within  a  neglibile 
propagation  delay  (less  than  3.5  msec  for  a  1000  km  link).  The  continuous  tone 
used  for  monitoring  the  channel  is  at  a  different  frequency  than  the  data  for 
ease  of  gating  and  synchronizing.  For  constant  bit  rate  systems  the  essential 
feature  of  this  scheme  is  the  ability  of  the  transmitter  to  discern  the  beginning 
of  the  burst.  For  variable  bit  rate  systems  the  strength  of  the  channel  as 
gauged  by  the  transmitter  in  addition  to  the  burst's  start  are  important  for 
adaptive  bit  rate  determination. 

The  above  description  is  the  rudimentary  layer  upon  which  additions  and 
refinements  could  be  made. 


3.2  Improvements 

The  advantages  of  the  analysis  herein  stem  basically  from  considering  a 
more  complex  MBC  model  than  found  in  the  literature.  Most  researchers  [9,  21] 
while  acknowledging  the  significant  contribution  to  communication  of  the  rare 
yet  robust  overdense  burst  have  tended  to  concentrate  on  the  underdense  burst 
phenomenon  because  its  simpler  exponential  form  is  more  mathematically  tractable. 
Able  [11] ,  however,  does  analyze  the  overdense  trail  by  a  piecewise  linear  (first 
order)  approximation  resulting  in  conservative  channel  duration  estimation.  We 
choose  a  second  order  approximation  to  improve  the  fit  of  the  overdense  tran¬ 
scendental  power  time  function.  This  yields  better  results  for  channel  duration 
for  the  overdense  case.  For  the  exponential  case,  researchers  to  date  have 
assumed  a  fixed  decay  time  constant.  Hampton  [21]  in  his  analysis  of  an  MB 
broadcast  system  used  a  constant  decay  rate  which  was  the  conditional  expected 
time  constant  given  a  fixed  burst  duration.  This  approach,  however,  implies 
that  the  decay  time  constant  is  dependent  upon  an  average  of  a  function  of  the 
electron  line  density  whereas  it  is  known  that  the  time  constant  is  independent 
of  the  electron  line  density.  To  enhance  the  accuracy  of  our  model  we  chose  to 
adhere  to  empirical  results  and  assume  the  decay  time  constant  to  be  a  random 
quantity  independent  of  the  electron  line  density. 

Some  researchers  [34]  have  further  simplified  the  MBC  model  by  assuming  a 
fixed  electron  line  density  for  all  exponential  bursts.  This  is  contrary  to 
experimental  observation.  Here,  as  in  the  case  of  the  exponential  decay  time 
constant,  we  incorporate  this  random  parameter  into  our  model. 
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In  analyzing  throughput  we  compare  two  systems:  constant  bit  rate  system 
and  variable  bit  rate  system  and  derive  the  average  throughput  for  both.  For 
constant  bit  rate  we  derive  the  expression  for  the  optimal  average  throughput. 
This  optimal  average  throughput  is  achieved  by  finding  the  "best"  bit  rate  to 
be  used  for  the  communication  system  so  as  to  maximize  throughput  for  the 
underdense  burst.  This  is  not  to  be  confused  with  the  optimal  bit  rate  for  a 
given  burst  that  maximizes  the  number  of  transmitted  bits  for  the  burst  under 
consideration. 

In  summary  we  note  that  the  bedrock  of  a  communication  system  using  the 
MBC  is  its  model.  To  enhance  the  analysis  we  used  a  model  with  as  few  sim¬ 
plifications  as  possible  while  maintaining  mathematical  tractability . 

3.3  Sample  System 

For  quantitative  appreciation  of  the  results  we  need  to  assume  some  values 
for  a  communication  system.  The  set  of  these  assumed  parameters  values  will 
comprise  what  will  heretofore  be  refered  to  as  our  sample  system.  This  sample 
communication  system  represents  practical  or  average  values.  Chapter  2  contains 
the  definition  of  the  parameters. 


L  -  1000  km  link  distance.  From  section  2.2 

we  have  for  the  next  two  parameters : 
sec($)  -  18.5  and 


*r-R. 

Pt 

Gt 


Gr 


-514  km. 

-  1000  w 

-  10  dB. 

-  10  dB.  And  for  antenna  beanrwidth  we  use 
-  45.84°  (.8  rad) 


X.  -6  m.  (50  MHz) 

D  -  8  m2/s 

Pboax  “  Max.  BER  for  constant  bit 

rate  system. 

Using  BPSK  modulation  we  get: 

g-HPbo-x)  -  9 
W,  -  4  •  10-21|^80^)23-2.5 


W/Hz 
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See  discussion  above  Eq.  (2.4-2). 

N0  -  4.89  10-2°  at  50  MHz. 

For  comparison  the  BER  for  variable  rate  system,  Pb, 
is  set  equal  to  Pi _ 

R^y  -  1000  b/s  (minimum  variable  rate  allowed) 


3.4  Constant  Bit  Rate  and  Constrained  BEJL  System 

3.4.1  Channel  Duration 

Channel  duration  or  burst  duration,  tB  is  the  time  for  which  the  power 
received  from  the  burst  exceeds  a  given  minimum  value  ^min  •  1 1  i  S  also  the 
duration  of  data  transmissions  for  a  given  burst  where  all  bits  possess  a 
probability  of  error  less  than  some  maximum  allowable  BER,  Ptaax.  P^  is  determined 
thru  the  modulation  function  by  the  following  parameters:  Transmission  rate 
(non- time -varying) ,  power  spectral  density  of  the  received  noise,  and  the  maximum 
allowable  bit  error  rate.  The  burst  duration  tB  is  determined  by  and  Cu 

which  incorporates  system  parameters  such  as  transmitter  power,  link  geometry 
effects,  carrier  wavelength  etc. 

3. 4. 1.1  Burst  Duration  Statistics  for 

Constant  Bit  Rate  --  Pnderdense 

The  channel  duration,  tB,  -  burst  duration  -  is  the  time  of  transmission 
for  a  given  burst  i.e.  the  usable  portion  of  the  burst,  assuming  a  constant 
transmission  rate,  R,  some  given  maximum  allowed  BER,  and  power  spectral  density, 
N0  ,  see  Fig.  11. 

For  such  conditions  for  the  underdense  burst  we  have: 


where 


and 


(3.4.1  .1-1) 


2.5179581  10 


Wo*g'V,MI) 

_  _  _  %  a  (  33  3S»  \ 

PTCrC,\  exp[-kWJ 


*T 


or  equivalently 
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(  P  bn  ax  ) 


'  R 


(3.4.1 .1  -2) 


We  now  average  over  the  ensemble  of  underdense  bursts  (averaging  over  q  and  b) 
to  yield  the  average  burst  duration: 

The  Average  Burst  Duration--  Pnderderdense : 


-  2  B 

t.  -  ■ 


p  -  1  1  -  x 


[1  -x(l-ln.v)] 


(3.4.1  .1-3) 


where 


f  Qmln  ^ 

l  <7(/ ; 


y," 

'N0-Rg\Pbma,y 

j 

CuRl 

and 


B- 3.16610 


. 3  \2sec2$ 


Ru  “  10'4  p  -  1  .6 


r,  R  _  rWmll» 

K  N  p  v.  ) 


^  ofif  (^dma*) 


For  a  given  communication  system  CU(  J? ,  tf0,  s  and  gr()are  specified.  RN 

is  the  normalized  (  with  respect  to  the  highest  possible  rate  for  an  underdense 
burst)  bit  rate. 

Fig.  14  plots  the  above  the  above  for  our  sample  system  (Sec.  3.3)  for  various 
link  distance  as  a  function  of  R^.  (R^R/R^  .R^x-86.5  kb/s) 

The  vay  Unfit  2f  il  given  hll 


1 


(p-  l)2  (1  -x)2 


,«*l 


where 

-  2  -  x(3  +  y2)  +  x2(l+y2) 


(3.4.1 .1-4) 
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c2  -  1  -  x(2+ln2(x))  +  x2 
y  -  1  -  ln(x) 

x  from  previous  equation 


and 


For  B  Rayleigh  Jfc  -  -  -  1  ,  For  B  exponential  k  -  1 


The  nth  moment  of  tB  Is  given  by : 


-  2*B*  n' 


(P-1)'  I"* 


,ln'(x) 


l-I(-l)'^ 

1-0  1 


The  PDF  and  CDF  of  tB  are: 


^i.Oi  )“<?'(  1  ~  xF,(bmin)- f^exp^-£^J/t(b)db 


(3.4. 1  . 1  -5) 


(3.4.1  .1  6) 


(3.4.1  .1-7) 


where 

x  and  p  are  from  above  eq.  (3.4. 1.1-3) 

P'  ~(P  -D/2 

^min  “  ^1  /  (  (l/p')  ln(l/x)  ] 

Q'  -  l/(l-x) 

See  Fig.  15  for  a  plot  of  Eq.  (3. 4. 1.1-6)  (where  B  is  assumed  Rayleigh), 
and  an  exponentially  disributed  burst  duration  is  drawn  for  comparison. 
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For  the  case  of  x«l  l.e.  ?^n«  Cy  q02  .  essentially  Ignoring  (considering 
all  physical  trails)  we  have  for  the  solution  of  the  above  integrals (  boln  -  0 

): 


B  assumed  exponential: 

/,f(/,)-2£-*0(z)  (3.4. 1.1-8) 

D 

and 

1  “  2*,(z)  '>-°  (3. 4. 1.1-9) 

wherez  *•  ,  /  4t=-  ■  t , 

V  s 

where,  K0  and  iq  are  the  modified  Hankel  functions  of  order  0  and  1 

respectively . 

B  assumed  Rayleigh: 

F ,  (f,)“  1  “2  f  bexp(-b2  -  I  ,/b)db  (3.4.1.1-10) 

*  Jo 

“  1  -27,(1,) 

where  f1(t1)  is  a  tabulated  function  [35].  See  Fig.  9. 

3. 4. 1.2  Burst  Duration  Statistics  for 

Constant  Bit  Rate  is,  Over  dense 

An  exact  solution  to  the  roots  of  the  power  received  equation  -  overdense 
-  is  not  known.  The  difference  between  these  time  instances  yields  the  duration 
of  the  channel.  We  can  approximate  the  power  function  in  a  piecewise  linear 
fashion  (comprised  of  two  straight  lines)  or  alternatively  use  a  second  order 
approximation.  The  second  order  approximation  used  here  improves  the  accuracy 
of  the  results. 
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The  results  ere  stated  without  the  lengthy  derivation.  For  derivation  see 
section  4. 3 


1 1  «  .95 aq 


9«...  10 “e/m 


where 

a  and  C0  are  from  Eq  (2. 1.5-2) 
qr  ~  (e/a)  (Pmin/C0  )2  -  qyo  RN02 
RN0  -  R/Roaxo  •  normalized  bit  rate. 

Rmaxo  -  C0[a  quo/e]  VNog'KPbcax) 

maximum  possible  bit  rate  for  any  overdense 
burst. 

q»ino  -  “ax  IqT.qLo)  qLO-lO1*  quo  -  1016  e/m. 
e  is  the  natural  logarithm  base. 


(3.4.1 .2-  1 ) 


Note  that  q,^  is  the  minimum  (physical)  electron  line  density  that  can  support 
an  overdense  burst.  qT  is  specified  by  our  choice  of  system  parameters  and  bit 
rate  and  may  not  correspond  to  a  real  level  of  overdense  electron  line  density. 
See  discussion  following  Eq.  (3. 4. 1.2-3).  See  Fig.  12.  but  note  that  in  the 
figure  qT<qLo  • 

We  can  also  relate  qT  to  q^  recalling  the  definition  of  Pmln  -  C„  q2min  ; 


(3.4.1  .2-2) 


where  the  parameters  Cu,  C0,  a  and  q^  are  defined  in  Eqs .  (2.1. 1-2)  and  (2. 1.5-2) 
and  (3, 4. 1-1)  respectively. 

Average  Burst  Duration  -  Overdense: 

The  average  burst  duration  for  the  overdense  case  is  found  by  averaging  over 
the  ensemble  of  bursts,  yielding 


where , 


(3.4.1  .2-3) 
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Q  mu*o  -  max{qT.ql0} 
Qr  m  Quo'  &  no 


C0Ja ■ q t 


■fe -N0- g'\Ph. 


) 


11,*,  is  the  normalized  bit  rate  with  respect  to  the  maximum  rate  possible  for  any 
overdense  burst,  Rn^o.  Note  that  the  above  can  be  expressed  as  two  functions 
one  for  qT  <  (or  equivalently  RN0  <  .1)  where  qalno  -  qj,0  and  the  other  for 
qT  >  ^lo  (°r  equivalently  R„0  >  .1)  where  q,,^  -  qr.  In  addition,  the  value  of 
Q  varies  with  the  domain  as  well  since  Q  is  a  function  of  qmino  -  see  Eq.  (2. 1.6-2). 
For  the  case  where  R^q  >  .1  (qT  >  q^)  we  have 
Q  “  -5(quo)  5  Rho  /  (l-R«o) 

For  the  case  where  R^  <  .1  (qT  >  q^)  we  have 
Q  “  .5(qL0)-s  / •  9 

For  many  practical  systems  we  generally  have 

qT<<clL0<<clU0  i-e-  laino  “  RlO  •  Q”(qLo)1/Z/2 
resulting  in: 

t„~  .95a  •  -JqvoQio  (3. 4. 1.2-4) 


and  for  qT  <  .01  q^  we  have  for  the  burst  duration  variance: 

<--95y^gi;2-.9Syl031  (3.4.1  .2-5) 

For  our  sample  system  definsd  in  Sec.  3.3  the  standard  deviation  of  tB  -  11  sec. 

Burst  Duration  pdf  and  cdf  -  Overdense: 
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CDF  of  CB: 


f° 

~  ^  l  mm 

1  <?' 

1-f— - 

V' 

1 

_ 

) 

where 

p— 1 . 5  and 


r.  9? 


^  1  min  m  ■9Saqnl,0^j  1 


T.,9S2|I.,95|(^ 

PDF  of  tB: 


max<.qT.qL0) 


P-  1 


(1*2)' 


■H] 


where 


and  all  other  parameters  are  defined  by  previous  equations. 


(3. 4. 1.2-6) 


(3.4.1  .2-7) 


3.4.2  Bits  £££  Gl^en  Burst 

Bits  per  given  burst  under  constant  bit  rate  and  constrained  BER  is  the 
number  of  bits  transmitted  for  the  burst  whose  q  and/or  B  are  known.  Since  the 
bit  rate,  R,  is  assumed  constant  the  number  of  bits,  NB,  is  simply  the  product 
of  R  times  the  burst  duration  tB  as  defined  in  previous  sections.  Clearly,  the 
level  of  power  P..t„  which  reciprocally  affects  tB  is  dictated  by  the  modulation 
function  and  in  general  increases  with  increasing  bit  rate  R.  As  a  result  we 
see  the  trade-off  between  the  bit  rate,  R,  and  transmission  duration,  tB,  such 
that  increasing  one  decreases  the  other.  This  suggests  the  existence  of  an 
optimal  bit  rate,  R*B,  (or  equivalently,  optimal  transmission  duration)  that 
maximizes  the  number  of  bits  transmitted  during  the  given  burst.  Clearly,  using 
any  bit  rate  other  than  the  optimal  R*B  will  result  in  under-utilization  of  the 
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burst.  Alternatively  stated,  given  a  bit  rate  one  should  transmit  only  for  as 
long  as  the  optimal  burst  duration  dictates .  The  optimal  bit  rate  for  the  given 
burst  is  then  used  to  calculate  the  maximum  number  of  bits,  NB*  ,  that  can  be 
transmitted  during  this  burst.  A  ratio  is  then  obtained  for  the  same  burst 
between  the  number  of  bits  transmitted  using  variable  bit  rate  (BER  constant  - 
P, „„.T)  and  the  NB*  using  constant  bit  rate.  This  ratio  has  been  shown  by  Abel 
[11]  to  be  around  2.5.  Using  the  "best"  (highest  q)  burst  for  overdense  and 
underdense  trails  bounds  can  be  found  for  the  maximum  possible  NB*  (constant  rate 
system)  and  maximum  transmitted  bits  using  variable  rate. 

Two  points  however,  must  be  kept  in  focus: 

1)  Since  R*  varies  from  trail  to  trail  and  a  priori  knowledge  of  the  trail 
behavior  is  impossible  to  obtain,  the  optimal  (constant)  bit  rate  for  the  system, 
i.e.  for  all  possible  bursts,  is  still  unknown.  Parenthetically ,  we  note  that 
attempting  to  find  the  burst's  shape  from  knowledge  of  its  beginning  is  both 
wasteful  of  valuable  transmission  time  and  unreliable  since  many  bursts  deviate 
from  ideal  behavior.  As  will  be  shown,  the  important  parameter  for  throughput 
determination  is  not  NB  --  the  number  of  bits  for  a  given  burst  --  but  rather 
the  number  of  bits  per  burst  averaged  over  the  ensemble  of  all  bursts  since  the 
latter  is  directly  proportional  to  the  average  throughput. 

2)  While  it  is  true  that  more  bits  can  be  transmitted  using  variable  bit  rate 
during  the  burst  rather  than  employing  constant  bit  rate,  the  situation  is 
somewhat  analogous  to  an  apples/oranges  comparison.  For  the  underdense  burst 
and  constant  bit  rate  all  bits  except  the  last  one  will  possess  a  BER  smaller 
than  Pv,TT  --  the  worst  BER,  whereas  for  variable  bit  rate  all  bits  will  possess 

the  same  BER  -  P, _  Consequently,  the  probabiliiity  of  having  an  error-free 

packet  of  data  is  higher  for  constant  bit  rate  than  for  variable  bit  rate. 

We  present  the  analysis  for  the  given  underdense  burst  by  first  considering 
constant  bit  rate  with  constrained  BER  and  then  considering  variable  bit  rate 
system  with  constant  BER. 

Pennant  Sit  rata  mta 

Given  an  underdense  burst  -  q  and  B  are  assumed  known.  PhriT  is  specified  for 
the  system  and  R  is  the  constant  bit  rate.  Let  NBbe  the  number  of  bits  transmitted 
on  this  burst.  All  bits  have  BER<PbMX.  All  bits  have  power  P>PBln.  From  the 
discussion  on  burst  duration  and  the  above  definition  of  NB  we  have: 

(3. 4. 2-1) 

Using  eq.  (3. 4. 1.1 -2)  for  tB  we  rewrite  the  above  as 
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N,- R  -  B  \n(R0/R)  (3. 4. 2-2) 

where  R<,  -  Cu  q2  /  N0  g'KPbmax) 

Note  that  R,  Is  a  constant  set  by  the  choice  of  system  parameters  and  physical 
characteristics  (1)  of  the  given  burst.  From  the  above  it  is  clear  that  for  a 
given  burst  a  different  choice  of  bit  rate  R  yields  different  number  of  transmitted 
bits,  Nb.  To  find  the  bit  rate  RB  that  maximizes  NB  we  simply  set  its  derivative 
with  respect  to  R  to  zero  and  solve 
dN  . 

—~-Q-\n(R<1/R)-\  (3. 4. 2-3) 

CZ  a 

or  R*B  -  R0/e  where  e  is  the  natural  logarithm  base.  Using  the  last  result  in 
eq  (3.4. 2-2)  yield  NB*  the  maximum  number  of  bits  for  this  burst. 

.  R0- B 

(3. 4. 2-4) 

where  -  Cu  q2  /  N0  g'^Pbmax) 

From  the  above  we  see  that  P*^,,  the  minimum  power  corresponding  to  R*B  is  equal 
to  P(0)/e.  This  implies  that  the  burst  duration  (transmission  time)  i.e.  the 
time  it  takes  to  drop  to  P*^,„  from  P(0)  is  equal  to  one  decay  time  constant  B. 
Equivalently,  then,  the  above  analysis  states  that  for  a  given  bit  rate  R  the 
optimal  time  of  transmission  T*B  must  be  set  to  B  to  yield  the  above  maximum 
transmitted  bit  N*B.  This  result  has  been  demonstrated  by  Abies’  work  [11]  as 
well . 

Yayjgble  Ml  Ra.t.6  Ssa&sa 

For  the  given  underdense  burst  we  have  Pb(t)-Pb«ii*x  f°r  bits.  Under  such 
conditions  the  variable  bit  rate  R(t)  is  given  by 

R(t)  -  Roe'6/8  ((3. 4. 2-5) 

where  R<,  -  Cuq2/g*1(PbCMlx) 

Note  that  the  power  to  transmission  rate  ratio  which  specifies  the  bit  energy 
is  constant  and  hence  the  BER  is  maintained  constant. 


Since  the  burst  is  assumed  to  be  known  a-priori  we  can  compute  the 
bit  duration  (inverse  of  bit  rate)  for  each  bit  exactly  such  that  the  integral 
of  power  received  over  the  individual  bit  duration  i.e.  the  bit  energy  is  the 
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sane  for  all  transmitted  bits.  See  Fig.  13. 

The  number  of  bits  for  variable  rate  is  given  by 

.V.„-  f* R(t)dt- *„•  B  (3. 4. 2-6) 

Jo 

Comparison  of  the  Two  Systems 

We  nov  compare  constant  bit  rate  system  using  the  optimal  R  (with  BER 
<  Pbo,^)  for  the  particular  burst  versus  variable  bit  rate  which  mimics  the  bursts 
power  behavior  (such  that  all  bits  possess  the  same  BER,  Pbm«x)  •  For  each  of  the 
above  we  found  the  number  of  bits  transmitted  N*B  and  NBV  respectively.  Taking 
the  ratio: 


Nbv/N*b  “  e 

We  see  an  improvement  of  2 . 7  for  the  utilization  of  variable  bit  rate 
over  constant  bit  rate.  We  stress  however  that  the  necessary  priori  knowledge 
(or  estimation)  of  burst  characteristics  is  unreliable.  In  addition,  although 
variable  bit  rate  results  in  a  higher  number  of  transmitted  bits  for  a  given 
burst  versus  constant  bit  rate  there  is  a  price;  all  bits  have  a  BER  —  P^.,  for 

variable  bit  rate  system  whereas  for  the  constant  bit  rate  system  P, _ is  a  worst 

case  (last  bit  BER) .  This  fact  portends  worse  probability  of  correct  packet 
transmission  for  variable  rate  system  than  constant  rate  system  under  the  above 
conditions . 

3.4.3  Throughput  f'..*  Constant  Bit  Rate 

We  are  given  a  communication  system  with  a  specified  maximum  allowed  BER 
P>_„.  constant  transmission  rate  R,  received  noise  power  spectral  density  N0  and 
general  system  parameters  (link  distance,  frequency,  etc.)  as  defined  by  Ca  or 
Co- 
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To  express  Che  throughput  T  of  the  system  we  introduce  NB1  as  the  number 
of  bits  transmitted  on  the  1th  burst  and  Mx  as  the  number  of  bursts  occurring 
in  the  period  of  x  seconds.  Since  the  bursts  are  independent  of  each  other,  and 
assuming  no  overlap  of  bursts ,  we  have : 

1  ^ 

(3. 4. .3-1) 

T  i- 1 

•AN, 

where  A  is  the  average  number  of  bursts/s  and  A/ais  the  expected  number  of  bits 
per  burst  as  averaged  over  the  ensemble  of  all  bursts'  profiles.  (The  subscripts 
U  or  0  are  added  to  distinguish  between  underdense  and  overdense  cases.) 

For  underdense: 

N.o~  J  J  N,{q.B)f9{q)f,{b)dqdb  (3. 4. 3-2) 

where  fT  ()  is  the  pdf  of  x. 

For  overdense: 

N»o-/wf(q)/,(q)dq  (3.4. 3-3) 


Note  that  the  above  three  equations  are  dependent  on  the  choice  of  bit  rate  R. 

3. 4. 3.1  Optimal  Ay?ri>gft  Thlfiu&frpyt  far  Constant 
Bit  Rate  System  --  Pnderdense . 

In  this  section  the  optimal  average  system  throughput  for  the  underdense 
burst  is  given.  This  is  the  maximum  average  number  of  bits  per  unit  time  for 
the  'best'  transmission  rate  R  (and  other  system  constants)  assuming  only 
underdense  bursts  are  present. 

This  optimal  throughput  is  found  by  first  deriving  the  dependence  of  the  average 
throughput  (Eq.  3.4. 3-1)  on  the  bit  rate  R,  and  solving  for  the  optimal  bit  rate 
R*  that  would  yield  the  desired  optimal  throughput  T*.  The  solution  is  done  with 
respect  to  the  normalized  bit  rate,  R#,  which  yields  a  result  independent  of  a 
particular  choice  of  system  values. 


For  the  given  system  P, - N0,  R  and  Cu  are  specified.  We  have  defined  R* 

the  normalized  bit  rate  as: 


R  /  N 

R  mxx  Qu 


where 


R 


max 


C  t/Qv 

Wo9'‘(  Porno*) 


(3. 4. 3.1  -  1) 


the  maximum  bit  rate  achievable  in  a  given  system  using  underdense  bursts. 

Omitting  the  cumbersome  derivation  we  present  the  following  results:  (for 
the  relevant  derivation  see  Sec.  4.2) 

The  average  number  of  bits  per  underdense  burst  N  tu - R  ~t  can  be  expressed  in 
terms  of  the  normalized  bit  rate  using  Eq.  (3. 4. 1.1-3)  as: 


—  2  B 


1  -R% 


(3.4.3. 1  -  2) 


where 

p'-(p  -  l)/2  ,  p-1.6  and  R„.T  is  defined  above. 


Note  that  in  the  above  expression  the  bursts  considered  were  those  with  an 
electron  line  density  q  between  q^  (as  defined  by  Pnin)  and  qu  (the  physical 
limit).  The  arrival  rate  for  such  condition  is  given  by  Eq.  (2.3-2)  where  q! 
is  set  to  q,,^,,  and  q2  is  qy.  Together  with  Eq.  (3.4. 3. 1-1)  we  have  for  the  arrival 
rate  of  underdense  bursts : 


where, 

,,  19.1667  •  tsint(£«Oz 

A  * - 

Qu 

Qu  *  1 0 14 


(3. 4. 3.1  -3) 


The  expression  for  throughput  using  underdense  bursts  is  given  by  the  product 
of  the  last  two  Eqs.  as  specified  by  Eq.  (3. 4. 3-1).  The  result  is  a  function 
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of  the  normalized  bit  rate  which  is  just  a  scaled  bit  rate.  In  the  product  we 
substitute  for  R,..  from  Eq.  (3.4. 3. 1-1)  and  for  Ct,  and  B  from  Eq.  (2. 1.1-2). 
For  a  given  set  of  system  parameters  we  thus  have  the  throughput  (using  underdense 
bursts)  for  any  choice  of  (normalized)  bit  rate  as: 


,  ,  /  33  339  \ 

-i9  prGTG,^$sin$sec2*exp[-;^7r;J- X 

(p  -  1  )N0  ■  g'l(P inax)  ■  D  -  Rt 
■J  ft  N  ~  R  N 


\-n 


L[1  -K  !"(*„)] 


(3.4.2. 1  -3) 


where 


P 


K 

Rt  and  4> 
Pi 

gt.Gr 

k 

D 

No 

6'HP^) 


is  an  empirical  constant  .  p'-(p-l)/2 
(taken  here  as  1.6,  Section  2.1.2) 
idealized  antenna  beamwidth. 
geometric  parameters  detrained 
by  link  distance  (from  section  2.2) 
transmitter  power. 

transmitter  and  receiver  antenna  gains, 
carrier  wavelength. 

diffusion  coefficient  of  the  atmosphere, 
power  spectral  density  of  the  noise. 
(Sec.  2.4) 

ratio  of  minimum  bit  energy  to  N0. 


Note  that  the  equation  is  written  in  two  parts;  the  first  incorporates  all  the 
parameters  of  the  particular  system  and  the  second  part  reflects  the  variation 
in  throughput  with  choice  of  bit  rate.  The  above  throughput  (normalized  with 
respect  to  its  maximum)  for  our  sample  system  (defined  in  Sec.  3.3)  is  plotted 
in  Fig.  16. 

We  would  like  to  find  from  the  above  equation  the  optimal  normalized  bit 
rate  R**  which  yields  the  maximum  throughput  T,,*  .  This  optimal  normalized  bit 
rate  for  any  given  system  is  found  numerically  to  be  (see  Fig.  16): 
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R**  -  -07 

or  equivalently  the  optimal  bit  rate, 

R*  -  .07  R^  (3. 4. 3. 1-5) 

See  Figs.  19  and  20  for  plots  of  the  above, 
resulting  in 


/,TGTG,J;2*sin(*)sec2(*)exp(-“^)\s 


9'l(.Pblna.)D 


where  all  the  parameters  are  defined  above. 


(3.4.3. 1  -6) 


The  above  expression  for  the  maximum  throughput  Tu*  allows  us  to  quantify  the 
affects  of  various  system  parameters  and  constraints  on  the  best  throughput  using 
a  constant  bit  rate.  As  an  example  we  plotted  the  variation  of  T„*  as  a  function 
of  wavelength  for  our  sample  system  (Sec  3.3).  See  Fig.  17  bottom  curve.  Using 

our  sample  system  with  a  wavelength  of  6  m.  (50  MHz)  we  get  R _ -  86.5  kb/s, 

R*-6.055  kb/s  and  Tu*-  214  b/s .  Additional  plots  are  given  by  Figs  21  and  22. 


3.4. 3.2  Opt fill  A gfiEflgft  fllEgfflrtiPUS  isx.  Constant 
Bit  Bata  SiattiB  zz.  azardsnaau. 

In  this  section  the  optimal  average  system  throughput  for  the  underdense 
burst  is  given.  This  is  the  maximum  average  number  of  bits  per  unit  time  for 
the  'best'  transmission  rate  R  (and  other  system  constants)  assuming  only 
overdense  bursts  are  present. 

This  optimal  throughput  is  found  by  first  deriving  the  dependence  of  the  average 
throughput  (Eq.  3. 4. 3-1)  on  the  bit  rate  R,  and  solving  for  the  optimal  bit  rate 
R*  that  would  yield  the  desired  optimal  throughput  T0*.  The  solution  is  done 
with  respect  to  the  normalized  bit  rate,  R^j,  which  yields  a  result  independent 
of  a  particular  choice  of  system  values. 
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For  the  given  system  Pw...  N0,  R,  a  and  C0  are  specified.  With  a  reference 
to  Eq.  (3. 4. 1.2-3)  and  the  accompanying  disscussion  we  rewrite  --  the  bit 
rate  normalized  with  respect  to  R,.^.  the  maximum  bit  rate  achievable  for  the 
given  system  using  overdense  bursts  --  and  related  parameters: 


i?™,,.-  max{qT.qL0} 


where,  q^  -  1016  and  -  101* 

a  and  C0  are  defined  in  Eq.  (2. 1.5-2). 
N0  is  defined  in  Sec.  2.4 


(3. 4. 3. 2-  1  ) 


We  note  from  the  above  that  for  qT  <  q^  (or  equivalently  <  .1)  where 
-  quj  and  the  other  for  qt  >  q^  (or  equivalently  RN0  >  .1)  where  q,^,,,,  -  qT.  In 
addition,  the  value  of  Q  (Eq.  (3. 4. 1.2-3))  varies  with  the  R^j  as  well  since  Q 
is  a  function  of  -  see  Eq.  (2. 1.6-2).  For  the  case  where  R„0  >  .1  (qT  > 

qnj)  we  have 

Q  “  -5(quo)  5  R«o  /  (l'R«o) 

For  the  case  where  R^o  <  •  1  (qj  >  quj)  we  have 

Q  -  •5(qLO)-5  /-9 

Using  the  above  relationships  in  Eq.  (3.4. 1.2-3)  to  express  the  average  overdense 
burst  duration  aa  a  function  of  R^j  we  get: 


For  Rho  5  0.1  : 
'iQioquoa-  Rn 


[  0.01  , 

~~p2  1  +  cos  (10tf*0)-cos  '(**„) 

K  NO 


For  Rjrj  >0.1  : 


(3. 4. 3. 2-2) 
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The  average  number  of  bits  per  overdense  burst  is  by  definition 

N.0-R~.  (3. 4. 3. 2-3) 

and  can  be  expressed  in  terms  of  the  normalized  bit  rate  using  Eq.  (3.4. 3.2-1) 
for  R  as  a  function  of  and  Eq.  (3.4. 3.2-2)  for  the  average  overdense  burst 
duration. 

Since  the  throughput  T0  according  to  Eq.  (3.4. 3-1)  is  the  product  of  the  arrival 
rate  and  the  last  equation  i.e. 


To-  Ao  '  N go 


(3. 4. 3. 2-4) 


we  need  to  express  the  arrival  rate  of  overdense  bursts  Aq  as  a  function  of  RN0 
We  start  by  employing  Eq.  (2.3-2)  with  q!  set  to  q^,,,,  (the  lowest  overdense 
q  as  determined  from  the  choice  of  i.e.  bit  rate)  and  q2  set  to  quo  (the 

physical  maximum).  The  dependence  of  q^,,,,  on  R^q  is  then  used  (disscussed  above 
after  Eq.  (3.4. 3. 2-1))  yielding  the  following  arrival  rate: 


(  .99- 

1 

A- 

R  NO  -  \ 

A0~l 

1 

) 

V  ^ 

^.01  •  A ' 

*N0>-  1  / 

(3. 4. 3. 2-5) 


where , 

v  3»sin  »(**»* 

Qco 

quo  -  1016  and  q^  -  101* 

By  direct  substitutions  from  the  above  equations  we  express  the  throughput 
using  overdense  bursts  as  a  function  of  normalized  bit  rate  and  system  parameters 
as : 


.33-  10 


i#  />rCTCjt^2»sin^sec3»  -  X3 


KU~l~yl°Kio~ 


i-*, 


N0  g-'(P>ma,)- D  Rr 

1  cos' '(10#*o)- cos' '(#*,<,) 

P  NO  >  ■  1 


1  -cos''(tf*o) 

K  NQ 


R  NO  -  ■  1 


(3. 4. 3. 2-6) 
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where , 


Rt  and  4> 


RT 

gt>gr 

X 

D 

No 


idealized  antenna  beamwidth. 
geometric  parameters  detrmined 
by  link  distance  (from  section  2.2) 
transmitter  power. 

transmitter  and  receiver  antenna  gains, 
carrier  wavelength. 

diffusion  coefficient  of  the  atmosphere, 
power  spectral  density  of  the  noise. 
(Sec.  2.4) 

minimum  of  bit  energy  to  N0  ratio. 


Note  that  the  equation  is  written  in  two  parts;  the  first  incorporates  all  the 
parameters  of  the  particular  system  and  the  second  part  reflects  the  variation 
in  throughput  with  choice  of  bit  rate.  The  above  throughput  (  multiplied  by 
2.55)  for  our  sample  system  (Sec  3.3)  is  plotted  in  Fig.  18. 

We  would  like  to  find  from  the  above  equation  the  optimal  normalized  bit 
rate  R^*  which  yields  the  maximum  throughput  T0*  .  This  optimal  normalized  bit 
rate  for  any  given  system  is  found  numerically  to  be  (see  Fig.  18): 


R»o*  -  1 

or  equivalently  the  optimal  bit  rate, 

R*  -  -1  Rmaxo  (3. 4. 3. 1-7) 

The  result  is  intuitively  obvious  since  R^o-.l  is  the  maximum  bit  rate  that  still 
includes  all  overdense  trails  i.e  maximum  arrival  rate.  For  a  higher  bit  rate 
we  have  a  smaller  arrival  rate  that  reduces  the  throughput.  A  smaller  bit  rate 
will  not  change  the  overdense  arrival  rate  but  will  diminish  the  throughput. 

Using  the  optimal  bit  rate  for  overdense  bursts  we  get  the  maximum  overdense 
throughput  T0*  by  substituting  Eq. (3 .4 . 3 . 2-7)  into  Eq.  (3.4. 3. 2-6)  yielding: 


To - .3078-  10 


/>rGrG*tJ$sin$sec,$-  Xs 
A fog'(.P^o,)DRj 


(3. 4. 3. 2-8) 


where  all  the  parameters  are  defined  above. 
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The  above  expression  for  the  maximum  throughput  T0*  allows  us  to  quantify  the 
affects  of  various  system  parameters  and  constraints  (wavelength,  power,  link 
distance,  maximum  bit  error  rate  etc.)  on  the  best  throughput  using  a  constant 
bit  rate.  Using  our  sample  system  (Sec  3.3)  with  a  wavelength  of  6  m.  (50  MHz) 
we  get  R*  -  252  kb/s  and  T0*  -  8688  b/s. 


3. 4. 3. 3  Total  6SSXAZS.  Throughput  tat 

Constant  Bit  Rate  System. 

In  this  section  we  combine  the  result  of  the  last  two  sections  to  present 
the  total  throughpjt  and  optimal  throughput.  Since  we  used  differently  normalized 
variables  for  underdense  and  ox  :rdense  we  shall  now  employ  a  single  normalized 
variable;  the  one  corresponding  to  the  overdense  case. 


Recalling  previous  definitions: 
For  underdense  we  used 

2 


^  <?</ 


where 


R 


max 


C  uQu 

N0g'(P>na.) 


(3. 4. 3. 3-  1  ) 


the  maximum  bit  rate  achievable  in  a  given  system  using  underdense  bursts.  Cy 
and  N0  are  defined  in  Eq.  (2.11-2)  and  Eq.  (2.4-2),  respectively. 


For  overdense  we  used 


A  maxo 

where 

p  *  Cq^ a  q i/o 

ie-N0- Sf'‘(F*„„) 

where,  q^  -  1016  and  q^  -  101* 

a  and  C0  are  defined  in  Eq.  (2. 1.5-2). 
N0  is  defined  in  Sec.  2.4 


(3.4.3  3-2) 


From  the  above  we  can  write 
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-  s  -  6.446  •  •  sec* 

&  m«* 

and 

*N~  *•  *V0 


(3. 4. 3. 3-3) 

(3. 4. 3. 3-  4) 


Note  that  Rjj— 1 ,  the  maximum  underdense  bit  rate,  corresponds  to  R^q  -  1/s. 

Using  Eq.  (3. 4. 3. 1-4)  for  the  underdense  throughput,  with  R„  expressed  by  the 
last  equation,  and  Eq .  (3. 4. 3. 2-6)  for  the  overdense  throughput  we  get  for  the 
total  throughput  T : 


~r  ^  ,„-i»  frCrG,l;2*sm*sec2*- \s  _  ^  _ 

/  •  .509'  10  •  - ; - - - {T  ,.  *  T  0,  *  T  n?> 

No- 9~'{P  D-  R7  L  02t 


(3.4.3.3-S) 


where 


(  33.359  \^sRh0-sRn0  .  f  \ 

T  v  **  exp^~  j  ~i~  (s  Rtl0)p  *  *  P(sRvgy  ■  \n(sR  V0)J  ■  l/(^- - R, 


A'sec 


secpRlg  /  1 

0,‘  1401  V  Ri0~ 


1  -cos  '(10/fvo)-cos ~\Rmo)  •  6'(.1  -  j?no) 


T02  -  .649  sec  t  • 


■wo  L  \  ano 


—  -  1  -cos  '(#„0)  •  L'(R„0- A  ) 


s-6.446  e(”3”Mi«c!»,sec* 


/I  for  x>0\ 

"“’■(o  lor  x  <0j 

is  the  unit  step  function. 

All  other  parameters  are  defined  in  previous  sections. 


3-5  Variable  Bit  Rate  and  Constant  BER  System 
3-5.1  Channel  Duration 

For  a  variable  bit  rate  system  the  bit  rate  is  changed  so  as  to  maintain 
a  reasonably  constant  energy  per  bit.  Clearly  it  implies  a  constant  BER,  Pb, 
for  all  transmitted  bits.  Considering  for  example  the  underdense  burst,  if 
apriori  knowledge  of  the  burst  were  available,  we  could  divide  the  received  power 
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curve  into  equal  energy  portions  and  compute  the  resulting  individual  bit  durations 
before  ve  transmit.  Under  such  conditions  all  bits  will  possess  exactly  the 
same  energy  and  hence  the  same  BER.  The  total  number  of  bits  transmitted  for 
that  burst  will  simply  be  proportional  to  the  total  area  under  the  received  power 
-versus-  time  curve.  Since  we  don't  have  an  apriori  knowledge  of  the  burst  but 
rather  an  instantaneous  (within  a  negligible  delay)  indication  of  channel  strength 
we  can  only  evaluate  the  necessary  duration  for  the  current  bit  and  transmit  it 
as  such.  During  the  duration  of  the  transmitted  bit  the  power  received  drops. 
If  the  bit  duration  is  small  the  resulting  drop  in  bit  energy  is  negligible  and 
our  assumption  of  constant  bit  energy  is  not  violated.  As  time  progresses, 
however,  the  bit  duration  increases  (power  received  decays)  and  the  drop  in  bit 
energy  relative  to  the  assumed  constant  bit  energy  becomes  appreciable  and  finally 
intolerable.  This  fact  limits  the  maximum  bit  duration  which  in  turn  restricts 
the  duration  of  transmission  or,  equivalently,  channel  duration. 

3. 5. 1.1  Transmission  Duration  for 

Variable  £it  Rets  Underdense 

We  start  by  imposing  our  fidelity  constant,  k,  for  the  normalized  bit 
energy  --  the  ratio  of  the  real  bit  energy  to  the  desired  one.  For  our  derivation 
we  will  assume  k  >  .9  i.e.  no  more  than  10%  drop  in  bit  energy  shall  be  tolerated. 
We  further  define: 

Pb  Assumed  constant  BER  for  all  transmitted 

bits . 

Tbi  Duration  of  the  1th  bit 

R(t)  Instantaneous  bit  rate 

T0  Duration  of  the  first  bit,  1/Rg 

N  Last  bit 

T(m)  Total  duration  of  the  first  m  bits; 

also  starting  time  of  the  (m+l)^  bit: 

m 

i- 1 

Tt  Duration  of  transmission  (burst  duration)  - 

T(N) 

c  Max  (TbH/B)  for  a  given  k. 


In  order  to  maintain  a  constant  BER,  Pb,  we  have  for  the  variable  bit  rate 
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*(*)■ 


Cuq 1 


N0  gr'(/\) 


-l/J  -  D  . 

e  m  K  o  fi 


and  the  n01  bit  duration  is  given  by 

T(»-l) 

Tt.-T 0e  * 


where 

T0  -  1/Rq  -  N0  g-i  (Pb)/Cu  q2 

For  k  >  . 9  we  have  for  the  last  bit,  N: 

7(lA,<tl.5-V2.25-6(l  -*)]•  Bm  c  -  B 


Combining  with  eq  (3. 5. 1.1-2)  we  have 


T  _  T  o 
B  "  B  ' 


<  c 


Since  the  transmission  duration  is  given  by 

Tt  -  T(N-l)  +  TbN  we  have  from  the  last  3  eqns . 


(3.5.1 .1  -  1) 


(3.5.1 .1  -2) 


(3.5.1  .1-3) 


(3.5.1  .1-4) 


7,<5[c*ln 


(3.5.1  .1  -5) 


for  k  -  .9,  c-  .2  we  have  from  the  above  Tt/B  -  .2  +  In  (RoB/5). 

For  all  practical  systems  we  can  safely  use 

(3. 5. 1.1-6) 

where  Ro  is  defined  in  Eq.  (3. 5. 1.1-1) 


3.5.2  Throughput  fat  Variable  Bit  Rate  System 
c.  Pnderdense 

Using  the  results  from  the  previous  section  on  transmission  duration  we 
can  find  the  approximate  number  of  transmitted  bits  by  integrating  R(t)  from 
time  0  to  time  Tt: 


Nau-RoB-CuBqVg-KP^No  (3. 5. 2-1) 

where  R0  is  the  variable  bit  rate  at  time  0. 

and  since  B  and  q  are  independent  the  average  of  NOT  over  the  joint  space  of  q 
and  B  yields: 
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N  uv 


CuBq 2 


g'l(P>)N  o 


(3. 5. 2-2) 


The  derivation  of  the  above  equations  is  in  section  4.5 

We  shall  assume  that  the  bit  rate  shall  vary  between  R^nv  and  R,^,.  Bearing 
in  mind  that  for  a  given  bit  rate  R  we  have  a  corresponding  q  such  that 
Cuq2-N0  g'\P,)  R  qL<q<qu  (3.5. 2-3) 


As  a  result  we  have 

Q 


N0  •»•'(/»*) 


Cu 


■JT. 


Q  maxV 


H0g-'(Pb) 


J*. 


(3. 5. 2-4) 


(3. 5. 2-5) 


R^nv  must  be  chosen  such  that  q^y  >  qL,  the  physical  minimum. 

R^v  must  be  chosen  such  that  qmfT,g  <  qUt  the  physical  maximum. 

Note  that  because  our  empirical  arrival  rate  relationship  (Sec.  2.3)  is  valid 
for  q>1013  we  shall  set  R^y  such  that  q*^  >  1013. 

To  get  the  most  throughput  out  of  our  system  we  set  R^.Tg  to  its  maximum  value 
i.e.  corresponding  to  q-qu.  We  thus  have  set  q^y-qu  and  using  it  in  the  last 
two  equations  we  have 


Quay 
9  no  xy 


where 

^  naxV  “ 


RfimV  is 

c  uqj 

No  g~'(P>) 


given  by  Eq.(3.5.2-5) 


as 


(3. 5. 2-6) 


(3.5  2-7) 


where,  C„  is  from  Eq.  (2.1. 1-2) ,  N0  is  from  (2.4-2)  resulting  in  Eq .  (3 . 5 . 2- 13) . 

Note  that  RTT.q  -  R _  the  maximum  underdense  for  constant  bit  rate  system  if  we 

set  P>.-Pi _ _  a  direct  result  from  Eq.  (3.4. 3. 1-1) 
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where,  (Eq.  (3. 5. 2-7)  ,  Eq.  (2. 1.1-2)  and  (2.4-2)) 


RmaxV  -  2.518  •  10,# 


„  _  ^  ,  (  33.339  \ 

PTCTC,\  exp(-K,m^j 
(.063X*3-  I) g'V*)*? 


t, 

Rt  and  <(> 


P  T 

Gi>G* 

X 

D 

No 


g'HPb^x) 


Rainv/R««v  ;  Rmwv  defined  in  Eq.  (3. 5. 2-7) 
R^v  is  arbitrary  yet  if  resulting  R^,  is 
small  system  complexity  and  cost  increases, 
idealized  antenna  beamvidth. 
geometric  parameters  detrmined 
by  link  distance  (from  section  2.2) 
transmitter  power. 

transmitter  and  receiver  antenna  gains, 
carrier  wavelength. 

diffusion  coefficient  of  the  atmosphere, 
power  spectral  density  of  the  noise. 

(Sec.  2.4) 

minimum  of  bit  energy  to  N0  ratio. 


(3.5.2-  13) 


For  our  sample  system  (Sec.  3.3)  we  chose  R^y-lOOO  b/s  i.e  q^nv  -  1.1  1013- 
g*1(Pb)-9  we  get  Rj^-. 011558  and  Tuy  about  1150  b/s. 

If  we  assume  the  constant  bit  energy  of  the  variable  rate  system  tobe  the  same 
as  the  minimum  bit  energy  of  the  constant  rate  system  then  for  our  sample  system 
(  underdense  only)  we  improved  the  throughput  from  214  (Sec.  3. 4. 3.1)  to  1150 
-  a  ratio  of  5.3. 

3.5.3  Throughput  for  Variable  Bit  Rate  System 
z.  Overdense 

The  maximum  power  point  of  an  overdense  burst  Pmax  is  well  aproximated  by 
(Eq.  (2. 1.5-2)  neglecting  k) 

(3.5.3- 1) 

We  recall  that  for  our  variable  rate  system  the  rate  varies  as  the  power  scaled 
by  the  bit  energy  N0g_I(Pb) 

This  variation  in  time  can  be  approximated  as  a  triangle  whose  base  lies  on  time 
axis  from  the  origin  to  the  point  given  by  aq  and  whose  top  vertex  is  at 
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Pmax/Nog'^Pb)  • 

Since  the  number  of  bits  sustained  by  this  overdense  burst,  Nbv, 
equal  to  the  area  underneath  the  rate  time  function  we  have 


N 


$v 


1  P  m«»  _  _  1  C0a3'2  3/2 

2A/0  g-\P,)a'  ~2TeNog''(.P»)q 


is  theoretically 


(3.S.3-2) 


Averaging  over  the  ensamble  of  overdense  bursts  using  Eq.  (2. 1.6-2)  with  q^. 
no-<ko“10u  (qio/quo  “  -01)  we  get 


which  yields  for  average  number  of  bits 


1  C0a3/2 

4  feN0g-'{Pb) 


and  the  arrival  rate  is  given  by  Eq.  (2.3-2) 

A  3- <1lo\ 

Qlo  \  q  i/o ) 


(3.5. 3-3) 


(3.5.3-  4) 


(3. 5. 3-5) 


Combining  all  of  the  above  for  the  throughput  expression  (varible  rate  and 
overdense  bursts)  we  get 


1  94.686  •  10' 


No\g'(P>)- D- Rt 

For  our  sample  system  (Sec.  3.3)  we  get  a  throughput  of  26.7  kb/s. 


(3. 5. 3-7) 


If  we  assume  the  constant  bit  energy  of  the  variable  rate  system  to  be  the  same 
as  the  minimum  bit  energy  of  the  constant  rate  system  then  for  our  sample  system 
(Sec  3.3)  (assuming  overdense  only)  we  improved  the  throughput  from  8688  b/s 
(Sec.  3.4. 3.2)  to  26.7  kb/s  -  an  improvement  ratio  of  3. 
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3.5.4  Comparison  fif  Throughout  fsr  Const SHE  find 
Variable  Ii£  Rate  Systems 


In  this  section  we  shall  apply  the  derived  relationships  in  the  previous 
sections  to  a  practical  system.  The  various  parameters  of  interest  are  computed. 
The  throughput  of  the  system  for  the  optimal  constant  bit  rate  is  computed  and 
compared  with  that  of  throughput  of  the  system  using  variable  bit  rate.  This 
comparison  will  demonstrate  the  improvement  in  throughput  for  variable  bit  rate. 
We  start  by  copying  the  parameters  of  the  system  from  Sec.  3.3: 

L  -  1000  km  link  distance.  From  section  2.2 

we  have  for  the  next  two  parameters: 
sec($)  -  18. S  and 


RT~  R  ,  -514  km. 

PT  -  1000  w 

Gt  -  10  dB. 

Gjj  -  10  dB.  And  for  antenna  beamwidth  we  use 

£  -45.84°  (.8  rad) 


X.  -6  m.  (50  MHz) 

D  -  8  m2/s 

Pixmx  -  10"5  Max.  BER  for  constant  bit 

rate  system. 

Using  BPSK  modulation  we  get: 


g-l<Pta«)  “  9 

N0-4-  10'2l[80^]  +2.5 


W/Hz 


See  discussion  above  Eq.  (2.4-2). 

N0  -  4.89  10'20  at  50  MHz. 

For  comparison  the  BER  for  variable  rate  system,  Pb> 
is  set  equ;1  +0  P^j* 

R^y  -  1000  b/s  (minimum  variable  rate  allowed) 


Underdense  Comparison 


The  improvement  we  get  for  underdense  bursts  is  arrived  at  under  the  assumption 
that  the  optimal  rate  is  used  for  the  constant  bit  rate  system  and  for  the 
variable  bit  rate  system  some  R^y  is  assumed.  In  addition,  we  assume  that  the 
Pb-Pbmax  i.e.  the  constant  bit  energy  (variable  rate)  is  assumed  equal  to  the 
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minimum  allowed  bit  energy  (constant  rate  system) .  We  thus  divide  the  relevant 
equations:  Eq.  (3.5.2-12)  and  Eq.  (3.4. 3. 1-6) 
resulting  in  the  improvememt  I0 


where , 


R 


NV 


R  mtmV 

R  maxV 


RmaxV  -  2.51 8  •  10 


P ^ C^C|X^fixp 


(- 


33  339  'N 

X2I«C2J 


( .063X.2  3  +  Rj 


(3.5.2-  13) 


The  above  is  evaluated  using  the  values  of  our  sample  system  where 
R«v-R»inv/R»«  -  1000/86500  -  .01156 
resulting  in 

lu  -  5 . 38  or  about  7  dB 

ov$rden?e  cqippsrisQn 


The  relevant  equations  to  be  used  for  comparison  of  throughputs  are  Eq.(3.5.3-7) 
and  Eq.  (3.4. 3. 2-8)  resulting  in  an  improvement  of 

Io  “  Tot/T*o  “  3 


BAaigaalsn 

We  shall  compare  now  the  improvement  in  throughput  of  adaptive  rate  relative  to 
fixed  rate  systems  for  underdense  in  terms  of  its  components;  the  arrival  rate 
and  the  average  number  of  bits  transmitted  per  burst.  For  fixed  rate  (underdense 
bursts)  we  have  the  following  results: 

Fixed  Rate: 

Rbmz  “  86-5  kb/s 
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Since  the  normalized  optimal  bit  rate  is-. 07  we  have 
R*-  .07  Kmx  -  6.05  kb/s 
For  our  system  we  also  have 

B - .264 sec 
Cu  -  3.8  •  10'42 

From  eq  2.4-2  we  have: 

N0  -  4.89  IO-20 

For  this  rate  we  can  compute  some  parameters  for  general  interest: 

the  average  number  of  bits  per  underdense  burst  is  (from  Eq.  (3. 4. 3. 1-2))  1 . 846 
kb/s  per  burst 

q^n  -  .26457  10u  e/m  (  from  eq.  3.4. 1.1-1  set  R  to  R*> 

Which  yields  an  arrival  rate  of 

.116  bursts  per  sec  or  equivalently  an  average  inter-arrival  time-8.6  s. 

The  product  of  the  last  two  results  yield  a  throughput  of  214  b/s 

Variable  Rate 

In  general  the  ratio  of  the  bursts'  arrival  rates  of  variable  bit  rate  to  fixed 
bit  rate  system  is  given  by 

Av/A*  -  .36  {  (R^jt/Kainv)  5  -1  )  (Eq.  (3.4. 3.1-3)  and  Eq .  ( 3 . 5 . 2- 11) 
which  in  our  example  yields  2.9. 

The  average  number  of  bits  per  burst  using  variable  rate  is  from  Eq.  (3.5.2-10) 
3.328  kb/s  per  burst . 

The  ratio  of  improvement  for  the  average  number  of  bits  per  burst  is 
3.328/1.846  -  1.8 

whereas  the  improvement  in  arrival  rate  is 
2.9  yielding  a  total  improvement  of  5.3 

It  is  interesting  to  note  that  the  larger  improvement  is  due  to  the  arrival  rate. 
The  arbitrary  set  RminV,  the  minimum  rate  allowed  in  adaptive  rate  system,  sets 
a  lower  threshold  of  electron  line  density,  q^y,  whic  utilizes  bursts  that  are 
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precluded  for  the  fixed  rate  system  with  its  higher  electron  line  density  threshold 
<Wn-  It  should  be  noted  that  however  that  system  hardware  complexity  increases 
with  lower  R^y. 
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4  DERIVATIONS 

This  chapter  supplements  the  derivations  outlined  in  the  previous  chapter 
for  the  results  that  have  thus  far  been  presented. 


4.1  Link  Geometry 


With  reference  to  fig.  6  and  the  definitions  in  section  2.2  we  have  the 
following  relationships: 


V-L/4R, 

(4.1-1) 

n 

Y--V 

(4.1-2) 

x  -  2i,sinJ(v) 

(4.1-3) 

y-xsin(v)-2£,sin2(ip) 

Since  it's  always  true  (L  <  2000 

km)  that  (L/4R,)2 

«1  we  have 

(4.1-4) 

l  2 

sin2(tp)-  n>2-  — 

(4.1-5) 

By  the  Law  of  Cosine  we  have: 
ft2  +  xz-2xh  cos(y)  “ 

(4.1-6) 

We  substitute  for  x  and  y  and 

sin2(v)  from  eqns. 

(4.1-3),  (4.1-2), 

(4.1-5) 

respectively  in  eqn.  (4.1-6)  and  arrive  at  the  stated  result  in  section  2.2  after 
some  trivial  algebra. 

V(h+£z/8ff.)2+IJ/4  (2.2-1) 


From  the  figure  we  also  have: 

2  *5 

sec2(4)-— - — 

(h  +  y)2 

which  when  combined  with  eqs.(4.1-4),  (4.1-5)  and  (2.2-1) 
yields: 


3ec2(4)  -  1  ♦ 


(2.2-2) 


and  thus  completing  the  derivation. 
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4.2  Derivation  ffif  Average  Burst  BamiSIL 
Constant  Bit  Rate 


This  derivation  is  for  the  underdense  case.  We  start  with  equation 
(3. 4. 1.1-1)  as  the  expression  for  t8  Since  we  want  to  average  over  the  joint 
space  of  q  and  B  which  are  statistically  independent.  The  joint  pdf  is  given 
by  the  product  of  the  individual  pdfs.  Because  B  acts  as  a  coefficient  in  the 
expression  for  tB  the  two  integrals  are  separable.  (The  individual  pdfs  for  q 
and  B  are  from  equations  (2. 1.2-2)  and  (2.1. 3-2,3)  respectively. 

<5-2  m~P-Q  f‘(ln  (q/qm,jyqpdq  (4.2-1) 


The  above  integral  can  be  solved  by  employing  the  following  substitutions: 

z  -  (p-1)  In  (q/qnan) 

q  -  q«in  e*/<P-*> 

dq  -  [qain/Cp-l)]  «*/<*rl>  dz 

lnn(q/qaiin)-[z/(p-l)  ]n 

Reduces  the  above  to  a  single  integral  of  the  form 


2 «g«  i  ,*"«>">  , 

t ,  - - -  z  e  dz 

(p-  1)‘  1  ~r  J0 


(4.2-2) 


where 


~  N0  - R  g  '{Pbnaxy 

l  q</  J 

CuQu 

which  can  be  solved  recursively  by  parts,  yielding 


t»m 


2*B* 

(p-1)* 


n  • 

1  -x 


1  -x 


ln‘(x)\ 


Setting  n-1  and  n-2  in  the  above  yields  the  necessary  equations  for  the  expected 
value  and  variance  of  tB  although  some  wading  through  algebra  is  necessary  for 
the  derivation  of  the  variance. 


4.3  j&£3£  BaERtian  Constant  Bit  Rate  ji  Qverdense 

Since  the  roots  of  the  overdense  power  received  equation  cannot  be  solved 
for  directly  we  use  a  second  order  approximation.  In  fact  each  side  of  the 
ov'  rdense  'mountain'  is  fit  with  a  parabola.  The  parabolas  are  expanded  about 
the  maximum  power  point.  For  the  left  parabola  (which  approximate  the  left  side) 
we  use  the  smaller  root  and  subtract  it  from  the  larger  root  of  the  right  parabola 


(which  approximates  the  right  side  of  the  power  equation).  This  difference  is 
the  burst  duration  for  the  overdense  trail.  We  start  by  copying  eq.  (2. 1.5-2) 


/>(<)- C0<./(*+t)  In 


(S?) 


tiaq-k 


(2.1 .5-2) 


where , 


ptctg,\2 

C0-3. 16628610  3Jd— —  /- . 

P  T 


k -  . 105625/D 


a  -  7.14314310 


,lT  X.2sec2$ 


We  need  to  solve  for  the  roots  t2,  tx  of  the  above  eq.  when  P  is  set  to  Pmln  since 
their  difference  is  our  transmission  interval  -  burst  duration.  To  do  so  we 
first  modify  the  above  eq.  as: 

PN-yJr- ln^J  k/aq<r<  1  (4.3-2) 

or  equivalently 

rln(r)-P2-0  (4.3-4) 

where 

PH  -  P  /  Co  (aq)W2 
r-(k+t)/aq 

The  above  will  yield  the  roots  r2,  rx  for  P-P„,ln 
Note  that  the  burst  duration  is  defined  as 

(4.3-4) 


We  are  now  ready  to  approximate  r2-rx  where  rx  is  the  'left'  root  and  rx  the 
'right*  root.  Note  that  the  above  eq.  has  a  maximum  PH  of  17(e)1'2  at  r-l/e. 


To  approximate  rx  we  expand  ln(r)  about  r-l/e  and  substitute  the  first  two  terms 
for  ln(r)  in  the  last  eq  resulting  in:  r2  -  (2/e)r  +  (l/e)P„2  -  0 
Solving  and  using  the  smaller  root  we  thus  approximate  rx  as: 


r  ,  - -- (  1  -  V  1  -  3/»S) 


(4.3-5) 


To  approximate  r2  we  expand  rln(r)  about  r-l/e  and  substitute  the  first  three 
terms  in  eq  (4.3-3)  and  solve  for  the  larger  root. 
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ra---(l*V2(l -•/»;)) 
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r,-r,  Vl  -ePl-0.88jl-ePl 


(4.3-6) 

(4.3-7) 


A  comparison  of  the  above  with  an  iterative  solution  of  the  real  (without 
approximation)  roots  of  the  eq  (4.3-3)  and  their  difference  allows  us  to  replace 
the  .88  constant  in  the  last  eq.  with  .95 

r2-r!  -  .95  root  [1-e  P  sub  N  super  2] 

Since  the  power  for  which  we  want  to  find  the  roots  is  Pmtn  from  eqs .  (4.3-3)  and 
(4.3-4)  we  have  the  desired  approximation  for  tB 


.95a 


q- 


\_ 

q 


(4.3-8) 


which  is  the  result  stated  in  (3.4. 1.2-1). 

Note  that  q>qaiino  “  max  l^T.  qLo)  •  If  (which  is  set  by  Pmin  and  system  parameters 
reflected  in  Co  and  a  is  less  than  qLo  than  the  minimum  overdense  electron  line 
density,  qmino.  is  set  to  the  physical  minimum  overdense  electron  line  density 
q^.  On  the  other  hand  if  qT  is  bigger  than  qLo  than  the  overdense  trails  with 
electron  line  density  between  qLo  and  qT  have  zero  burst  duration  and  are  not 
"seen"  by  our  system. 


For  most  practical  systems  which  utilize  underdense  bursts  Pmin  is  such  that: 
qT  «  qLo  I  • e  •  q^ino  _<?Lo  and  CB  is  approximated  as  aq,  (q>qLo  )  and 
E( tB)  -,95a  E(q)  -  a(qu0  QLo')V2  ’  a  10». 


Derivation  of  CDF  and  PDF: 

From  eq .  (3. 4. 1.2-1)  we  note  that  tB  is  monotonically  increasing  with  q 
such  that  for  any  tB<tj  we  have  correspondingly  a  q<q:.  For  tB  -  tj  we  have  from 
eq.  (3. 4. 1.2-1)  : 


where 


x-aqT/2-(e/2)(Pmln/C0)z 


From  the  above  discussion  there  is  a  minimum  tlt  tlmln  corresponding  to  qffilt)0  and 
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a  maximum  tx,  corresponding  to  quo,  the  max  overdense  electron  line  density. 

We  can  now  relate  the  following  probability  events: 

Prob{(,</,)-Prob{q„llw^q<q1>  t  Imln  <  t ,  <  t ,  „„  (4.3.9) 


We  get  the  CDF  of  q  using  eq  (2. 1.2-2)  such  that  the  right  hand  side  of  the  last 
eq.  becomes 

Q*  U*<W>  /  <h)^) 

Eq  (4.3-9)  then  becomes 


■{■-w 


(4.3-10) 


Substituting  eq  (4.3-8)  into  the  above  eq.  we  get  the  desired  eq .  (3.4. 1.2  -6). 
To  get  eq.  (3.4. 1.2 -7)  we  simply  take  the  derivative  of  eq.  (3.4. 1.2-7) 


4.4  Optimal  Average  throughput  for 

constant  bit  rate  ^  Onderdense 

Our  first  objective  is  to  derive  eq.  (3. 4. 3. 1-2)  from  section  3.4.2  we 
know  the  number  of  bits  per  given  underdense  burst  for  a  given  R  and  a  set  of 
specified  system  parameters  is  given  by 

Ngu  ”  RtB  (4.4-1) 

Averaging  NBU  over  the  ensemble  of  underdense  bursts  (q  and  B  independent)  is 
given  by 

NtumRt «  (4.4-2) 


rewriting  the  result  in  Eq.  (3.4.1. 1.-3)  which  is  derived  in  Sec.  4.2 

-  2  B  1 

<*_FrT'iTT'[1"X(1“lnx)]  (3. 4. 1.1-3) 

where 

C  vlv 


and 
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we  substitue  for  x  and  R  in  terms  of  their  R„  definitions  and  write 


(Rn)p’  -  x  p'  -  (p- 1)/2  (p  -1.6) 

R  “  RN 

Substituting  these  into  eq .  (4.4-2)  yields  the  desired  eq. 

_  n  ~5  i 

in«)j 

P  1  *  K  N 


(3. 4. 3. 1-2) 


(3.4.3. 1  -2) 


This  eq.  specifies  for  a  given  choice  of  R  the  average  number  of  bits/burst. 
The  rest  is  in  Sec.  3.4. 3.1 


4.5  Transmission  Duration  &  Throughput 

for  Ml  Eat &  System 

With  reference  to  section  3. 5. 1,1  we  define  again  the  following  terms: 

Pb  Assumed  constant  BER  for  all  transmitted 

bits . 

Tbi  Duration  of  the  ith  bit 

R(t)  Instantaneous  bit  rate 

T0  Duration  of  the  first  bit,  1/R^ 

N  Last  bit 

T(m)  Total  duration  of  the  first  m  bits; 

also  starting  time  of  the  (m+l)th  bit: 

m 

i  •  1 

Tt  Duration  of  transmission  (burst  duration)  - 

T(N) 

c  Max  (TbH/B)  for  a  given  k. 
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We  wish  to  maintain  the  actual  bit  energy  above  a  certain  fraction  say  90% 
of  the  assumed  bit  energy.  This  will  insure  a  relativly  contant  BER  corresponding 
to  the  assumed  bit  energy.  Since  the  last  bit,  the  Nth  bit,  has  the  longest 
duration  (power  is  monotonically  decreasing)  the  loss  in  bit  energy  relative  to 
the  assumed  level  is  the  largest,  n  equation  form: 

Pb  actual  ^  ^  ^b  aaauaad 

i.e. 


rT  ** 

I  P(UN- 
J  0 


l))e'"4dt>Jfc-  P(7-(JV-  !))•  r„„ 


which  implies 

1  -  kT'bH  >  exp(-T'bl))  where  T'b„  -  TbH/B 
and  assuming  T'b„  <  1  we  use  the  first  four  terms  in  the  series  expansion  of 
the  exponential  (about  0)  and  get  the  root  of  T'bN  yielding 


The  quadratic  equation  in  T'b„  yields 
T\N-T-j-<\.S-42.2S-6(\-k).c 


From  the  above  definitions  of  the  bit  duration  we  have 

Tb,/B  -  [l/Ro  B]  exp  (T(N-1)/B)  <  c 

yet  by  definition  T(N-l)  -  Tt  -  Tb„ 

we  have  then  Tt  <  B  [c  +  In  (cRqB)] 

which  is  the  result  stated  ineq.  (3. 5. 1.1-5) 

We  now  evaluate  the  number  of  bits  for  a  given  burst  using  variable  bit  rate  as 

J  *(<)d*»*oij[l-exp^—  JJ 
0 

Using  the  maximum  Tt  from  eq.  (3.5. 1.1-4)  i.e.  letting  NB  -  maximum  NB  we  have 
NBu  “  Ro®  +  1  *l/c 
For  K- . 9  c- . 2 

For  all  practical  system  R0B»4 
such  that 

Nbu  -  RqB  -  cu  B  qVq'^Pb)  N0 

as  stated  in  section  3.5.2  Eq.  (3. 5. 2-1) 

The  average  throughput  is  proportional  to  the  ensemble  average  of  NBU.  Since  B 
and  q  are  statistically  independent  we  have 


Nmu- 


Cu 


g~\P>)N  o 


Bq 2 


which  is  Eq.  (3.5.2. -2) 
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5  AUTOAMATIC -REPEAT -REQUEST  (ARQ)  TRANSMISSION  over  MBC 

5.1  Performance  MRAiaCftR 

Three  performance  measures  were  used  to  characterize  a  meteor  burst  channel: 
duty  cycle,  throughput,  and  waiting  time.  Duty  cycle  is  a  function  of  system 
parameters  only,  but  throughput  and  waiting  time  are  strongly  dependent  on  the 
chosen  communication  protocol. 

5.1.1  Throughout  Rate 

We  will  examine  two  ARQ  protocols  for  point-to-point  data  transmission 
over  a  meteor  burst  channel,  similar  to  those  described  in  [36].  The  main 
difference  is  in  the  channel  behavior,  since  it  is  not  continuous  but  rather 
intermittent.  When  stop-and-wait  strategy  is  used,  transmitting  terminal  has  to 
wait  for  the  receiver  response  not  only  because  of  the  processing  in  the  receiver, 
that  is  done  in  the  idle  time  between  trails.  The  main  cause  of  the  delay  is 
the  fact  that  a  receiver  must  wait  for  the  next  available  trail  to  send  back 
its  acknowledgment,  after  it  is  done  with  checking  the  validity  of  received 
packet.  In  a  selective-repeat  strategy,  the  communication  is  full-duplex  only 
during  the  burst  duration.  Forward  channel  is  used  for  the  transmission  of  new 
and  repeated  packets,  and  along  feedback  channel  goes  the  information  about  the 
numbers  of  packets  that  have  to  be  repeated.  In  the  idle  time  between  the  trails, 
communication  stops  in  both  channels. 

In  a  stop-and-wait  strategy,  the  transmitting  terminal  sends  one  data 
packet  per  burst.  The  opening  of  a  channel  (i.e.  the  beginning  of  a  burst)  is 
indicated  either  by  a  probe  signal  from  the  receiving  terminal  at  the  beginning 
of  a  session  or,  later  on,  by  acknowledgment  signals.  Data  packet  is  N-bits  long, 
consisting  of  k  information  bits  and  r  bits  for  error  detection.  After  sending 
a  packet  the  transmitter  waits  for  the  receiver  response.  Upon  examining  the 
incoming  packet  the  receiver  sends  back  an  s-bit  positive  (ACK)  or  negative  (NAK) 
acknowledgment  signal  during  the  next  available  burst.  The  transmitter  repeats 
the  current  packet  if  it  received  the  NAK  signal  or  proceeds  to  a  new  packet  of 
data  in  case  of  the  ACK  signal.  The  repeated  or  the  new  packet  are  sent  immediately 
upon  reception  of  the  ACK/NAK  signal. 

The  throughput  rate  of  such  a  strategy  is  defined  as  the  ratio  of  information 
bits  to  total  number  of  bits  [37] 
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N-r 

(N  ♦  T  +  s)  •  E  ' 


(S. 1.1-1) 


where  t  denotes  the  number  of  idle  bit  times  between  the  last  bit  of  a  data 
packet  transmitted  and  the  first  bit  of  an  ACK/NAK  signal  received.  G  represents 
the  expected  number  of  transmissions  of  one  packet  accumulated  to  the  moment  the 
transmitter  proceeds  to  a  new  packet  of  data. 


The  probability  A  that  a  given  transmission  of  a  packet  at  the  same  time 
happens  to  be  the  last  one  is 


&  -  Pc  +  P. 


(5.1.1  -2) 


where  we  assumed  that  ACK/NAK  signal  is  received  correctly  with  probability  equal 
to  1.  (In  practice,  this  assumption  can  be  closely  achieved  by  heavy  coding  of 
acknowledgment  signals.  These  signals  are  short  and  even  with  long  codes  they 
would  occupy  negligible  amount  of  a  burst  time.  Furthermore,  they  are  sent  at 
the  beginning  of  the  burst  when  bit  error  probability  is  at  its  minimum.  If  the 
burst  time  is  critical,  the  pilot  signal  itself  can  carry  acknowledgment 
information.  )  P.  is  the  probability  that  a  received  packet  contains  an  unaecectable 
error  pattern.  For  an  average  linear  (n,k)  code  it  is  bounded  by  [38] 

P.  <  2'(n  -  (1  -  Pe) 

or,  for  a  double-error-correcting  primitive  BCH  code  by 

P  <  2'<n  -  k) 

which  means  that  P,  can  be  neglected  as  compared  to  Pc,  the  probability  that  a 
packet  is  correct,  for  any  practical  packet  length.  The  same  conclusion  is  valid 
for  the  case  when  cyclic  redundancy  check  (CRC)  polynomials  are  used  for  error 
detection.  Thus,  we  can  safely  continue  the  derivation  with  Eq.  (5. 1.1-2)  reduced 
to 


A  -  Pc.  (5.1 .1  -3) 

Now,  the  probability  that  we  need  exactly  n  -  i  transmissions  of  a  packet  to 
receive  it  correctly  is 
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0  -  A,  fjd-A,). 
*-0 


(5. 1.1-4) 


and  the  expected  number  of  transmissions  becomes 


£  -  i  ft(l-At) 

1-  1  *-0 

-i  fid-'.). 


(S.l  .1  -5) 


Averaging  E  over  all  possible  values  for  underdense  and  overdense  electron  line 
densities  qx  and  q2,  respectively,  and  decay  constant  b,  we  obtain  from  Eq. 
(5. 1.1-5) 


£  -  t 

i- 1 


(5.1  .1  -6) 


since  the  probabilities  Pe  in  different  trails  are  statistically  independent  and 
have  identical  averages .  Performing  the  summation  of  a  geometric  series  in  Eq . 
(5. 1.1-6),  the  average  number  of  packet  transmissions  becomes 

£  -  J-.  (5.1 .1  -  7) 

P< 

On  the  other  hand,  the  probability  that  a  packet  is  received  correctly 
during  a  particular  try  may  be  calculated  as  the  total  probability 

Pc  mP{C  |  underdense  transmission}  x  P(underdense  transmission )  + 

+  P{C  |  overdense  transmission)  x  P{overdense  transmission} 

-  PcttP .  +  P"P.-  (5.1.1  -8) 


The  probability  that  a  packet  is  transmitted  over  an  underdense  trail  is  the 
ratio  of  the  expected  number  of  underdense  trails  to  the  total  expected  number 
of  trails  (during  some  observation  time  TD>) 
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Similarly,  the  probability  that  a  packet  is  transmitted  over  an  overdense  trail 
is 


(5.1  .1  -  10) 


The  new  statistical  parameter  tr  is  defined  in  Eq.  (5. 1.2-3).  Averaging  Pc  from 
Eq.  (5. 1.1-8)  over  qx,  q2 ,  and  b, 

p  -  P  P  *  P  P  , 

*C  *  CU  *  <1  1  CO*  o' 

and  substituting  the  result  into  the  Eq.  (5. 1.1-1),  the  throughput  rate  follows 
as 


1 


i 


TSu 


N-r 
N  *  r  +  s 


r  p  p  *  p  p  i 

l  *  cu  '  m  *  co  *  o  J 


(5.1  .1  -11) 


Dividing  the  numerator  and  denominator  in  Eq.  (5.1.1-12)  with  the  bit  rate  R  we 
obtain 


Tp-Tr 


T p*  T ,dl,  +  T , 


ip  p  +  p  p  -j 
l*  tu1  «  ico/oJ‘ 


(5.1 .1  -  12) 


where  Tp,  Tg,  and  Tr  are  the  durations  of  a  packet,  ACK/NAK  signal,  and  parity 
bits,  respectively. 

Idle  time  Tldl,  is  a  random  variable, 

—  *  “  Tp  , 

with  x  the  time  distance  between  consecutive  bursts.  With  respect  to  the 
probability  density  function  for  x  [39] 
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we  find  the  average  value  for  Tldl.  as 


Tia„  -  f(x  -  T p)f(x)dx  -  I,  -  7> 
y  0 


(5.1  .1  -  13) 


Combining  Eqs.  (5. 1.1-9),  (5.1.1-10),  (5.1.1-12),  and  (5.1.1-13),  the  final 

expression  for  throughput  rate  is 


Tr 


Tsu  =  t,*T, 


l,  -  t,  - 

—  p  +  —  p 
t  cu  i  e 

‘11/  lIO 


(5.1  .1-14) 


when  stop-and-wait  strategy  is  used. 

In  a  selective- repeat  strategy,  the  transmitter  sends  as  many  data  packets 
as  the  duration  of  a  current  burst  allows.  The  transmitter  determines  the  burst 
duration  with  the  help  of  a  feedback  channel,  which  is  continuously  used  for 
reception  of  acknowledgment  information.  If  the  receiver  buffer  is  infinite  and 
if  we  assume  no  errors  in  the  feedback  channel,  the  throughput  is  equal  to  the 
probability  that  the  packet  is  received  correctly  [40] , 

Ts,-Pc,  (5.1.1-15) 


representing  the  average  number  of  packets  communicated  (i.e.  transmitted  and 
received  correctly)  per  transmission.  The  throughput  is  obtained  as  the 

capacity  of  an  M-ary  erasure  channel,  where  M  is  the  number  of  possible  choices 
for  a  packet  [41] .  Hence  the  selective-repeat  strategy  represents  the  optimal 
form  of  ARQ  when  the  receiver  buffer  size  is  unlimited. 

We  can  consider  Pc  to  be  a  discrete  random  variable  taking  values 

Pe  -  Peu  ,  during  underdense  bursts 

-  Pco  ,  during  overdense  bursts 

-  0  ,  elsewhere  . 

Hence,  its  expected  value  will  be 
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E{  Pe }  -  Peu  P{Pe  -  Pcu)  +  PC0  P(Pe  “  Peo>  ■ 

Now,  the  probability  that  Pc  is  equal  to  Pcu  is  the  ratio  of  the  total  underdense 
burst  time  to  total  observation  time 


P{PC=P  cu> 


T_£  i 

>,u  BU  <BU 

— K  — 

T  o  /iu 


(5.1  .1  -  16) 


Total  burst  time  is  a  product  of  the  expected  number  of  bursts  and  their  average 
duration.  The  mean  time  between  trails  is  always  greater  than  the  average  burst 
duration,  in  underdense  and  in  overdense  case.  Similarly,  the  probability  that 
Pc  is  equal  to  Pco  is  found  to  be 


P<P'=P  „> 


1  BO 
*  10 


(5.1  .1  -  17) 


Using  Eqs .  (5.1.1-16)  and  (5.1.1-17),  the  expression  for  the  average  throughput 
of  a  selective- repeat  system  becomes 


Ts. 


ip  +  ‘jo 

i  c“  t 

1 IU  ‘10 


(5.1  .1  -  18) 


The  question  remaining  to  be  solved  is  the  evaluation  of  average  values 
for  Pcu  and  Peo.  The  beginning  of  the  packet  can  fall  into  time  instants  separated 
by  a  packet  length  Tp,  counting  from  the  beginning  of  a  trail.  In  that  sense, 
Pcu  is  found  as  the  total  probability 

Pcu  -  I  n~7T~\  p“<~kT»')-  (5.1.1-19) 
t-o  1 1  >v'  1  p  J 

where  we  assumed  that  packet  positions  are  uniformly  distributed  inside  the 
trail.  Notation  [x-Jsignifies  the  integer  less  than  or  equal  to  x.  In  that  context, 
the  ratio  tBU/Tp  represents  the  number  of  packets  that  will  fit  into  the  burst 
duration.  Peu(kTp)  is  the  probability  that  a  packet  in  the  k-th  time  slot  will 
be  received  correctly.  Bar  on  probabilities  Pcu  and  Pco  indicates  averaging  over 
electron  line  densities  and  a  decay  constant.  Expression  identical  to  the  one 
in  Eq.  (5.1.1-18)  describes  the  situation  with  the  overdense  trails,  if  tBU  is 
replaced  with  tM. 
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5.1.2  Waiting  Tims 

Waiting  time  to  deliver  a  data  packet  is  defined  as  the  time  elapsed  from 
the  moment  a  transmitter  has  a  packet  ready  for  sending,  to  the  time  instant 
the  receiver  starts  sending  back  a  positive  acknowledgment  signal. 

Underdense  and  overdense  trail  arrivals  are  independent  Poisson  processes, 
with  respective  parameters 

\u  -  and  k0  -  (5.1.2-  1) 

1 IV  1 10 

The  random  process  consisting  of  all  trail  arrivals  is  also  Poisson  with  the 
parameter 

k  ‘  ku  +  k0.  (5. 1.2-2) 


what  says  that  the  mean  time  between  bursts  in  this  new  process  is 


1  hu  (io 
k  tw  *  1 10 


(5.1  .2-3) 


When  trail  arrivals  are  modeled  as  Poisson  random  points  tlt  the  waiting  time 
t*  for  stop-and-wait  strategy  is  given  by 

tw  “  Xn+l  “  t-n+i  -  t0  , 

where  x„+1  is  the  time  distance  from  the  fixed  beginning  of  the  observation  time 
t0  (i.e.  the  instant  a  transmitter  has  a  packet  ready  for  sending,)  to  the  (n  + 
l)th  random  Poisson  point  to  the  right  of  t0.  Namely,  if  the  number  of  transmissions 
required  to  receive  a  packet  correctly  is  equal  to  n,  we  must  wait  to  the  (n  + 
l)th  trail  for  the  receiver  confirmation.  Random  variable  xI)+1  has  Erlang  prob¬ 
ability  density  function  [39] 


/..,(*)  - 


and  averaging  over  x  leads  to 
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n  *  1 
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Averaging  further  over  n,  a  discrete  random  variable  defined  by  Eq.  (5. 1.1-4), 
the  expected  waiting  time  becomes 


tv  =  (1  *E)  I, 


(5.2.2-  4) 


The  minimum  possible  waiting  time  is  achieved  in  the  case  when  every  packet  is 
correctly  received  during  its  first  transmission.  Expected  number  of  transmissions 
E  is  then  equal  to  1,  and  the  resulting  lower  bound  on  waiting  time  to  receive 
a  correct  packet  in  a  stop-and-wait  strategy  is 


(5.1  .2-5) 


a  function  of  meteor  burst  parameters  only. 

5.1.3  Duty  Cycle 

The  duty  cycle  of  a  meteor  burst  communication  system  is  defined  as  the 
percentage  of  time  that  the  received  signal  is  above  some  arbitrary  threshold, 
what  is  equivalent  to  the  condition  that  the  trail  electron  line  density  be  above 
some  minimum  value  q,^,,.  It  can  be  found  as  the  ratio  of  the  total  burst  time 
within  an  observation  time  TD,  to  the  observation  time  TD  itself.  On  the  average, 
the  total  burst  time  is  the  product  of  the  expected  number  of  trails  during  the 
observation  time  and  their  mean  durations.  Taking  into  account  both  types  of 
trails 


duty  cycle  - 


1  BU 

f|U  ^10 


(5. 1  .3-  1  ) 


where  mean  values  for  underdense  and  overdense  burst  durations,  tBU  and  tgo, 
contain  the  information  about  the  receiver  threshold  (i.e.  qmin.)  Eq.  (5. 1.3-1) 


assumes  that  trails  do  not  overlap.  If  that  happens  to  be  the  case  it  is  an  upper 
bound . 

5.2  Meteor  Burst  Communication  System 

Digital  data  are  transmitted  over  a  meteor-burst  channel  using  binary 
phase-shift  keying  modulation,  at  a  constant  data  rate  R.  The  transmitted  signal 
is 


s(l)  =  ±'j2PTcosa)0t ,  0<t<Tb. 

where  PT  is  the  transmitter  power  and  f0  -  2n/Lu0 is  the  carrier  frequency.  Tb  - 
1/R  is  the  bit  duration. 

Multipath  phenomena  are  encountered  in  meteor -burst  communications  only 
at  high  data  rates;  there  appear  very  few  underdense  trails  with  electron  line 
density  high  enough  to  provide  sufficient  scattered  signal  power, and  overdense 
trails  become  dominant  propagation  medium.  At  low  data  rates  underdense  trails 
are  much  more  numerous,  they  act  as  small  coherent  sources  and  scattered 
continuous -wave  signals  exhibit  good  space  and  frequency  correlation.  The  phase 
of  the  received  signal  is  assumed  constant,  representing  only  the  gross  time 
delay  due  to  the  average  distance  which  the  signal  must  propagate.  Such  delay 
amounts  only  to  a  shift  in  time  origin. 

The  signal  received  from  an  underdense  or  overdense  trail  is  of  the  form 

r(()  =  ±Acos<x>0t  +  n(t ) 

=  ±v'2/pJcosuo0(  *  n(t),  0  <  t  <  T  b , 

where  A  is  a  random  amplitude  and  n(t)  is  a  Gaussian  noise  with  one-sided  power 
spectral  density  defined  by  Eq.  (2.4-1),  BPSK  signals  have  the  same  energy  Eb  - 
P,Tb ,  and  the  optimum  decision  threshold  is  set  at  Vt  -  0  volts.  Since  the 
attenuation  of  the  signal  Is  always  positive,  the  optimum  decision  regions  are 
invariant  to  radial  scaling  of  received  signals.  Under  these  conditions,  a 
correlation  or  matched- filter  receiver  structure  is  still  optimum,  completely 
independent  of  the  fact  whether  amplitude  A,  or  probability  density  function  of 
amplitudes  p(A)  is  known  precisely  [42]. 


The  received  signal  amplitude  is 


Ik 


A  =  j2P~t. 


For  underdense  trails,  the  received  power  P,  is  equal  to  P,u  from  Eq.  (2. 1.1-2), 
a  function  of  two  random  variables ,  qj  and  b  .  Transforming  their  density  functions 
from  Eqs  .  (2 . 1 . 2-2)  and  (2 . 1 . 3-  3)  we  obtain  the  probability  density  for  underdense 
signal  amplitudes  A„ 
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(5.2-  1  ) 


that  is  valid  in  the  interval 


\2Cuqrwnu  <  Au  <  v2 Cuqr 


The  integration  limits  in  Eq .  (5.2-1)  are 
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For  overdense  trails,  the  received  power  P,  is  equal  to  Pso  from  Eq .  (2. 1.5 -2), 
a  function  of  one  random  variable  q2 .  Transforming  its  distribution  from  Eq . 
(2. 1.6-2),  probability  density  for  overdense  signal  amplitudes  A0  becomes 
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valid  in  the  interval 
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Both  probability  density  functions  are  too  complicated  as  to  gain  some  useful 
insight  into  the  detection  process. 
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5.2.1  Pncoded  Transmission 

5 . 2 . 1 . 1  Transmission  over  Pnderdense  Trails 

The  probability  of  the  i-th  bit  in  a  packet  being  in  error,  expressed  in 
terms  of  received  power  P,u  from  Eq .  (2. 1.1-2),  is 

1 

2 

(5. 2. 1.1-1) 
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when  BPSK  modulation  is  used.  We  took  the  bandwidth  to  be  equal  to  the  bit  rate 
R  (1  bit/sec/Hz)  .  The  time  within  the  trail  has  discrete  values  t  -  i/R.  To 
simplify  the  analysis,  we  assumed  that  Psu  remains  constant  during  the  bit  interval 
and  equal  to  the  power  at  the  end  of  the  interval.  This  is  a  'worst  case'  approach 
on  a  small  scale,  since  the  actual  S/N  ratio  for  every  bit  will  be  better  than 
the  one  used  in  calculation.  The  probability  that  the  packet  N-bits  long  is 
correct  is  then 
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CU 


(5.2.1  .1-2) 


with  the  average  value 


Pcu  “  j  J  Pcu' /*(&)•  nq^dq.db, 

b-Oq,-q  mum 


(5.2. 1  .  1  -  3) 


In  addition,  the  variable  from  Eq .  (5.1.1-19)  is  evaluated  as 


I- 1 


1  -fftul  J  +  kT  P 


and  its  average  is  calculated  in  the  same  way  as  in  Eq .  (5. 2. 1.1-3). 
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5. 2. 1.2  Transmission  over  Overdense  Trails 

The  bit  error  probability  expressed  in  terms  of  received  power  P,0  from  Eq. 
(2. 1.5-2)  is  now 


i 


with  all  the  underlying  assumptions  stated  for  Eq .  (5. 2. 1.1-1).  Accordingly,  the 
probability  of  a  packet  being  correct  in  these  conditions  becomes 
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CO 


n 

i- 1 


(5.2.1  .2-2) 


with  the  corresponding  average  value 

Q  no 

p  co  “  J  Pco- f(qz)ctq2.  (5. 2. 1.2-3) 

9  2  "  ^  nuo 

5-2.2  Coded  Transmission 

For  coded  transmission,  convolutional  codes  were  chosen  because  of  their 
superior  performance  compared  to  block  codes,  for  the  same  implementation 
complexity  of  the  encoder-decoder  [43]. 


We  did  not  consider  the  class  of  burst-error  correcting  codes,  since  very 
few  statistical  data  are  available  on  the  length  of  meteor-burst  channel  memory. 
Although  exponential  decay  of  received  power  toward  the  end  of  the  burst  in 
underdense  case  leads  to  conclusion  that  bit  errors  become  more  frequent,  the 
noise  model  from  Eq.  (2. 1.4-1)  shows  that  for  narrowband  modulation  the  noise 
may  be  considered  white,  and  noise  samples  in  successive  bit  intervals  are 
consequently  independent. 


For  higher  data  rates  or  for  spread  spectrum  systems  noise  can  not  be  taken 
as  white  any  more,  since  its  power  spectral  density  will  change  inside  the 
signaling  bandwidth .  Noise  samples  become  dependent,  and  this  noise  model  predicts 
the  appearance  of  channel  memory.  If  experimental  data  confirm  these  conclusions, 
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use  of  interleaving  will  become  necessary  to  "randomize"  transmission  errors  and 
maintain  the  performance  of  convolutional  codes,  which  are  not  best  suited  for 
burst  errors. 

We  considered  both  rate  1/n  and  rate  k/n  convolutional  codes.  Rate  1/n 
codes  provide  for  superior  error-correcting  performance,  at  the  cost  of  sub¬ 
stantial  lengthening  of  the  original  packet.  For  this  reason,  codes  with  n  greater 
than  2  are  at  disadvantage  in  meteor  burst  communications,  since  the  coded  packet 
length  becomes  prohibitively  large  and  will  often  exceed  burst  duration,  resulting 
in  the  loss  of  many  bits  as  their  time  slots  fall  beyond  the  receiver  threshold 
point.  Rate  k/n  codes  are  less  powerful  but  more  compact,  and  their  error  - 
correcting  capability  is  not  easily  overwhelmed  by  the  strong  limiting  factor 
of  relatively  short  burst  durations,  particularly  when  k/n  ratio  approaches  to 
1.  In  our  work,  code  rates  of  3/4,  7/8,  and  15/16  were  considered. 

When  decoding  convolutional  codes,  the  error-correcting  capability  is 
difficult  to  state  pre  sely.  With  maximum  likelihood  decoding,  a  code  can  correct 
t  errors  within  3  to  5  constraint  lengths,  where  t  is  defined  as 

.  -  . 

L  2  J 

in  terms  of  the  free  distance  d£.  The  exact  length  depends  on  actual  error 
distribution,  and  for  a  particular  code  and  error  pattern,  transfer  functions 
methods  are  used  to  obtain  bounds  [43],  In  our  case  the  problem  is  further 
aggravated  by  the  fact  that  signal- to-noise  ratio  is  not  constant  throughout  the 
packet,  but  rather  changes  as  the  received  power  varies  with  time. 

One  way  to  overcome  these  difficulties  is  to  use  the  notion  of  cutoff  rate 
R0  for  binary  transmission 

R0=  1  -  iog2(l  +  Z?)  bits/symbol  , 

in  order  to  decouple  the  influence  of  the  coding  channel  characterized  by  R0 
from  the  coding  technique  [44],  For  example,  in  case  of  rate  1/2  code  with  K  - 
7  and  hard  decision  decoding,  the  decoded  bit  error  rate  is  bounded  by  [44] 

P  b  coded  *  1 8  D 1 0  -  1 05 . 5  D  12  *  702  D  1 4  +  58  1 6 . 5  D  " 6  +  . . .  (5.2.2-  1) 


where 


D  =  2 yfPtil-K)  .  (5. 2. 2-2) 

and  Pb  is  the  channel  bit  error  rate. 

For  rate  k/n  codes,  we  used  the  following  bound  on  decoded  bit  error 
probability  [45] 
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(5. 2. 2-3) 


where  df  is  the  free  distance  of  the  code,  E,  -  Ebr  is  the  symbol  energy,  r  is 
the  code  rate,  E*,  is  the  bit  energy,  and  N0  is  the  noise  power  spectral  density. 
Code-specific  coefficients  Wj  are  summarized  in  tables  for  different  values  of 
the  code  constraint  length  [46], 

Inf inite- level  quantization  results  in  2.2  dB  improvement  over  two -level 
quantization.  It  can  be  substituted  by  the  3-bit  soft  decisions  of  the  channel 
output,  resulting  in  approximately  a  2  dB  gain  over  the  hard  quantized  binary 
symmetric  channel. 

Bounds  on  decoded  bit  error  rate  are  derived  for  constant  signal- to-noise 
ratio  and  arbitrarily  large  path  memory,  i.e.  the  depth  of  the  input  bit  history 
stored  by  the  decoder.  But,  it  was  shown  in  [47]  that  a  fixed  amount  of  path 
history,  namely  4  or  5  times  the  constraint  length,  is  sufficient  to  limit  the 
degradation  from  the  optimum  decoder  performance  to  about  0.1  dB  for  the  BSC 
and  Gaussian  channels.  In  a  meteor  burst  channel  signal - to-noise  ratio  is 
constantly  changing,  but  most  of  the  time  it  does  not  change  substantially  during 
the  fixed  decoding  delay,  and  we  took  it  to  be  constant  for  the  bits  along  one 
path  memory  length.  Furthermore,  the  decision  about  a  specific  bit  is  brought 
by  analyzing  the  trellis  branches  consisting  of  bits  that  all  have  equal  or 
higher  signal-to-noise  ratios  than  the  bit  in  question. 
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5. 2. 2.1  Transmission  over  Pnderdense  IraUs 

Parameter  D  from  Eq.  (5. 2. 2-2)  in  this  case  becomes 

D  =  2jPbu(\-Pbu)  , 


when  we  use  the  values  for  Pbu  from  Eq .  (5. 2. 1.1-1).  After  calculating  the  values 
for  Pbu  cod#d  according  to  Eqs  .  (5. 2. 2-1)  or  (5. 2. 2-3),  the  probability  that  encoded 
packet  will  be  correctly  received  is  found  as  the  product 


/V 

-  n 


i  -  p 


bu  coded 


(*)]■ 


Its  average  value  is  then 

«  Quu 

^cucod.d  =  /  /  ^cucod.d-/6(b)’  f{q,)dq,db.  (5. 2. 2. 1-1) 

5 . 2 . 2 . 2  Transmission  over  Overdense  Trails 

Following  the  same  steps  as  in  Sec.  5. 2. 2.1.  we  obtain  the  parameter  D  in 
the  form 


D  =  2jPbo{\~Pbo)  . 


the  probability  that  a  coded  packet  is  correctly  received  as 
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and  its  corresponding  average  value 


co  coded  J  ^  co  coded  ■  f(q2)dq2. 

2  "  Q  rnino 


(4.2.2-  1  ) 
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6  CONCLUSION  AND  IDENTIFICATION  OP  FUTURE  WORK 


In  this  work  a  complex  stochastic  model  for  the  meteor  burst  channel  was 
used  to  serve  as  the  bedrock  for  the  analysis  of  two  different  communication 
systems.  The  first  system  assumed  we  are  transmitting  with  a  constant  bit  rate 
while  the  BER  is  time  varying  yet  constrained  to  be  below  some  maximum  level. 
The  second  communications  system  employs  a  variable  bit  rate  whose  adaptive 
behavior  is  controlled  so  as  to  maintain  constant  BER.  Both  systems  operate 
under  the  control  of  a  communication  feedback  protocol  that  provides  knowledge 
of  channel  arrival  (the  instant  at  which  the  channel  first  becomes  available) 
and  monitor  the  strength  (power  received)  of  the  channel.  The  research  herein 
focused  on  throughput  and  related  parameters.  The  average  throughput  was  shown 
to  depend  upon  the  arrival  rate  of  the  bursts  and  the  expected  number  of  bits 
per  burst  as  averaged  over  the  ensemble  of  (random)  time  functions  that  describe 
its  behavior.  In  deriving  the  expressions  for  the  above  it  was  necessary  to 
derive  the  relationship  between  the  number  of  bits  per  given  burst  and  bit  rate. 
This  in  turn  requires  understanding  of  channel  duration  -  the  available 
transmission  time  -  for  a  given  channel  (burst) . 

For  constant  bit  rate  system  the  channel  duration  decreases  with  increasing 
bit  rate.  Since  the  number  of  bits  transmitted  on  the  burst  is  given  by  the 
product  of  burst  duration  and  bit  rate  a  clear  tradeoff  exists  between  bit  rate 
and  transmission  time.  This,  naturally,  implies  the  existence  of  an  optimal  bit 
rate  that  would  maximize  the  number  of  transmitted  bits  for  the  burst.  Having 
chosen  an  optimal  bit  rate  for  the  burst  at  hand  we  must  repeat  the  process  for 
the  next  burst  whose  behavior  is  different.  In  fact,  to  optimize  the  throughput 
we  must  change  the  bit  rate  from  one  burst  to  the  next.  Having  to  change  the 
bit  rate  on  a  per  burst  basis,  however,  requires  knowledge  of  the  burst  time 
behavior.  The  problem  with  this  last  requirement  is  twofold;  first,  a  priori 
knowledge  of  the  burst  is  not  available  and  second,  estimation  of  the  burst  based 
on  assumed  ideal  behavior  is  unreliable  since  the  bursts  often  deviate  from  such 
assumed  bahavior.  To  counter  this  limitation  we  derived  the  relationship  between 
the  constant  bit  rate  and  the  expected  average  number  of  bits  per  burst  since 
the  latter  is  directly  proportional  to  throughput.  We  then  find  the  constant 
bit  rate  (for  the  system)  which  will  maximize  the  throughput.  This  bit  rate  is 
fixed  over  all  bursts  and  therefore  it  is  a  system-optimal  bit  rate.  Having 
found  the  maximum  average  number  of  bits  per  burst  we  multiply  by  the  arrival 
rate  of  burst  and  get  the  maximum  throughput  expression  in  terms  of  system 


81 


parameters  and  constraints.  The  significance  of  such  expression  is  its  utility 
in  analyzing/improving  existing  systems  and  in  designing  future  constant  bit 
rate  meteor  burst  systems. 

In  variable  bit  rate  systems  we  change  the  bit  duration  (on  a  bit  by  bit 
basis)  so  the  bit  energy  remains  constant  corresponding  to  our  desired  fixed 
BER.  The  bit  energy  is  given  by  the  time  integral  of  the  power  received  function 
over  the  duration  of  the  bit.  Finding  the  individual  bit  durations  is  tantamount 
to  dividing  the  power  received  curve  into  equal  strips  each  of  area  equal  to 
the  bit  energy.  As  discussed  earlier  a  priori  knowledge  of  the  burst  time 
function  is  not  available  and  estimation  is  at  best  problematic.  This  precludes 
exact  solution  of  the  bit  duration.  The  best  we  can  do  is  determine  the  bit 
duration  based  on  the  bit  power  at  the  beginning  time  of  the  bit.  This  may  yield 
an  optimistic  (too  short)  bit  duration  if  the  power  decays  during  the  duration 
of  the  bit.  The  energy  per  bit  for  such  a  case  would  be  below  the  desired  fixed 
one.  If  the  power  is  monotomically  decreasing  there  will  be  a  bit  whose  bit 
energy  relative  to  the  assumed  bit  energy  will  be  unacceptable.  The  restriction 
on  maximum  bit  duration  in  turn  constrains  theoretically  yet  not  practically, 
as  has  been  shown,  the  total  time  of  transmission  for  the  burst.  The  number  of 
bits  transmitted  during  a  burst  can  be  approximated  by  the  area  under  the  power 
received  for  the  interval  given  by  transmission  time.  The  throughput  is  found 
by  multiplying  the  arrival  rate  by  the  average  number  of  bits  per  burst. 

As  it  turns  out  the  throughput  for  variable  rate  system  is  about  six  dB 
ff  higher  for  a  conservative  sample  system  .  The  interesting  point  however  is  the 

fact  that  the  improvement  of  throughput  for  variable  over  constant  bit  rate 
system  (using  underdense  bursts  only)  was  in  large  part  due  to  the  arrival  rate. 

*•  For  constant  bit  rate  system  the  high  optimal  bit  rate  precluded  many  bursts 

(those  with  electron  line  density  below  the  minimum  line  density  threshold)  from 
being  useful  resulting  in  a  low  arrival  rate.  For  variable  rate  systems  since 
the  bit  rate  can  be  quite  low  (though  limited)  implies  a  lower  level  of  cutorf 
for  useful  bursts  and  higher  arrival  rate.  The  tradeoff  lies  in  the  fact  that 
a  low  threshold  rate  for  the  adaptive  rate  increase  system  complexity. 

The  Improvement  in  throughput  however  must  be  kept  in  perspective.  Our 
analysis  assumed  that  the  constant  BER  under  variable  rate  system  is  same  as 
the  worst  BER  for  the  constant  bit  rate  system.  Consequently,  given  N  transmitted 
bits  the  probability  of  all  being  error  free  is  lower  for  the  constant  bit  rate 
then  for  the  variable  bit  rate  system. 
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The  higher  throughput  for  overdense  relative  to  underdense  bursts  is  a 
natural  outcome  of  its  longer  average  burst  duration.  The  tradeoff,  however, 
is  in  the  much  longer  inter-arrival  time  of  the  overdense  bursts  compared  with 
underdense;  a  ratio  of  15:1.  We  see  then  that  for  short  urgent  messages  where 
waiting  time  is  crucial  if  we  are  given  a  constant  bit  rate  system  it  is 
advantageous  to  operate  the  system  at  the  optimal  rate  corresponding  to  underdense 
burst  utilization. 

Two  transmission  protocols  were  investigated  in  automatic-repeat-request 
strategy:  stop-and-wait ,  when  one  data  packet  is  transmitted  per  burst,  and 
selective-repeat,  when  entire  duration  of  the  burst  is  used  for  transmission. 

Fig  23  shows  the  throughput  rate  for  stop-and-wait  strategy  as  a  function 
of  packet  length,  for  the  uncoded  case  and  rate  1/2,  3/4,  7/8,  and  15/16 
convolutional  codes.  In  coded  as  well  as  in  uncoded  transmission  there  exists 
an  optimum  packet  length  for  a  given  data  rate .  Improvement  in  throughput  is 
obtained  in  the  coded  case,  but  only  for  code  rates  close  to  1.  High  rate 
convolutional  codes  are  less  powerful  than  low  rate  ones ,  but  they  do  not 
excessively  increase  the  original  packet  length  and  still  can  correct  some  of 
the  transmission  errors.  Low  rate  codes  are  at  disadvantage  because  of  the  strong 
limiting  factor  of  relatively  short  burst  durations.  Fig.  24  shows  the  waiting 
times  necessary  to  receive  a  specific  data  packet  correctly.  For  a  given  packet 
length  there  always  exists  an  optimim  data  rate,  in  coded  and  in  uncoded 
transmission.  High  rate  codes  yield  lower  waiting  times  with  respect  to  uncoded 
case  for  all  data  rates,  and  with  respect  to  high  rate  codes  at  low  data  rates 
when  throughput  is  at  its  maximum.  Low  rate  codes  show  better  performance  in 
terms  of  minimizing  the  waiting  time  at  higher  data  rates,  but  then  the  throughput 
becomes  prohibitively  low. 

Throughput  performance  for  selective -repeat  strategy  is  summarized  in  Fig. 
25.  Here  again  high  rate  codes  are  superior  to  low  rate  ones  and  to  uncoded  case. 
There  is  no  optimum  packet  length  and  short  packets  yield  higher  throughputs, 
with  the  ultimate  maximum  obtained  for  packet  length  of  1  bit.  The  practical 
lower  limit  on  the  packet  length  is  imposed  by  necessary  overhead  bits  and  user 
considerations . 

In  this  analysis  many  parameters  considerations  and  variables  interplay. 
Amongst  them  complexity  of  model,  empirical  values,  transmitter  power  antenna 
design  link  geometry  and  geographic  location,  frequency,  noise  behavior  modulation 
bit  rate  and  BER.  All  of  the  above  have  been  integrated  and  the  derived  closed- form 
relationships  are  stated  in  the  most  general  form  incorporating  all  possible 
variables  to  facilitate  their  usage  as  design  and  analysis  tools.  Furthermore, 
these  expressions  define  and  better  delineate  the  tasks  for  future  analysis. 
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Some  of  the  issues  that  need  be  investigated  are  coding,  packet  design,  modulation 
protocol  development  and  networking.  Each  of  these  requires  extensive  research 
and,  clearly,  there  is  no  'ultimate'  meteor  burst  system  but  rather  application 
dependent  design.  With  faster  and  cheaper  hardware  for  control  and  storage  and 
renewed  focus  on  this  complex  channel  it  is  clear  that  we  are  in  the  exciting 
phase  of  using  this  natural  phenomenon  for  communication. 
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7  TABLES 


Table  1  —  Order-of-Magnitude  Estimate*  of  the  Properties  of  Spormdic  Meteor*  (1 ) 


Meteor 

Particles 

Maaa 

(g) 

Radius 

Number  of  This 
Maas  or  Greater 
Swept  Up  by  the 
Earth  Each  Day 

Electron  Line 
Density  (electrons 
per  meter  of  trail 
length) 

Particles  pass  through  the 

atmosphere  and  fall  to 
the  ground 

104 

8  cm 

10 

Particles  totally  die- 

103 

4  cm 

102 

— 

integrated  in  the 

102 

2  cm 

10s 

— 

upper  atmosphere 

10 

0.8  cm 

104 

1018 

1 

0.4  cm 

105 

1017 

iO'L, 

0.2  cm 

10s 

10l« 

nr2  . 

0.08  cm 

107 

JO15 

10*3 

0.04  cm 

10® 

1014 

io-4 

0.02  cm 

10® 

1013 

io-3 

80  jim 

1010 

1012 

Approximate  limit  of 

io-« 

40  jire 

1011 

1011 

radar  measurements 

io-7 

20  jim 

1012 

1010 

io-* 

8  jim 

? 

f 

Micrometeorities  (Particles 

io*9 

4  jim 

Total  for  this 

Practically  none 

float  down  unchanged  by 

10-i° 

2  jim 

group  estimated 

atmospheric  collisions) 

10*u 

10*12 

0.8  jim 

0.4  jim 

as  high  as  1021* 

Particles  removed  from 

10’13 

0.2  jim 

— 

— 

the  solar  system  by 
radiation  pressure 

Table  2 
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ABSTRACT 

A  model  It  described  for  the  enderdente 
meteor  burst  phenomenon  thet  tikes  Into  account 
Its  random  decay  time-constant  and  derives  the  CDF 
and  PDF  for  the  burst  duration  as  well  as  the 
Joint  PDF  for  burst  duration  and  decay  time- 
constant.  In  addition  the  time  varying  bit  error 
rate  BCR  and  Its  average  Is  calculated  for  binary 
FSK  and  coherent  IPSK. 

INTRODUCTION 

The  development  and  understanding  of  a  model 
for  the  beha vlorof  a  meteor  burst  Is  crucial  In 
the  study  and  use  of  this  phenomenon  as  a  viable 
communication  channel,  ref.  (1-7,  11).  For  1 
communication  system  engineer  It  Is  Important  to 
have  an  accurate  understanding  of  the  time  varying 
bit  error  rate.  BER,  characteristics. 

It  Is  known  that  the  SNR  of  the  underdense 
trail  decays  eaponentlal ly  In  time  with  a  random 
Initial  amplitude  and  a  random  decay  time-constant 
which  In  turn  affect  the  BER  variation  In  time. 
To  date,  however.  Investigators  have  used  either  a 
constant  -  (worst  case)  •  SNR  based  on  fixed  decay 
rate  In  the  calculation  of  the  bit  error  rate, 
ref.  (11,  13)  or  employed  a  SNR  based  on  a 
conditional  expected  value  of  decay  rate  given  a 
fixed  burst  duration,  ref.  (10). 

This  paper  takes  Into  account  effects  of  the 
random  decay  time-constant  In  Its  derivation  of 
BER  expression.  Using  empirically  reasonable 
distributions  for  Initial  SNR  amplitude  and  decay 
constant  we  derive  the  cumulative  distribution 
function,  CDF,  of  the  burst  duration,  tg,  as  well 
as  the  Joint  PDF  of  decay  time-constant  and  bur»t 
duration.  The  SNR  Is  then  expressed  In. terms  of  a 
minimum  SNR  level  n  ,  decay  time-constant, 8,  and 
burst  duration,  tg.  This  SNR  expression  Is  used 
In  obtaining  the  appropriate  expression  for  BER 
and  average  BER  In  time  for  thg  two  cases  of 
binary  PSK  and  coherent  BPS*. 

THE  MODEL  FOR  METEOR  BURST 

Meteor  bu-sts  are  charaeterlted  as  overdense 
or  underdense.  Since  the  underdense  trail 
phenomenon  Is  the  dominant  contributor  to  the 
overall  throughput  we  shall  consider  Its  model, 
ref.  (5). 

•Acknowledgement.  This  work  Is  supported  by 
AFOSR  Grant  No,  ?7747  and  ay  SCS  Telecom,  Inc. 


According  to  ref.  (3)  wt  can  write  the  SNR  for 
underdense  trill  is: 

sn  .  t-t/i 

(fa^tftjd-tli/feoa’dNF 

Where, 

Pj  Is  the  transmitter  power 

Gg  Is  the  receiver  antenna  gain 

Gf  Is  the  transmitter  antenna  gain 
X  Is  the  wavelength  of  the  carrier 
0  Is  the  diffusion  coefficient  of  the 
atawsphere  (7.5  ar/s) 

g  Is  the  angle  of  Incidence  of  the 
transmitted  plane  wave 

B  Is  the  angle  between  the  trail  and  the 
great  circle  path  from  receiver  to 
transmitter 

Rj  Is  the  nominal  distance  between  the 
trail  and  the  stations 

ra  Is  the  nominal  Initial  radius  Of  a 
trail  (0.65m) 

<»,  fs  the  effective  echoing  area  of  the 
electron  (10‘ZB  mz) 

HF  Is  the  background  noise  floor 
B  Is  the  decay  time-constant  of  the  burst 

a  Is  the  angle  between  Incident  E-vector 

at  trail  and  the  direction  from  the 
tra 1 1  to  the  receiver 
R  Is  the  electron  line  density 

Further  simplifications  are  according  to  ref  (10). 
NF  Is  computed  using 

NF  •  kleH  Ft 

when  k  Is  Bolt? man's  constant,  T.  Is  ( ? 7 0  k),  W  Is 
the  receiver  bandwidth  and  Fj  Is  the  external 
noise  figure  •  assumed  to  be  due  to  galactic  noise 
only. 

n  and  p  are  assumed  to  be  uniforms! ly  distributed 
over  (-■/?,  »/7)  and  hence  <a  s  In  >  *  <  s  in  B>  •  . 

Also  It  Is  assumed  fori  (the  angle  of  Incidence 
for  the  transmitted  wave)  that  all  trails  are 
located  at  1001m  altitude  and  halfway  between  the 
transmitter  and  receiver. 

for  athamatlcal  tractablllty  we  rewrite  eg  (1)  as 
fol lows : 


SM(t)  •  AqVt/B 

where  q  tlx  electron  11m  density  varies  between 
Ri  ini  4U  the  alnlaua  and  aexlaua  « 1  ectroa  llm 
dtniltlti.  Me  dlflOtl  by  n  the  alnlaua  SNA  level 
iWthst  it  find  tht  karst  end  (duration)  sad 

corresponds  to  A  . 

0 

♦  ata 


SftxOg) 


•  A*  (J)  yields  an  expression  for  the  Initial  SNA 
torsi.  SNA(t*o). 


which  whon  substituted  Into  eq  (?)  provides  ut 
with  the  slgnal-to-nolse  ratio  In  teras  of  the 
alnlaua  SNR  level.  .  burst  duration,  tg.  and 
decay  tlae  constant,  (. 


SNR  ( t ) 
COF  AND  POF  OF  t. 


(t,  -  t)/8 


Throughout  this  p»ptr  we  shill  use  the  POT  of  q 
is: 

<  i.t  .  -3  V,,4N  (6) 


Where 

a  .  %  \  and  r”1 

'  VA  \ 

This  yelld  the  CDF  of  q 

r, <,)<•-?>  0) 

In  addition  we  shall  assume  the  PDF  of  B  as 
exponentially  distributed 

f,(S)  •  <J^1  e',/*1  0L  -  «  9,9 

or  Rayleigh 

i|H)  •  !,  *  ■•*/*<>*  »  ;  o  (9) 

r 


and  8,  -  .35  sec.  8,  Is  the  expected  value  of  8 
and  lls  value  used  is  from  the  most  recent 
observations  by  AAOC.  (The  most  extensive 
program  for  the  study  of  aeteor  burst  phenomenon 
Is  currently  conducted  by  the  Air  Force  (RADC), 
ref  (6)). 

Using  eg.  (7)  to  represent  SNR  behavior  In  l\mt) 
and  egs.  (3)  end  (7),  setting  q  •  nL  11  ,c,n  **  ’ 
shown  that  the  Joint  POF  of  t8  and  0  Is  given  by 
eq.  (10) 


1 

Vi  ltB.»>  •  15  •  * 

B  (10) 

and  by  Integrating  over  I  we  get 

•i,  'ls'  *  I  h  •  **  ft  !*•  *  (11) 

We  are  left  now  to  substitute  eq.  (8)  In  eqs. 
(10)  and  (II)  resulting  In  the  joint  RDF  of  t| 
and  8  (8  -  exponents  I ) 

1  *  ;•,) 

ft,  ,  (,a.i’  *  hf  I  • 
the  PDF  of  tB  (8  -  exponential) 

'1,‘V  *  t  V,}  lB  >°  (,3> 

where 

and  K,(.)  Is  the  aiodl  fled  llankel  function  of 
order  I. 

and  the  COF  of  tg 

Ft  (tB)  •  1  *  i*,(0  t,>0  (ia) 

B 

eqs.  (13)  and  (Id)  are  compared  graphically  to 
an  exponential  distribution  with  E(tg)  •  .37 / 
sec.  and  depicted  In  fig.  1. 

Slail  1  »r! y  we  repeat  the  process  by  substituting 
eq.  (9)  (8  Rayleigh)  into  eqs.  (10a)  and  (11)  to 
provide  us  with  the  Joint  RDF  of  tg  and  8 

,  •  ?l'j  * 

Vs  (,s.a>  •  y  *  (15) 

PUF  for  t  (8  Rayleigh) 

-  •  '!J 

f.  («_)  •  1  ,  /”  «  1  J  1  *>  (16) 
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and  CDF  for  t. 


ij*7  ,  ’ll] 

r,  !•,)  •  I  -  /“  Sr  *  0!  *  vs  (17) 

a  o  o 

eqs.  (16)  and  (17)  are  compared  graphically  with 
an  exponential  behavior  where  t(tg)  •  .58  sec. 
and  are  depicted  In  Fig  7. 


Cq.  (I/)  (or  the  COF  «r  t.  assuming  ■  at  Ray* 
Idlgh  1l  Mil  apprqxlaaated  by 
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where 


end 


(16) 


U  2/3 
T  •  I  (?) 


Cj  •  .1,  Cj  ■  2.4  and  C„  •  3.18 


t  attuned  it  Rayleigh  (0  •  .28) 

k*  (i)  •  j  j-,  /  /  .He  ( frm  . 

.  |  at  t. 

•*P  <-  J  •  (*))  dt,d* 

Eqt.  (21a)  and  (21b)  for  *  2  ara  depleted  la 
Ft 9.  3 
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In  general  the  bit  error  rata  li  equal  to  a 
function  «(.)  of  the  SNR: 

BE*  ■  y  (SNR  *  W/R) 

where,  W  It  the  bandwidth  of  the  receiver  and  R 
I*  the  data  rate.  W/R  •  1.25  at  uted  In  the 
COMET  system.  SNR  It  defined  by  equation  (1). 

In  our  dltcuttlon  two  cates  are  contldered 
blnar/  FSK  and  coherent  BPSK. 

For  TSK  we  have 


aca  (t) 


Ml.-O/l 

1  ) 


(HI 


and  for  BPSK 

ICR  (t)  •  J  «rf«  (>(l.M  »R) 

-  |  trfc  (>/!.»  •  u  exp  (((,-0/(2  1))  (20) 

Multlpljrlny  eqt.  (19)  and  (20)  by  eqt.  (12)  and 
(IS)  and  Inteyrattny  over  t»  and  8  yield  the 
averaye  bit  error  rate  for  FSK  and  BPSK. 

For  FSK:  B  attuned  Is  exponential  (B|_  ■  .35) 


Eqt.  (22a )  and  (22b)  for  •  2  ara  daplcted  In 

fig.  4 


RESULTS 

Eqt.  (13)  and  (16)  for  the  probabtlty  densities 
of  t  n  under  the  attuaptlont  that  B  It 
exponential  and  Rayleigh  respectively  provide  os 
with  the  quantitative  Infornatlon  to  evaluate 
events  of  burtt  duration  accurately.  Perhaps 
wore  Important  their  depiction  suggests  a  good 
aggreenent  with  an  exponential  behaviour  for  tg 
-  the  currently  prevailing  dogma  •  particularly 
for  the  ease  when  B  Is  Rayleigh.  The  case  for 
which  B  Is  assumed  Rayleigh  seems  correct  for  It 
yields  lero  probability  for  small  B  •  fast 
decays  -  which  aggrees  with  our  physical 
Intuition. 

The  results  for  bit  error  rate  and  average  bit 
error  rate  are  graphed  In  Figs.  3  and  4  and  the 
expressions  In  eqs.  (21)  and  (22)  can  be 
modified  to  test  the  effects  of  other  modulating 
schemes  and  effects  of  different  average  decay 
time  constants  and  minimum  signal  to  noise 
ratios. 

It  Is  Important  to  note,  however,  that  assuming 
the  same  fixed  xMifmum  allowable  bit  error  rate 
for  both  modulating  schemes,  FSK  and  PSK, 
Implies  a  higher  minimum  SNR.n,  for  FSK  then 
PSK  and  a  shorter  burst  durltlon  for  FSK 
resulting  In  a  worse  throughput  perforMnce  for 
FSK. 


(1,-1,  , 

Kl  (0*  j  //  j  «»P  (-(•  J  e  f  *  jJ  *  41,4* 

•*l .ts  (21a) 

B  assumed  la  Rayleigh  (o  *  .?R1 

•(>  (()•  j  //  tip  (-  |  (ana*  1  ^*a  y  •  ?j))  dt,dt 

••1.2$  (21b) 

Fir  coherent  BPSK:  B  assumed  Is  exponential 

.  ,  .  _ _  11.-0/21 

hr  in  •  J  j,j  /  /  j  »rtc  i/nnr, »  1  )  • 

•*p  <-‘n  *  (22a) 


In  addition,  the  graph  for  BPSK  (the  modulating 
scheme  used  In  the  COMET  System)  shows  lower 
average  bit  error  than  FSK  as  expected.  The 
lower  average  BCR  encountered  when  B  assumed 
Rayleigh  as  opposed  to  B  assumed  exponential  can 
be  explained  by  noting  the  Bz  In  the  exponent  of. 
eqs.  21b  and  22b. 

CONCLUSIONS 


A  mode)  has  been  described  for  a  meteor  burst 
communication  channel  that  takes  Into  effect  the 
Statistical  variation  In  the  decay  rate  of  the 
stgnal-to-nolse  ratio  of  the  channel.  Both  the 
PDF  for  t,,  the  burst  duration,  and  Joint  PDF 
of  tg  and  B,  decay  constant,  were  derived  and 
are  consistent  with  empirical  observations  made 
by  RADC  (6).  The  bit  error  rate  for  two 


modulating  schemes  FSK  and  .  .*>  were  calculated 
by  avenging  over  the  effect  of  decay  rate  and 
little  difference  Is  found  whether  the  Initial 
assumption  for  the  distribution  of  the  decay 
time  constant  was  Rayleigh  or  exponential.  The 
expressions  found  for  average  BER  can  be 
modified  to  test  for  other  coding  techniques. 
The  high  BER  suggests  the  need  for  further  study 
of  the  channel  and  the  use  of  coding  for 
Improving  the  BER. 
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ABSTRACT 

In  this  papar  wa  prasant  the  results  of  an 
analysis  of  a  meteor-burst  channel  showing  the 
expected  throughput  and  tha  probability  of 
coapletlng  a  given  sassage  under  two  different 
scenarios.  In  the  first  scenario,  the 
transmitter  Is  presuaed  to  have  knowledge 
(obtained  by  appropriate  probe  signals)  of  the 
entire  duration  of  a  given  meteor-burst  (HB). 

In  the  second,  only  tha  starting  time  of  each  MB 
Is  known.  Numerical  results  based  on  the 
analysis  are  provided,  and  it  Is  concluded  that, 
through  the  use  of  appropriate  signal  design  and 
error-correction  coding,  sufficiently  reliable 
communications  over  tha  MB  channel  is 
achievable. 

I.  INTRODUCTION 

Tha  possibility  of  using  Ionization  trails 
created  by  meteors  entering  the  earth's 
atmosphere  to  provide  beyond  llne-of-slght 
communications  has  been  known  for  a  long  time. 
Several  experimental  [1],  (2)  and  practical  [3J 
systems  have  been  built  using  this  node  of 
propagation.  The  physics  of  the  phenomena 
generating  these  Ionization  trails  and  the 
resultant  reflection  or  radio  waves  has  been 
described  In  detail  elserfiere  in  tha  literature 

[«].  (53. 

In  this  paper,  we  examine  the  performance 
of  two  protocols.  A  fixed  length  packet  is 
assumed.  The  first  protocol  attempts  to 
optimally  use  all  available  bursts.  The 
terminal  to  receive  data  is  assumed  to  be 
continually  broadcasting  ao  as  to  probe  for  a 
channel  opening.  The  terminal  sending  the 
message  data  begins  transmitting  as  soon  as  It 
hears  the  probing  signal.  As  a  result,  the  time 
delay  between  the  opening  of  the  channel  and  the 
start  of  the  data  transmission  Is  at  most  equal 
to  the  one-way  propagation  time.  Once  the 
channel  closes,  the  probing  signal  disappears 
and  data  transmission  ceases;  the  search  for 
another  channel  opening  now  begins. 


The  second  protocol  examined  utilizes  the 
trail  duration  less  efficiently.  However,  It  la 
simpler  to  implement  and  It  reduces  the  exposure 
of  terminals  to  detection  by  decreasing  the 
transmission  requirements.  In  this  protocol, 
the  message  sending  terminal  probes  for  a 
channel.  Channel  openings  are  detected  when  the 
sending  terminal  receives  a  response  to  Its 
probe.  This  results  in  a  delay  almost  equal  to 
the  two-way  propagation  time  prior  to  the  start 
of  data  transmission.  However,  for  typical 
trail  durations  (.2  to  1  second),  the  two-way 
propagation  time,  which  Is  usually  less  than  ID 
msec,  can  be  considered  negligible  in  most 
circumstances.  Using  this  protocol,  only  one 
packet  is  transmitted  per  burst. 

We  assume  pscketlzed  data  transmissions 
using  binary  FSK.  Noncoherent  matched  filter 
detection  Is  used,  and  s  rats  -1/2  Read-Solomon 
cods  Is  saployed.  A  stellar  analysis  of  tha 
expected  performance  of  a  communication  system 
attempting  to  utilize  meteor  bursts  for  packet 
communication  was  presented  In  (5).  However, 
aoat  of  our  models  differ  from  that  of  [5],  and 
we  are  able  to  provide  closed-form  analytical 
solutions  for  three  of  our  four  key 
expressions.  These  are  the  throughput 
expressions  for  each  of  two  modes  of  operation 
of  the  meteor-burst  channel,  as  well  as  the 
probability  of  successfully  completing  a  packet 
for  one  of  the  two  aodes.  The  probability  of 
successful  completion  for  the  other  mode  of 
operation  is  presented  as  an  infinite  suxnatlon, 
similar  to  an  expression  given  In  [5]. 

II.  ANALYSIS 

In  this  section,  expressions  for  both 
throughput  snd  probability  of  eorrsctly 
transmitting  a  g'ven  messagt  art  presented.  The 
supporting  analysis  leading  up  to  thass  rssulta 
can  be  found  In  (6).  As  Indicated  above,  two 
different  modes  of  system  operation  are 
considered.  In  the  first  mods,  It  la  assumed 
that  the  receiver  continuously  transmits  a 
sounder  to  the  transmitter.  Any  time  the 
transmitter  rscsivss  the  sounder  etgnel,  it 
knows  that  the  meteor  burst  (MB)  channel  Is 
present  and  hence  starts  to  transmit  its 
message.  As  soon  as  reception  of  the  sounder 
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signal  oeaeea,  the  transmitter  stops 
transmitting.  The  transmission  of  a  given  packet 
•  might,  howwver,  apan  many  MB's,  and  it  la 
aaauaad  tltat  a  oar  tain  aaownt  of  overtimed  la 
neoeaeary  for  tha  transmitter  to  aoquire  tha  HS 
ohannal  on  aaeh  and  ovary  aaparata  burst, 
apeoiriealiy ,  it  la  aaauaad  that  t.  aaoonda  ara 
naeaaaarjr  for  aaoh  acquisition.  U  aaaanoa 
than,  tha  firat  aoda  of  oparatlon  corresponds  to 
tha  tranaaittar  knowing  praoiaaly  aaoh  inatant 
or  tlaa  whan  uio  MB  ohannal  la  praaant. 

With  thla  acanarlo,  tha  following 
mathematical  nodal  la  uaad't  Lot  n  Pa  tha  total 
number  of  floor  burat  a  that  ooour  in  aof  tlaa 
Interval  I  aaoonda,  and  lot  B  Pa  tha  nunbar 
of  burata  In  T_  aaoonda  that  allow  axaotljr  J 
pack at a -par -burat  to  ba  trananlttad.  Than 
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l  P(«j |n)P(n) 
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whara  n  la  taken  to  ba  a  Poisson  randoa  variable 
with  probability  mass  function  glaan  by 

h 

'  T.  T, 


X(?J 
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In  (2),  T.  Is  tha  averege  lntarval  betwaan 
bursta.  tlaarly ,  thla  aodel  la  only 
approxlmata,  slnca  the  suuatlon  in  (I)  allows 
an  infinite  nuabor  of  bursts  to  occur  in  a 
finite  Interval  of  tlaa.  However,  if  tha 
average  Interval  between  burata  Is  auch  lass 
than  tha  tlaa  duration  T  .  the  approximation 
should  yield  reasonable  results. 

If  we  define  P  as  the  probability  that  a 
given  burst  lasts  for  a  duration  equal  to  that 
of  exactly  J  paokats,  then 
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To  eveluate  P(H  ),  the  distribution  of  the 
duration  of  a  burst  oust  be  known.  Towards  that 
end,  assume  the  probability  density  of  a  burst 
duration  is  given  by 
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where  T  la  tha  avarage  length  of  a  burst, 
this  model,  it  Is  straightforward  to  ehow 
that  Nj  la  also  Poisson.  That  la. 
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With 
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To  as  great  an  extent  as  Is  possible,  we 
use  tne  foraulatlon  and  notation  of  £ 5 J . 


and  where  T.  la  tha  duration  of  a  packet. 

TO  obtain  tha  axpaotad  nuabar  of 
Information  bits  that  oan  ba  transmitted  over 
tha  HB  channel  In  T  aaoonda,  whloh  wa  refer  to 
as  the  throughput2,  let  t.  be  defined  as 

B_  *  I  J  W,.  IT) 
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where  I  la  the  nuabar  of  information  bits  par 
packet.  That  Is,  kB^  Is  tha  actual  nuabar  of 

transaltted  bits  In  one  realisation  of  tha 
process.  Then 
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Consider  now  the  second  protocol,  whereby 
the  transnltter  makes  no  attempt  at  monitoring 
tha  presence  of  the  channel  other  than  the  start 
of  the  burst.  However,  each  tine  a  burst  la 
detected,  an  entire  packet  la  transaltted. 

Then,  with  I  again  equal  to  the  number  of 
information  bits  per  packet,  the  number  of  bits 
transmitted  across  the  channel  in  a  time 
duration  of  TQ  seconds  is  given  by 

"B  •  I  Bp,  (9) 

where  Np  Is  defined  as  the  number  of  bursts 
whose  duration  is  greater  than  or  equal 
to  t_ *T _ .  Notice  that  the  parameter  t.  Is 
larger  In  this  second  mode  of  operation  than  It 
was  In  the  first  mode. 

With  this  model,  the  expected  number  of 
bits  that  can  be  transaltted  over  the  channel 
can  be  shown  to  be  given  by 
I  T 

E  «*P  |-(Tp  ♦  V'Tj  •  l,0) 


While  the  throughput  expressions  derived 
above  give  some  perspective  regarding  the 
performance  or  a  MB  channel,  at  times  a  more 
meaningful  performance  measure  Is  the 
probability  of  successfully  completing  a  message 
within  a  specified  period  of  time.  Considering 
first  the  mode  of  operation  wherein  the  reoalver 


1  Note  that  this  definition  ignores  the 
effect  of  thermal  noise.  The  throughput 
expressions  given  here  can  be  modified  as  In  £9] 
to  Include  errors  due  to  noise  by  assuming  a 
specific  IRQ  model.  However,  noise  effects  are 
included  In  this  paper  rfien  the  probability  of 
correct  transmission  Is  considered. 


continuously  sounds  ths  ohannal,  suppose  a  (Ivan 
packat  la  repeated  for  aa  ssny  tlaaa  aa  poaalbla 
over  a  given  duration  of  tlae,  aajr  T_  aaoonda. 
Lot  L.  bo  daflnad  aa  the  grestest  Integer  laaa 
than  4r  aqual  to  T-/T-,  and  lot  P_M  b«  tha 
probability  of  oorr tolly  ooapleting  tha  message^ 
on  at  laaat  on#  of  L.  successive  tranaalaalona 
of  tha  sea#  packat. 

If  fp  la  daflnad  aa  tha  probability  of 
aaklng  an  error  In  a  apaolflo  packat,  given  that 
It  haa  baan  received  (l.a.,  given  that  tha 
ohannal  waa  praaant  long  anough  for  tha  packat 
to  ba  ooaplatad),  than  It  can  ba  shown  that 
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Note  that  In  deriving  (11),  it  waa  assumed  that 
arrora  In  all  packata  occur  with  the  aaaa 
probability.  This  la  not  true  In  reality,  since 
the  strength  of  a  MB  decays  approximately 
exponentially.  However,  If  we  uae  the  smallest 
algnal-to-nolae  (SNR)  In  any  packet  for 
coaputlng  the  probability  of  error  or  all  the 
bits  In  tha  packata,  than  (11)  will  be  a  worst- 
case  result. 


Finally,  for  tha  second  protocol,  It  can  be 
shorn  that  _ 
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The  results  derived  up  to  this  point  are  In 
teras  of  tha  parameter  Pp,  tha  probability  of 
aaklng  an  error  In  the  picket.  Hence  an 
explicit  expression  for  P.  is  needed  In  order  to 
evaluate  the  overall  system  performance. 

Towards  this  end,  asauae  that  binary  FSK  with 
noncoherent  detection  la  tha  modulation  foraat 
and  that  the  overall  packet  Is  to  be  broken  up 
Into  smaller  groups  or  Information  digits  which 
are  to  be  encoded  with  a  forward-error 
correction  (FEC)  coda  and  then  transmitted  over 
tha  MB  ohannal. 

As  indicated  above,  because  the  SNR  of  a 
given  MB  decays  approximately  exponentially  In 
tlaa,  the  probability  of  error  due  to  thermal 
noise  at  the  receiver  will  vary  from  syabol-to- 
syabol.  Therefore,  to  keep  the  analysis 
tractable,  the  probability  of  error  results  will 
bo  ooeputed  using  tha  SNR  that  Is  present  at  the 
end  of  a  MB  (l.a.,  using  the  saallest  SNR  of  the 
burst),  so  that  tha  final  results  will  really  be 


1  For  the  purposes  of  this  section,  a 
message  will  ba  taken  to  be  a  single  packet. 


upper  bounds.  Also,  again  because  of  tha 
exponential  decay  In  SKR,  It  is  reasonable  to 
assuae  that  tha  ohannal  will  be  quits  bursty  in 
nature.  Hence  the  FEC  a  chase  should  bs  one  that 
performs  well  In  tha  presence  or  bursty  errors, 
and  for  the  purposes  of  this  paper,  a  Reed- 
Soloson  (RS)  code  will  be  used. 

To  be  speclfle,  assuae  the  paokst  oonsleta 
of  I  bits  of  information,  and  these  I  blta  are 
broken  up  into  J  groups  of  aK  bits  par  groie. 

An  (N,K)  RS  code  la  then  used  to  enoode  each' 
group,  where  ■  -  2*  -  1 .  To  reoelvs  the  packet 

correctly,  all  J  oode words  must  be  correctly 
received.  The  probability  of  thla  occurring  la 

.  03) 

where 

pw  ■  *  <’  *  yN~*.  o*> 

I  -E*l 

E  Is  the  error-correction  capability  of  the  RS 
code  and  P  is  the  error  probability  of  a  IS 
symbol.  It  Is  given  by 


s 

where 


(1 


1 


V 


p  .  1  I 

rb  2  • 

is  the  bit  error  rate  Tor  noncoherent  binary 
FSK. 


(15) 

(16) 


III.  NUMERICAL  RESULTS 


Figures  1  and  2  contain  curves  showing  the 
probability  of  completing  a  message  for  tha 
first  and  second  protocols,  respactlvaly  (l.a., 
equation  (11)  is  plotted  in  Figure  1  and 
equation  (12)  is  plotted  In  Figure  2).  Tha 
abscissa  in  both  cases  Is  TD,  tha  total 
observation  time  during  which  the  aessage  Is  to 
be  completed.  The  curves  in  each  figure  are 
parameterized  by  the  si gnal-to-nolse  ratio, 
which  in  turn  Implies  a  value  for  T,,  the 
average  time  between  bursts.  For  the  results 
presented  In  Figures  1  and  2,  the  following 
pairs  of  SNR  and  Tj  were  used: 


SNR(dB) 

T 

7  .5 
8 

10 

12 


T.  (seconds) 
28 
29 
3' 

39 

50 


The  message  duration  (packet  size)  was 
taken  as  2(00  Information  bits,  and  those  2*00 
bits  were  divided  Into  32  codewords.  Each 
codeword  corresponded  to  a  (3'.'5)RS  code. 
Finally,  TL  and  Tp  were  each  taken  to  be  0.2 
seconds,  and  tg  was  taken  to  be  .03  seoonds  for 
the  first  protocol  end  .06  seconds  for  ths 
second  protocol. 

Pecause  the  value  of  Tj  Increases  ae  the 
required  SNR  Increases,  at  some  point  requiring 
a  larger  snr  becomes  self-defeating,  since  the 


nusbir  of  aeteor-bursta  in  a  gtvan  Tp  second 
interval  beoo—a  too  aaall  to  aalntaln  a 
specified  degree  of  reliability.  Thi»  la 
ol early  indicated  in  oaed  of  tna  two  figures, 
where  it  la  aoon  that,  for  tha  paraaeter  valuta 
onoaan,  an  Mil  of  •  SB  provldvs  uia  boat 
roaulta. 

If,  than,  an  6  dS  MR  threshold  la  uaad, 
and  If,  for  eiaaplt,  it  la  doalrad  to  have  a 
probability  of  aeasage  coaplatlon  of  0.9,  froo 
Figure  1  It  la  aaan  that  tha  observation  tloa 
■uat  bo  about  throo  olnutaa  whan  tha  firat 
protocol  la  uaad.  Similarly,  if  tha  aacond 
protocol  la  ohoaan,  figure  2  lndloataa  that  the 
obaarvatlon  tloa  ahould  bo  ineraaaad  to  about 
four  and  one-half  olnutaa  in  ordar  to  achlava 
tho  aaoa  porfomanoa. 

19.  conclusions 

An  analysis  of  two  different  protocols  for 
uaa  over  an  MB  channel  has  bean  presented.  It 
was  found  that  with  sufficient  repetition  of  tha 
eassaga  and  with  a  sufficient  amount  of  error- 
correction  coding,  reliable  comuni  cat  Ions  could 
be  achieved.  This  Indicates  that,  at  least  In  a 
back-up  node,  the  MB  channel  la  quite 
appropriate  for  tha  transolsslon  or  relatively 
short  oesaages.  In  fact.  It  ahould  be 
eophaslted  that  the  results  presented  here  are 
pesalolatlc  results.  In  the  sense  that  they  were 


derived  under  tha  asauaptlono  that  tha 

Instantaneous  SNR  on  tha  ohannel  waa  at  oil 

tines  equal  to  the  lowest  MR  that  waa  d eased 

scoop table. 
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Figure  8  -  Observation  time  TD,  seconds 
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ABSTRACT 

In  this  paper  ve  analyze  an  accurate  channel 
nodal  baaed  on  the  classical  underdense 
aquations  and  recent  enplrical  data.  The 
analysla  results  In  analytical  expressions 
for  channel  duration  and  throughput  to  be 
used  as  design  and  analysis  tools  in  constant 
transmission  rate  communication  systems.  In 
addition,  ve  derive  the  optimum  constant 
transmission  rata  for  a  given  communication 
system  where  all  bits  are  constrained  not  to  I 
exceed  a  given  maximum  bit  error  rata.  This 
optimum  bit  rata  yields  the  maximum  average 
throughput. 

I  IHTRODOCTIom  ■= 

For  a  communications  system  engineer  the 
statistical  behavior  of  meteor  burst  channel 
duration  is  crucial  in  understanding  through* 
put  performance.  Ve  start  by  employing  a 
complex  stoclifcstlc  modal  for  the  channel 
where  both  the  underdense  electron  line 
density  and  the  decay  time -constant  are 
random.  Clearly,  for  a  given  set  of 
communication  system  parameters,  noise  and 
maximum  allowed  bit  error  rate  (BEK) ,  choice 
of  bit  rate  is  equivalent  to  s: reifying 
received  signal  power  threshold  abi  a  which 
ve  operate.  The  expressions  for  channel 
duration  and  ita  statistics  as  defined  by  the 
maximum  allowed  BEK  are  derived  in  terns  of 
the  system  (fixed)  transmission  race.  Ve  then 
derive  the  expression  for  throughput  in  terms 
of  bit  rate  and  average  it  over  the  ensemble 
of  underdense  trails.  The  optimal  bit  rata 
that  maximizes  throughput  is  then  found; 
demonstrating  the  .improvement  over  existing 
sub-optimal  systems. 
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2  MATHEMATICAL  MODEL  AND  CONCEPTUAL  FRAH EVOR 

In  this  section  the  basic  equations  tha 
model  the  MB  channel  are  presented.  1 
addition,  assumptions,  constraints  and  emplr 
ical  factors  are  atatad. 

2.1  Tower  Received  and  Related  Parameters 

The  power  received  from  the  trail  model  use 
hers  la  from  the  classic  papers  by  Eahelma 
and  Manning  [2]  and  Sugar  [1]. 

2.1.1  ieuBr.Hts.siv.B4  lira  »t>  SnflRiAaai*  1r«l 

Meteor  trails  with  fewer  then  1014  electron 
per  meter  of  length  are  referred  to  a 
underdense  and  the  equation  governing  th 
power  received  from  the  trails  is  given  by: 

f(t).  ^C,C.r?s.B»(«)A»«,(-^i)  ,  ^*'0  y 

16*’ArA ,( A,  ♦  A *)( ) -ces’/sln’#)*  v  i’set*#  j 

where.  (2.1.1 -I) 

Ft  Transmitter  power. 

Gy.Gg  Transmitter  and  receiver  antenna  gains, 
r,  Classical  radius  of  an  electron, 
e  Angle  between  the  electron  field 

vector  E  and  R,. 

q  Electron  line  density  of  the  trail, 
i  Carrier  wavelength, 

r*  Initial  radius  of  the  trail  (.65  m.). 

e  Angle  of  lncldence/reflectlon  of  the 

transmitted  plane  wave. 

Xj  Distance  from  transmitter  to  trail. 

Kg  Distance  from  receiver  to  trail. 

D  Diffusion  coefficient  of  the  atmo¬ 
sphere  .  V 

A  Angle  between  the  trail  and  the 

propagation  plane  formed  by  Rf  and  Kt. 

A  minimal  set  of  assumptions  is  applied  lr. 
the  literature  [3,4,5]  for  analysla: 


1. )  Th#  trail  occur*  at  Midpoint  between 
racalvar  and  transmit  tar  l.a.  R^-R, 

2. )  Th*  burst  occurs  at  an  altitude  of  100 

km. 

3. )  Tbs  trail  Is  travailing  at  a  plan*  . 
parpandlcular  to  tbs  plana  formed  by  R,  and  j 
Ry  l.a.  »*t/J  .  In  addition,  v*  sat 

Under  the**  assumptions  and  after  substitut¬ 
ing  for  th#  physical  constants  vs  rearrange 
th*  above  aquation  to  arrive  at  the  power 
received  equation  for  the  underdens*  case. 

/•(O-CeflV*"  (2.1.1 -2)  • 


where, 


C„  •  2.S179581 • 10 


^TCrC,iasxp(-gg^) 


3.166-10 -.Olgli 


It  should  be  emphaaized  that: 

1. )  P(0),  the  naxlnuM  power  froa  th* 
underdens*  trail  Is  proportional  to  q*  ;  for 
a  specific  trail  q  la  fixed. 

2. )  Cy,  th*  constant  of  proportionality, 
Incorporates  the  effects  of  link  geoaetry, 
transmitter  power,  antenna  gains  and  carrier 
wavelength.  For  a  given  communication  systea 
Cy  Is  constant. 

3. )  Th*  decay  tin*  constant  B  given  above  Is 
found  to  be  randomly  varying  froa  trail  to 
trail.  Th*  expression  in  th*  classic  aodel 
corresponding  to  a  decay  tla*  constant  Is 
taken  as  th*  average  decay  tlae  constant  in 
accordance  with  observations.  To  date, 
however,  researchers  have  assumed  all 
underdense  trails  to  share  the  same  decay 
tine  constant  to  fascllltat*  their  analysis. 
It  Is  in  this  added  complexity  to  the  modal 
that  our  treatment  deviates  from  published 
analysis  and  enhances  its  accuracy. 

4. )  The  decay  time  constant  B  and  th* 
electron  line  density  q  are  and  statistically 
Independent  { 1 ]  ( 3 ] . 

3.)  For  a  given  comaunlcatlon  system  (i  e.  C, 
fixed)  knowledge  of  the  electron  line  density 
,q,  and  decay  tla*  constant  ,B,  completely 
specify  th*  underdense  trail  behavior  lr. 
tla*. 


Based  on  recent  empirical  data  [6]  the  pc 
for  q  Is: 

/,(<?)-$<?*'  :  10u(2. 1.2-1 


where, 

q  -(p-i 


-fer 


p  “  1.6 


q_,.  Is  th*  ainlanJM  electron  line  density  q 
trail  xust  posses  to  be  ’seen'  by  e 
communication  systea  l.a.  Cy  g*..,.  Is  t 
minimum  power  level  ?_,.  detected  by  t 
communication  systea.  In  the  context 
constant  transmission  rate  systems  ( 
disscussed  later  )  P.,.  Is  the  power  lev 
that  corresponds  to  th#  maximum  allowed  b 
error  rat#  (for  a  given  bit  rate  and  noise) 


V#  start  by  stating  that  the  expected  val 
of  th*  decay  tlae  constant  B  Is  we 
approximated  by: 

£(B)eg-3.l6610'J- (2.1.3- 


(2.1.3- 


The  random  decay  tlae  constant  B  is  ass  us 
to  be  Rayleigh  or  exponential  which  has  be 
shown  to  be  consistent  with  recent  experla< 
tal  observation  |6). 

XayieigJi; 


/.(b)--V 


(2.1.3- 


Slnce  the  electron  line  density  and  the  dec 
time  constant  are  statistically  lndepende 
the  joint  pdf  Is  given  by  their  product. 


In  this  section  we  present  the  reletlonsh 
between  link  distance,  L,  and  station 
trail  distanca,  Rj  (-Ry)  and  angle 
Incidence/reflection  ♦ 


RT’Rlt->j(h*li/BR,y*L7/4 


(2.2- 


secJ(#)-  1  * 


(2.2-: 


where. 

h  trail  altitude  (100  km) , 
t  great  elrcla  distance  between  stations. 
Kg  diatanea  froa  transmitter  to  trail. 

Kg  diatanea  froa  racaivar  te  trail. 

Jt,  radlua  ot  tha  aartfa,  6400  km. 


2.)  Trail  Arrival 


tfa  lat  A'  (q)  b*  tha  number  of  tralla/a  with 
an  aloctron  llna  density  batvaan  q  and  q+dq 
a/a.  Using  Ualtzan  (71  we  write 

9 


where  q_..<o<q _ -10n  a/a  and  {  la  tha 

ldaallsad  baaawldth  of  tha  antenna  and  all 
othar  parameters  have  baan  previously 
defined.  Tha  number  of  tralla/a  with  q  In  a 
given  range  (qj.q,].  A,  la  found  by 
Integrating  tha  above  density  over  q 
yielding  the  following 


(2.3~  2) 

«J|  L  flaj 

*•*  Imn  SMCUtl  Penalty  of  Received  balsa 

For  signal  frequencies  appropriate  to  the  MAC 
the  received  noise  power  spectral  density, 
H,.  Is  a  function  of  the  galactic  noise 
picked  up  by  the  receiver  antenna  and  tha 
receiver  thenaal  noise.  Using  Abel  (4)  ve 
have: 


(2.4-1) 


where,  Boltzmann  constant,  1 ,3805  •  lO'nJ/K 


2"«"  290  K  F  —  receiver  noise  figure 

ig  -  power  loss  between  the  antenna  and  receiver 

1«  wavelength  in  meters 
In  our  analysis  we  shall  use  L,  and  F  as  1.3 
and  2.3,  respectively . 

2.5  Modulation  and  Bit  Error  kata 

Tha  BER,  denoted  here  as  Pk,  Is  related  to 
tha  bit  power  P(t) ,  transmission  rata  R  and 
power  spectral  density  of  the  received  noise 

H,  as 

where  g(  )  la  determined  by  tha  modulation 


technique . 

ot  altemtlvely,  wa  can  aolva  for  tha  power 
^(O-  Ng- R  ■ff*'(/,a(0)  (2.S-2) 

where  g‘l  la  the  Inverse  of  g. 

3  ANALYSIS  AND  RESULTS 

3.1  Introduction  and  Approach 

Our  goal  la  to  estimate  the  optimal  average 
throughput  for  a  given  communication  system 
a  system  with  a  constant  transmission  rati 
and  a  time  varying  bit  error  rate  which  1: 
constrained  not  to  exceed  some  maxlmui 
allowable  value.  For  a  constant  bit  rat< 
system,  since  the  power  la  monotonlcall; 
decreasing  In  time,  the  probability  of  bl 
error  increases  with  time  In  a  fashlo 
dictated  by  the  relevant  modulation  function 
Ha  sea  the  need  then  to  impose  a  celling  o 
tha  BCR.  Ha  thus  specify  a  maxiaua  allovabl 

BER,  p, _ _  such  that  all  tranmltted  bit 

shall  have  a  BER  less  than  or  equal  to  P, _ 

This  constraint  on  the  probability  of  bl 
error  for  constant  transmission  rate  system 
Implies  tha  existence  of  a  minimum  powei 
level  P.,.  such  that  all  transmitted  bit) 
have  a  received  power  greater  than  or  equa: 

to  F^,.  In  short.  Pfc(t)  -  P, _ when  P(t)  ■ 

P— -■  Tha  time  for  which  tha  power  recelvai 
froa  the  meteor  burst  exceeds  a  prescribe 
threshold  (P...)  Is  defined  as  the  channe 
(burst)  duration  t,.  This  random  burs 
duration  Is  crucial  In  determining  the  numbe 
of  bits  transmitted  during  a  given  burst  an 
ultimately  In  estimating  throughput.  If  w 
choose  to  transmit  at  a  higher  rate  th 
corresponding  has  to  be  Increase, 

resulting  In  shortening  the  effectlv 
duration  of  transmission.  This  clear  trade 
off  between  transmission  rate  and  duration  o 
transmission,  t»,  suggests  an  optimal  cholc 
for  transmission  rate  (or  equivalently  a 
optimal  t()  so  as  to  maximize  the  throughput 
The  link  between  the  number  of  bit 

transmitted  per  burst  and  throughput  la  dor 
through  the  average  arrival  rate  statistic 
and  trail  shape  statistics  (variation  In 
and/or  B).  The  Judicious  choice  of  transmls 
slon  rate  for  che  given  system  provides  us 
therefore,  with  the  optimal  throughput. 


p  -  1. 


3.1.1  fs.gtat.gl 


Both  station*  hava  receiving  and  transmitting 
capabilities.  Th*  station  which  acts  as 
transmitter  of  data  listens  continuously  for 
a  continuous  ton*  (or  a  frequently  transmit¬ 
ted  probe)  vhlch  Is  being  sent  from  the 
receiver.  Upon  detection  of  this  tone  both 
th*  presence  and  th*  strength  of  the  channel 
Is  known  to  the  transmitter  which  Immediately 
commences  transmission.  Th*  situation  Is 
completely  symmetricel  with  respect  to  th* 
receiver. 


3.2  Improvements 

The  adventsges  of  the  analysis  herein  stem 
basically  from  considering  a  complex  stochas¬ 
tic  KBC  model  vher*  both  the  initial  power 
and  decay  time  constant  are  considered 
random.  The  expressions  derived  will  enable  a 
system  designer  to  estimate  th*  best 
throughput  by  employing  the  optimal  transmis¬ 
sion  rat*. 

3.3  5«roL»  iaai 

For  quantitative  appreciation  of  the  results 
we  assume  some  practical  values  for  a 
communication  system. 

L  Varies.  Sets  R»  and  «(Sec.  2.2) 

Pt  1000  w,  Gj— 10  d&,  D«  8  m V* 

{-«.*<•;  l  varies;  SB  from  Sec.  2.4 
P, _ -10-»  Using  BPSK  w*  get  s-1(P, _ 1-9 

3.4  sonmns- Bit,  totv.  SanmtlngA  BLima 
^3.4.1  ghtnael  JEamign 

Channel  duration  or  burst  duration,  t.  is  th* 
-time  for  which  th*  power  received  from  th*  ! 
burst  exceeds  a  given  minimum  value  P.l0  . 

3. 4. 1.1  Burst  Duration  Statistics 

Assuming  a  constant  transmission  rate,  R, 
some  given  maximum  BER,  and  power  spectral 
density,  N,  w*  get  from  Eq.  (2. 1.1-1). 


(3. 4. 1.1-  1) 


where 


<-v 


and 


C„  -2.S179581  10 


-w  ^rCrCsi*«xp(-~~) 


Rr 


or  equivalently  a*  a  function  of  R 


t,  ■*  B  •  In  j 

'  Co 

'(^e - 

)q 

tfe  now 

average 

over 

the 

underdens*  bursts  (avereglng  over  q  and  b)  to 
yield  th*  average  burst  duration: 

The.  Aver  a  a*  Burst  Duration--  Underderdensa : 

C3-41-1'3) 

where  ti! 

T_/rq„.y1  g-y^n  ■ 

V  Q O  )  C yOp  J 


g-3.16610'*  — -c-  t  „  .l0i4 

0 

— )*:  R - Cl/<?* 

For  a  given  communication  system  C,.  R,  H 
F bu>  *  *r*  specified,  t,  l*  tt 

normalised  (  with  respect  to  th*  hlghei 
possible  ret*  for  an  underdens*  burst)  bj 
rata. 

Fig.  1.  plots  the  above  for  our  sample  ayst< 
as  a  function  of  R,. 

Dll  &u  moment  si  £»  il  riven  by: 

(p-1)*  1-*  \  &  *  i'.  J 

(3. 4. 1.1-4) 

Hu  I££  (md  ££f  fii  Si  are: 

(3.4.1.1-S 

1  -  *F,(bmin)- j^exp^^-ii  )/,(6)dt 


where  x  and  p  are  from  above.  p‘-(p  -l)/2 
b*u  “VI  WP’>  l«U/x)  ):  <?'  -l/(l-x) 

Fig.  2  depicts  th*  pdf  for  B  assume 
Rayleigh,  not*  th*  close  aggreement  with  tl 
exponentially  distributed  t,  (dashed  curve 
th*  current  dogma  [ 8 ] . 

3.4.2  Throughout  for  Constant  Bit  Rate 

Let  N’V  be  the  number  of  bits  transmitted  c 
the  1“  burst  and  u.  as  the  number  of  burst 
occurring  in  the  period  of  r  seconds.  Sine 
the  bursts  are  independent  of  each  ocher,  a- 
assuming  no  overlap  of  bursts,  we  have: 

1  p. 

T  ”  lim,..- A/,,  -  A  ■  N,  (3.4.2-; 

T  i- 1 

where  A  is  the  average  number  of  bursts/s  a 
V,  is  the  expected  number  of  bits  per  bur 
as  averaged  over  the  ensemble  of  all  burst 
profiles  l.e.  over  q  and  B. 

3. 4. 2.1  Optimal  Average  Throughput 

The  maximum  throughput  is  found  by  firs 
deriving  the  dependence  of  the  avers; 
throughput  (Eq.  3. 4. 2-1)  on  the  bit  rate  1 
and  solving  for  the  optimal  bit  rate  R*  th; 
would  yield  tha  desired  maximum  throughpi 
T*.  The  solution  is  done  with  respect  to  tl 
normalized  bit  rate,  Rj,  which  yields 
result  independent  of  a  particular  choice  i 
system  values . 

For  a  given  system  P,^,,  N0,  R  and  Cv  ai 

specified.  Tha  average  number  of  bits  p; 
undardenae  burst  can  be  expressed  i 

terms  of  the  normalized  bit  rate  using  Ec 


(3. 4. 1.1-5)  «•: 

(3.4.2.1-1) 

where  *11  parameters  have  been  previously 
defined.  Not*  that  in  the  above  expression 
the  bursts  considered  were  chose  with  qM1,<  q 
<qn.  The  arrival  rat*  for  such  condition  is 
given  by  Eq.  (2.3-2)  where  qt  is  set  to  q.1B 
and  q,  is.  qu-  Together  with  the  definition 
of  t£,Eq.  (3. 4. 1.1-3} .  we  have  for  the 
arrival  rate  of  underdense  bursts: 

(3. 4. 2. 1-2) 
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The  expression  for  throughput  using  under- 
dens*  bursts  is  given  by  the  product  of  the 
last  two  Eq*.  **  specified  bv  Eq.  (3. 4.2-1). 
The  result  is  a  function  of  the  normalized 
bit  rat*.  In  the  product  we  substitute  for 
r  froa  Eq.  (3. 4. 1.1-3)  and  for  and  » 
frSs  Eq.  (2. 1.1-2).  For  a  given  set  of 
systea  paraaeters  we  thus  have  the  throughput 
f or  any^cho ice  of  (normalized)  bit  rate  as: 

frC,G.t*4sln4scc»4*xp(-— -jr,)** 
T, - 3.056-  W" - (p-UW.-ff-'lf—  .):D-Tt 


1-Ai  1 

13.4.2.1-3} 

where  the  parameters  are  defined  above. 

Not*  that  the  equation  is  written  in  two 
parts;  the  first  incorporates  all  the 
parameters  of  the  particular  system  and  the 
second  part  reflects  the  variation  in 
throughput  with  choice  of  bit  race.  (Fig.  3) 

.«*  would  like  to  find  from  the  above 
•quaclon  the  optimal  noraallzed  bit  rate^Ra 
which  yields  the  maximum  throughput  T„ 

This  optimal  normalized  bit  rate  for  anv 
given  system  is’found  numerically  to  be(Fig.  3) 

or  equivalently  the  optimal  bit  rate, 

R*  -  .07  R  (3. 4. 2. 1-4)  See  Figs.  4  ,  5. 

resulting  in 


r; -3.433- 10  * - -  —,7; - 7Z - - 

(3.4.2. 1  -S) 


All  the  parameters  are  previously  defined. 

The  above  expression  for  the  maximum 
throughput  T„*  allows  us  to  quantify  the* 
affects  of  various  system  parameters  and 
constraints  on  the  best  throughput  using  a 
constant  bit  rate.  As  an  example  we  plotted 
the  variation  of  Tn*  as  a  function  of 
wavelength  and  link  distance  for  our  sample 
system  See.  figs .  8  and  7 . 


4  CONCLUSION 

We  note  from  the  result  for  burst  durati 
statistics  that  if  the  curve  for  the  pdf,  F 
2,  is  approximated  by  an  exponential,  t 
expected  value  needed  for  good  aproxlmati 
is  qlose  to  the  expected  decay  time  constan 
This  can  also  be  seen  from  Fig.  1.  (W»“  -0 
curve  2)  which  is  somewhat  more  conaervati- 
.then  previously  believed  (81.  (9j.  The  qu* 
titativ*  result*  for  channel  duration  a 
throughput  lets  us  choose  our  bit  race  a 
evaluate  the  tradeoffs  with  respect  to  oth 
communication  system  parameters  such  as  P 
P. _ etc. 
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